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PREFACE 


La théorie des équations intégrales, née d’hier, est d’ores et déja 
classique. Elle a fait son entrée dans plusieurs de nos enseignements. 
Nul doute que — peut-étre a la faveur de nouveaux perfectionnements 
— elle ne s’impose bientot 4 la pratique courante du calcul. C’est une 
fortune rare parmi les doctrines mathématiques, si souvent destinées 
a rester des objets de musée. 

Ce sort exceptionnel est, cependant, a notre avis, conforme a la 
logique. 

A mesure que, en Analyse, problemes et méthodes tendent a perdre 
leur caractere formel et a dépasser le cercle des cas d’intégrabilité 
proprement dits, il semble bien que lVintégration et non plus la 
différentiation, doive apparaitre comme l’élément simple — comme 
Voutil le plus usuel, parce que le plus puissant et le plus maniable 
— du calcul infinitésimal. L’intervention des équations intégrales dans 
létude des problemes de la Physique mathématique est, au fond, une 
phase de cette évolution. 

Il nous parait souhaitable que celle-ci ait, dés a présent, sa 
répercussion sur lenseignement. En tout cas, la belle méthode que 
Von doit a M. Fredholm, marque un tel progrés qu’il importe de la 
rendre accessible non plus seulement aux futurs docteurs et a ceux 
qui poursuivront les examens d’ordre élevé, mais a tous ceux qui, a 
quelque degré que ce soit, étudient les mathématiques supérieures. 

Une telle nécessité a été ressentie un peu partout, et, a l’étranger, 
d’excellents exposés, — tels que l’élégant traité de M. Bocher, pour 
n’en citer qu’un — ont été consacrés a la méthode qui nous occupe. 

MM. Heywood et Fréchet, en abordant a leur tour le méme sujet, 
ont visé a etre clairs, élémentaires et pratiques. Ils ont retenu de la 
théorie tout ce qui a déja acquis sa forme définitive et pratiquement 
utilisable, et se sont bornés a cet ensemble, déja singuliérement fécond 
a lui seul et suffisant dans tous les cas usuels. Ils n’ont pas séparé 
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cette théorie des applications qui en sont la raison d’étre et l’origine 
méme, et qu’ils ont passées en revue avec grand soin. Grace a l’ceuvre 
ainsi congue, nos étudiants pourront, tout en restant surs de ne pas 
étre entrainés dans des difficultés inutiles, posséder aisément le nouvel 
instrument analytique. 

Ils devront ce résultat a une collaboration internationale a laquelle 
on ne saurait trop applaudir. M. Heywood a été, il y a quelques 
années, notre hdte et l’auditeur assidu de nos cours parisiens. II s’en 
est souvenu en prétant cette fois son concours & un de nos jeunes 
mathématiciens dont le nom et le talent sont déja assez connus pour 
que nous n’ayons pas a le présenter au lecteur, en vue d’une ceuvre 
qui sera bonne et utile. C’est la une heureuse initiative : puisse-t-elle 
trouver de nombreux imitateurs! 


Jacques HADAMARD. 


L’EQUATION DE FREDHOLM 
ET SES APPLICATIONS A LA PHYSIQUE MATHEMATIQUE 


INTRODUCTION 


1. Caractere des questions traitées dans ce livre. — Les 
méthodes que nous allons développer s’appliquent surtout aux 
problemes de Physique mathématique qui concernent les équations 
aux dérivées partielles du type elliptique (voir plus loin n° 1, p. 14), 
dont le plus connu est le probleme de Dirichlet (n° 4, p. 17). 

Dans la plupart des cas déja étudiés, les problemes en question se 
ramenent a une équation de Fredholm, c’est-a-dire 4 une équation de 
la forme 


b 
(1) els) f K(s,1) o(t) dt = F(6) 
ott K(s,t)("), f(s) sont des fonctions connues. On cherche une fonction, 


y(s), qui satisfasse a cette équation. Le paramétre \ est introduit pour 
faciliter la discussion. Nous considérerons avec (1) l’équation associée 


b 
(2) b(s) — i K(t, 8) W(t) dt = f(s). 


Le Deuxiéme Chapitre sera consacré a la résolution de cette 
question, qui est effectuée par plusieurs méthodes distinctes donnant 
y(s) sous des formes différentes. 


(') On appelle cette fonction K le noyau. 
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La premiére méthode, celle d’itération, due 4 Neumann et Volterra, 
nous donne une série entiére en , dont les coefficients sont des 
fonctions de s. Elle est valable pour les valeurs de \ plus petites en 
module qu’un certain nombre fixe. 

La seconde méthode, celle de Fredholm, donne pour y(s) le rapport 
de deux fonctions entiéres en \, c’est-a-dire une fonction méromorphe 
(voir p. 5) de A. Ces fonctions sont obtenues sous la forme de séries 
entieres dont les rayons de convergence sont infinis. Le numérateur du 
rapport a pour coefficients des puissances de \, des fonctions de s qui 
résultent de l’intégration de certains déterminants. Le dénominateur 
est indépendant de s. La solution ainsi déterminée est unique. Il y a 
une difficulté pour les valeurs de \ qui sont poles de cette fonction 
méromorphe; la méthode montre que dans ce cas une solution n’existe 
pas en général, mais la méthode donne la solution dans les cas 
exceptionnels ot elle existe. 

Enfin une application convenable de la méthode — développée 
surtout par Hilbert et Schmidt — donne la solution cherchée sous 
une troisieme forme, celle d’une série de fonctions fondamentales. 
Ces fonctions sont dans les cas ordinaires les solutions de l’équation 
homogene 


b 
p(s) — af K(s, t) y(t) dt = 0. 
Cette équation n’est satisfaite en général que par 
p(s) = 0, 


mais il existe une suite de nombres — constantes caractéristiques (ou 
nombres fondamentauz (')) 


M1; dQ, ose An; 


pour chacun desquels cette équation a une solution finie, soit, 


y1(s), p(s), eth Yn(s), 


(1) D’aprés M. Goursar. 
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Ce sont les fonctions fondamentales. La solution de |’équation (1) 
s’obtient alors sous la forme d’une série 


p(s) = S an Yn(s). 


Cette méthode a été appliquée jusqu’ici surtout dans le cas ou 
K(s,t) est une fonction symétrique de s, t. 


K(s,t) = K(t, s); 


(dans ce cas les deux équations associées coincident) ('). 

Il faut remarquer qu’il n’y a rien de nouveau dans la notion de fonc- 
tion fondamentale, qui a eu son origine avec les séries trigonométriques 
de Fourier, et qui a occupé presque tous les grands géométres qui ont 
depuis étudié la Physique. Notamment M. Poincaré, a qui est dt le 
terme fonction fondamentale, a publié plusieurs mémoires profonds, 
dans lesquels il a traité des questions de la théorie du potentiel au 
moyen de ces fonctions. 

M. Hilbert (7) a appelé V’équation (1), une équation intégrale de 
seconde espéce, en réservant la désignation : équation intégrale de 
premiére espéce (3) pour l’équation 


b 
(3) [ Kod emat = se) 


L’équation 
b 


(4) ns) o(s) +2 [ K(s,) elt) at = f(s) 


a 


(') Voir la note B, p. 161 pour le cas général. 

(7) Nachrichten... zu Géttingen, 1904, n° 1. 

(3) La théorie de ces équations est beaucoup moins simple que celle des 
équations de deuxiéme espéce. En opérant comme au n°6, 2°, p. 43, on 


verrait facilement que méme dans le cas simple ot le noyau est de la forme 
q=p 

>> aq(s)Gq(t), Péquation de premiére espéce n’a pas, en général, de solution. 
q=1 
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qui se raméne immédiatement & (1) quand h(s) n’a pas de zéros dans 
Vintervalle d’intégration, est dite de troisiéme espéce(') dans le cas 
contraire. 

Nous ne nous occuperons pas en général de ces deux derniéres 
équations. 


2. — Le but principal de ce livre est d’exposer cette théorie 
au point de vue des applications physiques. Nous avons pour cette 
raison énoncé au Premier Chapitre un certain nombre de problemes de 
Physique qui conduisent a une équation de Fredholm. Ces problemes 
reviennent en général a chercher une fonction analytique a l’intérieur 
d’un domaine fermé, qui satisfait 4 une équation aux dérivées partielles, 
et qui se réduit a une suite de valeurs données a la frontiére de ce 
domaine. 

Dans le Troisieme Chapitre nous appliquerons 4 la résolution de 
ces problemes les résultats du Deuxieme Chapitre. 


3. Cas de plusieurs variables. — L’équation (1) définit une 
fonction y(s) qui dépend d’une seule variable s : ’extension au cas de 
plusieurs variables est immédiate, et il n’y a aucune modification a 
faire dans les démonstrations pourvu qu’on astreigne la frontiere de 
ce domaine d’intégration a des conditions de régularité convenables. 
En prenant par exemple le cas de deux variables, on a 


(81, 52) ff 8 (81, 82; t1, ta) y(tr, ta) d(t1, t2) = f(s1, 82). 


L’intégration s’étend a un domaine (D) dont un point est déterminé 
par les coordonnées (t;,¢2) : (D) est aussi le domaine ot sont définies 
les fonctions y(s1,s2), f(s1, 82). Il nous arrivera pour abréger de dire que 
(si, 82), f(s1,s2) sont deux fonctions du point M dont les coordonnées 
sont (s1,s2), et que K(s1,s9;t1,t2) est une fonction des deux points 


anes HILBERT, Nachrichten... zu Géttingen, 1906; E. PICARD, Comptes 
Rendus, 28 fév. 1910. 
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M(s1, 82), P(t1,t2) : on écrira léquation sous la forme 
(1" lM) A J K(M,P) p(P) dap = FM) 


ou dwp désigne un élément d’aire de (D). Dans cette formule pour 
simplifier l’écriture nous convenons que le signe { désigne cette fois 
une intégration double. 


4. Cas réel et cas complexe ('). — Nous ne nous occuperons pas 
des variables complexes; mais nous pouvons faire les remarques suivantes. 

La méthode d’itération et la méthode de Fredholm s’appliquent au 
cas général ott les fonctions K(s,t), f(s) et le chemin d’intégration sont 
complexes. 

La méthode et les résultats de Hilbert et de Schmidt ne s’appliquent 
pas toujours a ce cas général. 

Lorsque les fonctions K(s,t), f(s) sont réelles, et le chemin d’intégration 
réel, la solution et les fonctions fondamentales sont réelles. Les constantes 
caractéristiques sont réelles dans le cas symétrique, mais elles ne sont pas 
réelles en général. 


Dans la suite de Introduction nous énongons quelques définitions 
et propositions qui seront utiles plus tard. 


Rappel de quelques définitions 


5. — On dit qu’une fonction d’une ou de plusieurs variables 2, y,... 
est holomorphe pres du point (2o,yo,-.-), lorsqu’au voisinage de ce 
point on peut la représenter sous forme d’une série uniformément et 
absolument convergente de puissances entieres croissantes de x— x9, y— 
yo,---. ol la série reste absolument convergente quels que soient z,y,..., 
on dira que la fonction est entiére. 


(')On pourra sans inconvénient laisser de c6té lors d’une premiére lecture 
tous les paragraphes imprimés en petit texte. 
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Lorsqu’une fonction est holomorphe pres de tout point intérieur a 
un domaine D, on dit qu’elle est holomorphe dans D. Si la fonction, sans 


étre holomorphe dans D, peut étre considérée comme le quotient a 
2 


de deux fonctions holomorphes dans D, on dit qu'elle est méromorphe 
dans D. Une fonction méromorphe d’une variable est en général 
infinie pour certaines valeurs de \ appelées poles de la fonction. On 
démontre que les points représentatifs de ces poles sont en nombre fini 
dans toute aire bornée du plan complexe. Leurs valeurs sont racines 
du dénominateur fo et pres de lune quelconque d’entre elles, \’, on 
peut écrire la fonction sous la forme 


Ar AY 
pte ty + 9A) 


(5) £0) = Go yF 


ot: les A sont des constantes et ot y(A) est une fonction holomorphe 
pres de \’. La fraction rationnelle f —y est appelée la partie principale 


de f. 


5bis| Théoréme. — Soient f(s,t), g(s,t) deux fonctions bornées 
et intégrables dans le carré C : a < (s,t) < b, et qui sont continues 
dans C sauf peut-étre en certains points disposés de fagon qu’il n’y 
en ait jamais qu’un nombre fini ayant méme abscisse s ou méme 
ordonnée t. (Pour abréger, nous dirons d’une fonction qui posséde 
toutes ces propriétés qu’elle est continue presque partout dans C). Je 
dis que la fonction 


b 
(6) fey / feud 


est une fonction continue dans C. 

Il suffit de prouver que, quel que soit le point (s,t) fixe dans C, a 
tout nombre « > 0, on peut faire correspondre un nombre 7 tel que les 
inégalités 

Is—s|<n,  lé-tl<n 
entrainent dans C 
|F(s’,') — F(s,t)| <e. 
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Or s’il n’en était pas ainsi, on pourrait trouver un nombre « > 0 
tel que quel que soit n, il existe des valeurs sy, ty, dans (a,b) pour 
lesquelles 


1 1 
|s snl <_, |t tn] <—, |F(Sn,tn) — F(s,t)| > € 


Mais 
b 
F (sn, tn) — F(s,t) = | T(8n,T) g(r, tn) — g(t, t)| dr 
b 
+ J alrst)[(sn.7) fle.) ar 


Soient maintenant M, N les bornes supérieures de f et de g, on 
aura 


b 
(7) |F(5nstn) — F(s,t)| <M ; la(r, tn) — g(7,t)| ar 
+N f |flon7)- f(s,7)| dr. 


Considérons le coefficient de M; il n’y a qu’un nombre fini q de 
points de discontinuité d’ordonnée t; soient 71,72,...,7, leurs abscisses. 


En dehors des intervalles (x + 


é : , 
—..— }, la fonction g est continue sur 
16MNq )’ a 

la droite d’ordonnée t par rapport a l’ensemble de ses deux variables 
et par conséquent uniformément continue. On peut donc prendre r 


indépendant de 7 et assez grand pour que 
a 
4M (b — a) 
quand n>r, si 7 n’est intérieur 4 aucun des intervalles précédents. 
Il en résulte que la partie de la premiere intégrale de (7) relative 


|9(7,tn) — g(7,t)| < 


aux intervalles (x + } est inférieure a : 


E 
16MNq 


S 
ON x q = 
aunt EN SE 
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et la partie restante inférieure a 


E & 
b = ; 
Moog ae 
D’ot ’ 
e 
M tn) - t)|dr <M M= = =- 


pour n> r. On aurait de méme pour le second terme de (7) une 
inégalité analogue pour n>r’. Donc pour n>r +r’ 


|F (sn, tn) — F(s ES = 


i) 
ee 


contrairement a lhypothése. 


Fonctions orthogonales 
6. Inégalité de Schwarz. — Soient f(s) et y(s) deux fonctions 


continues dans l’intervalle (a,b) et \ un paramétre numérique, on a 
évidemment : 


b 
ies [ [ f(s) + Ay(s)]? ds 
- [ver Pasta fs) f(s ols)ds+% Tg (s)]? ds. 


En écrivant que le dernier membre est un trindme en \ jamais 
négatif, on obtient l’inégalité de Schwarz 


b 2 b b 
(8) / Hs) sh a <|[ (rte)? asx [ Iy(s)]? as. 


On obtient de la méme maniére 


(8’) [Ppn (s,t) ® dasa] <f fr (s, t)]| Pasar f° fe (s, t)]? ds dt. 
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7. Fonctions orthogonales. — On dit que deux fonctions continues 
f(s), o(s) sont orthogonales dans un intervalle (a,b) si ’on a 


6 
i f(s) p(s) ds =0. 


Elles sont normales si 


[isoPa =i [teres 


a 


Il est facile de voir que si des fonctions continues 


yi(s), p2(s), tees Yq(s), eid 


sont normales et orthogonales(!), ces fonctions sont linéairement 
indépendantes. Autrement dit, une identité de la forme 


€1 y1(s) +... + eg Yq(s) +... = 0, 


ot. les c sont des constantes et ot le premier membre est, soit une 
somme d’un nombre limité de termes, soit une série uniformément 
convergente, n’est possible que si ces constantes sont toutes nulles. En 
effet, en multipliant par y;(s) et intégrant, on a 


b b 
o=a | eils)ei(s)ds+.. +6 f [yi(s)]? ds +... 


b 
+04 f Pq(s) vils)ds+...= cj, 
a 


POUR FS 1a Won 


C) Comme exemple d’un tel systeme, on peut prendre : 
COS $ sin s cos ns 


1 
gi(s) = Vin’ p(s) = 7a y3(s) = aa -++)Pan(s) = Se Pan+i($) = 


sinns 


., avec a=0, b= 2r. 
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8. Normalisation. — Etant donné un systéme de fonctions conti- 
nues dans (a,b) et non identiquement nulles 


fils), fo(s), -.., Fn(s), «-- 


en nombre fini ou non, on peut toujours normaliser ce systeme, 
c’est-a-dire trouver un systeme de fonctions continues dans (a,b) 


y1(s), y2(s), ae 


qui sont normales et orthogonales et telles que toute combinaison 
linéaire des f(s) est une combinaison linéaire des y(s) et inversement. 

En effet, supprimons de la suite des f, celles qui sont une 
combinaison linéaire des précédentes et nous serons ramenés au cas ou 
les fonctions f sont linéairement indépendantes quand on en considére 
un nombre fini quelconque. 

Supposons alors démontré qu’on peut trouver p combinaisons 
linéaires de fi, fo,...,fp, Solent 91, ¢%2,--.,~p, qui sont normales et 
orthogonales entre elles. Ce sera aussi vrai quand on remplace p 
par p+1. En effet, formons 


(9) Pptl =C1i vit... + Cp Pp + Cp4i fpti 
ce sera une combinaison linéaire des f. Elle sera orthogonale a ¢1,..., %p 
si ’on a 


b b 
o= | ep+i(8) PR(8) d8= cy + ope | fp+i(s) pe(s)ds, (kK =1,...,p). 


D’ot 


P b 
Yor = epi[fp4i - Yvxts) f fp+1(8) px(s) ds]. 
1 Ja 
Le crochet n’est pas identiquement nul, sans quoi fp, serait 
une combinaison linéaire de ¢1,...,p~) et par suite de fi, fo,..., fp 
ce qui serait contraire a l’hypothése. Alors, on pourra évidemment 
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choisir cp; de fagon que yp}; soit normale, ce qui démontre notre 
remarque. Or celle-ci admise pour p fonctions est vraie pour une; elle 
est donc générale. Le théoréme est alors démontré, car les fonctions 
orthogonales et normales ¢,y2,... sont nécessairement indépendantes, 
et, non seulement elles sont bien des combinaisons linéaires des f, 
mais d’apres (9), les f sont aussi des combinaisons linéaires des y. 
De plus nous voyons que si les f sont linéairement indépendantes, le 
nombre des y sera égal a celui des f, en particulier s’il y a une infinité 
de fonctions f, il vy aura une infinité de fonctions y. 


9. Constantes de Fourier généralisées. — Considérons un systeme 
normé, c’est-a-dire un systéme de fonctions normales et. orthogonales 
entre elles 

1(s), 92(s), 93(s), 

Il est facile de voir qu’une fonction F(s) qui peut étre considérée 

comme une combinaison linéaire d’un nombre fini ou non de fonctions 


de ce systéme ne peut etre exprimée ainsi de plus d’une facon. D’une 
maniére précise si la série 


dy p1(s) + dz yo(s) +... +dn vn(s) +... 


est composée d’un nombre fini de termes ou converge uniformément, 
les coefficients supposés constants d,,d2,... sont bien déterminés par 
la somme F(s) de cette série. En effet, en multipliant par y,(s) et 
intégrant, on aura 


= [ "Hihentelds 


Par analogie avec ce qui se passe dans les développements trigo- 
nométriques, nous nommerons les d,, les constantes de Fourier de F(s) 
relativement au systeme des yp. 


10. Inégalité de Bessel. — Soit F(s) une fonction continue dans (a, b) 
et 


Yl, 2; 
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un systéme de fonctions continues normales et orthogonales dans (a,b). 
Posons 


b 
(10) ap= [ F(s) ¥(s) as 


On aura 


b b b 
i. [F(s) - yy Ap er()| di = i) [F(s)]? ds — 2 .s ap | F(s) p(s) ds 
b b 
+ s— of | ys (s) ds +2 ye Apr i p(s) Yr(s) ds 
b 
= | [F(s)]? ds — s ae. 
D’ot une formule importante 


Sa = [ [Fe yPas— f [F(s) — S= ap ¥p(s)] ds. 


En particulier on a l’inégalité de Bessel 


b 
(11) va < f [F(s)]? ds, 


ou les a, sont les constantes de Fourier (10) de F(s) relativement 


aux y(s). 
p=00 
Il en résulte que si les yp sont en nombre infini, la série S~ a7 sera 
b p= 1 


convergente et inférieure ou égale a | [F(s)]? ds. 
a 


11. — On peut en tirer avec M. Schmidt la conséquence suivante. 

Les constantes de Fourier d’une fonction continue de deux va- 
riables H(s,t) ot. l’on considére t comme constant, sont des fonctions 
de t : 3,(t). Je dis que si H(s,t) est une fonction symétrique, la série 


Soap By(t) => |r (yas H(s, t) pp(s) ds 
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est absolument et uniformément convergente quand ¢ varie dans (a,b). 
En effet, on a 


p=n+m b [p=n+m ‘p=n+m 
J= Xe [op Bp(t)] =) | Se Ap an) S- p(t) ov) ds 
pHn p=n 


p=n a 


puisque les » forment un systeme normé. D’apres les inégalités de 
Schwarz et de Bessel, appliquées au dernier membre, on a 


=n+m 
Jee f ‘> 3 x [ttn] a 


p=n 


Mais la fonction H(s,t) étant continue, la seconde intégrale a un 
maximum M lorsque ¢ varie dans (a,b). On a donc : 


p=n+m p=oo 
dle Sp(t)| < ISI < | M >> a3. 
p=n p=n 


Comme la série >a? est convergente, le second membre (et par 
conséquent le premier membre) peut étre rendu aussi petit que l’on 
veut, en prenant n assez grand et m quelconque, de sorte que la 
proposition est démontrée (voir la note ('), p. 102). 


CHAPITRE PREMIER 


PROBLEMES SE RAMENANT A L’EQUATION 
DE FREDHOLM 


I. — PROBLEMES DE POTENTIEL 


1. Fonctions harmoniques ('). — On dit qu’une équation linéaire 
aux dérivées partielles du second ordre (*) 


Vv av 
ye a ae dM a IV = f(#1,%2,...,%n) 


appartient au type elliptique, ou encore, est a caractéristiques (7) 
imaginaires, si la forme quadratique 


Sap XiXx 
ik 


(obtenue en remplacant, dans les termes du second ordre, chaque 


dérivée par un produit de deux indéterminées X;, X;) est 


. 
définie, c’est-a-dire se compose d’autant de carrés indépendants et 
tous de méme signe qu’il y a de variables. 

La plus simple et la plus connue des équations du type elliptique 
est léquation de Laplace ou équation des potentiels. 


1) Pv Ov &V 
Ox2 " Oy? ' Az? 


=0 ou AV=0 


!) HADAMARD, Propagation des Ondes, ch. I. 
g 
(7) HADAMARD, Ibid. Chap. VII, § 3. 


RECN! 
RECN! 


ROCCNS 
ROOCNe 
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On appelle fonction harmonique dans un domaine D & trois 
dimensions toute solution de cette équation analytique(') dans le 
domaine D. 

Dans le cas ot! D est illimité, il faut de plus que le rapport de V 
a =: R étant la distance de Vorigine, reste inférieur en module 4 un 


nombre donné (ainsi que le rapport de chacune des dérivées partielles 
OV OV OV ‘ 1 ) 
Ox’ Oy’? Oz R2/° 

2. Potentiel. — Les fonctions suivantes sont harmoniques en 
dehors des masses attirantes : 

1° Le potentiel d’une distribution spatiale 


@) w= fff © aeayee 


2° Le potentiel d’une simple couche étalée sur une surface S, 


(3) U(M) = i; oP as 


3° Le potentiel d’une double couche, 


woo ffan[i Qj s- [forsee 


ot. M désigne le point dont on considere le potentiel, P, un point 
attaché a un élément de masse, r = MP et ot ¢ est angle de PM et de 
la normale a la surface $ au point P, dirigée vers l’intérieur (fig. 1). 

On suppose que S a partout un plan tangent unique; la masse 
attirante est supposée étre située a une distance finie de lorigine. 
La densité p est une fonction finie et intégrable : elle peut étre 
discontinue. 


Cyn sufft de supposer l’existence de V et de ses dérivées premiéres et 
secondes. 
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A Vintérieur des masses attirantes, on a la formule de Poisson 


(opis) AV = —4r p. 


3. Continuité du potentiel.—Dans ces conditions : 

1° Le potentiel d’une distribution spatiale est partout continu, 
ainsi que sa dérivée premiére. Sa dérivée seconde a en général des 
discontinuités aux points ott p est discontinu. 

2° Le potentiel d’une simple couche est continu et nul a Vinfini. 
Sa dérivée premiére (continue ailleurs) a des discontinuités sur la 
surface S, sur laquelle les équations suivantes sont satisfaites, 


1f/oV ov 
1/o0V_ | Ov cos w 


Vo- Vm 
is Si pee ee 

lorsqu’en partant d’un point Q sur la normale, extérieur a la surface, 
on s’approche du point M de la surface le long de la normale. L’indice i 


Ici Vindice e signifie que ov est la limite du rapport 
nr 
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signifie qu’on prend la méme limite lorsque le point Q est intérieur 
a la surface. Le signe f désigne une intégrale double étendue a la 
surface S$, dont do est I’élément d’aire. ~ est l’angle de MP et de la 
normale a M. 

3° Le potentiel d’une double couche (continu ailleurs) a des 
discontinuités sur la surface S. Il est nul a Vinfini et satisfait aux 
équations suivantes : 


1 
(7) sia — We) = 2m p(M), 
1 Q 
(8) a(Wit We) = f oe) SSP as. 
2 S r 
4. Problemes relatifs aux fonctions harmoniques. — Le 


probleme de déterminer une fonction V harmonique dans un domaine D 
extérieur ou intérieur a une surface S et satisfaisant a diverses 
conditions sur la surface S de ce domaine a une importance capitale. 
Il est nommé suivant les différents cas : 
1° Probleme de Dirichlet. — On impose a la fonction V cherchée 
la condition d’étre égale sur S 4 une fonction donnée sur $ 


(9) Vs = f(M). 
2° Probleme de Neumann. — On donne sur S$ non plus les valeurs 
; OV 
de V mais celle de — 
On 
Ov 
(10) ang g(M). 
é ‘ ; OV : 
3° Probleme de la chaleur. — L’expression are est donnée 


sur 8, p étant une fonction donnée de M 


(11) im + p(M) Vg = A(M). 
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La condition 
OV 


1 On 
se raméne a la précédente tant que q 40. 

Le mot «probleme mixte» s’applique a un probleme qu’on rencontre 
en particulier en Hydrodynamique et correspond a q = 0 sur une partie 
de la surface, et p=0 sur autre ('). 

A chacune de ces conditions correspondent deux probléemes distincts 
— intérieur et extérieur. 


+ pV =h(M) 


5. Problemes généralisés. — Sont presque aussi importants que 
les précédents, les problemes qui consistent a trouver non une fonction 
harmonique, c’est-a-dire une solution analytique de 


AV =0 
mais une solution analytique de 
(12) AV = y(M) 


satisfaisant a l’une des conditions 1°, 2°, 3° sur une surface S. 
Ces problemes se ramenent aux précédents. Cherchons, en effet, 
une solution V de 
AV = (2, y, 2) 


dans un domaine D, qui, sur la frontiere S de ce domaine, satisfait a 
la condition 


V = f(M). 
Pour cela il suffit de trouver deux fonctions V, et V2 telles qu’on 
ait 
AVi = Ys 
ree 1 = (2, y, 2) 
AV2 =0 
et sur S: Vo = f(M)— Vi(M). 


(1) HADAMARD, Propagation des ondes. 
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La solution sera 
V=Vi4+ Vo. 


Or, grace A la formule de Poisson (2%), on a tout de suite une 
solution en V, 


Vi= ae 3 y(P) - dup 
ou r = MP. (Dans le cas du probléme plan, on remplacerait y par logr). 
Et la recherche de V2 revient au probleme de Dirichlet. 
Les deux autres problemes correspondant au probleme de Neu- 
mann, au probleme de la chaleur et au probleme mixte se traitent 
dune maniere analogue. 


6. Problemes de physique mathématique correspondant 
aux précédents. 1° Attraction newtonienne et électrostatique. — 
Chacun des cas du n° 4 correspond a un probleme de ces théories. 
On peut mentionner (!) celui de trouver le potentiel 4 l’extérieur ou 
a Vintérieur d’un conducteur sur lequel s’étale une charge donnée. Le 
probleme revient a chercher une fonction harmonique dans D et égale 
sur sa frontiére a une constante donnée. 

Le cas ot il y a plusieurs conducteurs auxquels on commu- 
nique différentes charges se réduit au cas précédent. On peut en 
effet considérer les différents conducteurs comme une seule surface 
multiconnexe. 

On peut mentionner aussi le probléme des marées (7). On considére 
la terre comme une sphere dont certaines parties, correspondant aux 
mers, sont conductrices. Cette sphere est placée dans un champ qui 
est produit par les attractions des astres, et par la rotation de la 


(1)O. D. Cuwo son, Traité de Physique (traduit par DAVAUX), 1910, Paris. 
Tome IV, Champ électrique constant, p. 129. 

(7) Darwin, Philos. Trans., 1879, pp. 447, 539; 1880 (IL), p. 713; 1881 (ID), 
p. 491; 1887, p. 379. 
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terre elle-méme (force centrifuge). La densité électrique sur la sphére 
correspond a la hauteur des marées ('). 

2° Magnétisme. — La fonction V peut également étre un potentiel 
magnétique, et l’on peut traiter le probleme de l’aimantation par 
influence (7). La Physique nous donne une relation entre la force 
magnétique extérieure et intérieure sur la surface du corps aimanté 


OV OV 
On; One 


(1 + 47k) = fonction donnée, 
k étant le coefficient d’aimantation. 

Ce probleme differe un peu des problemes des n°%*4 et 5, nous 
verrons dans le n°9 qu’on le résout de la meme maniere. 

3° L’hydrodynamique. — La fonction V est maintenant le potentiel 
permanent des vitesses (?). Le probleme du mouvement irrotationnel 
d’un liquide dans un vase fermé immobile que ce liquide est supposé 
remplir exactement correspond a la condition 


OV 
On 
Si les parois du vase se déplacent (le vase restant toujours 
complétement plein), on a 


oN = fonction donnée. 
On 


Le mouvement des parois est assujetti a la condition que le volume 
du liquide reste constant, de sorte qu’on a 


[Garo 
g dn 


(') Poincar, Journal de Mathématiques, 1897. 
(7) Mascart et JouBERT, Electricité et Magnétisme, t. I, § 270. 
(3) APPELL, Traité de Mécanique, t. III, p. 325, 327. 
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On obtiendrait un probleme différent, cas particulier du probleme 
de Dirichlet, en supposant que la direction de la vitesse en chaque 
point d’une surface fermée est normale a cette surface. Celle-ci devant 
étre alors une surface de niveau, cela revient a se donner 


V =constante sur la surface. 


Si enfin on a affaire a un liquide présentant une surface libre, on 
doit déterminer V en se donnant ses valeurs sur une partie d'une 


surface et celles de ov sur autre partie. Ce cas se ramene donc au 
probleme mixte. . 

4° EKlastostatique. — Considérons le cas d’une petite déformation 
d’un corps de forme donnée, et supposons que cette déformation dérive 
d’un potentiel V('). On a 


AV =9@ 
ot @ est la dilatation. Si des forces données agissent sur la surface on 
e Ov , , 
connaitra —. La fonction @(z,y,z) sera donnée par la force de masse. 
nr 
On a donc un cas du probleme généralisé. Par contre, le probleme 
est notablement plus difficile lorsqu’on ne fait pas sur la solution 


Vhypothese précédente. 
5° Théorie de la chaleur. — Le probleme d’un corps en équilibre 


(*) Soit (x,y,z), (© +€,y+n,2+6¢) les coordonnées d’un point P du corps 
avant et aprés la déformation due a l’action des forces de masse et des forces 
de surface. Notre hypothése s’exprime par les relations : 


OV OV OV 


$= 3, a Ou ae 


La dilatation @ est le rapport de l’accroissement de volume au volume primitif : 
ona: 


V. APPELL, Traité de Mécanique, t. III, p. 805. 
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de température avec rayonnement (') nous donne le probléme du n° 4 
avec la condition 


ov 


= fonction donnée. 
On; 


kV — 


C’est la loi de Newton. 


7. — Relativement a chacun des problemes des n°*5 et 6 se posent 
deux questions préliminaires. Est-ce qu’une solution existe? Si elle 
existe, est-elle unique? La méthode donne une réponse tres nette a 
chacune de ces questions. 

Il faut remarquer que la Physique résout ces questions dans 
les cas cités au n°6. Il est physiquement évident que les divers 
équilibres proposés seront réalisés, et qu’ils ne seront réalisés que 
d’une seule maniére. Mais cette réponse ne suffit pas au point de vue 
mathématique et nous aurons a la compléter. Pour cela, il nous sera 


souvent commode d’utiliser la notion de fonction de Green. 


8. Fonction de Green. — Soient M, P deux points du domaine D : 

be ae 
envisageons la fonction des deux points, — = —~ (*). 

Considérée comme fonction de P, elle est analytique en tout point 
sauf en M, et si M n’est pas sur la surface S, elle prend des valeurs 
déterminées en tout point de cette surface. 

Imaginons qu’on puisse construire une fonction @ harmonique dans 
le domaine D (limité par S et intérieur ou extérieur 4 S) et prenant 


1 1 4 ; 
sur la surface S les valeurs i ge M étant pour le moment un point 
r 


fixe. Ceci revient a la solution d’un probleme de Dirichlet de l’espéce 
indiquée au n° 4, p. 17. 


(1) BoussinEsq, Théorie analytique de la Chaleur, t. I. CHWOLSON, loc. 
cit., t. III, p. 303 et 331. 


1 1 
(7) Dans tout ce numéro, on devra remplacer — par log — pour obtenir le cas 


du plan; il faudra en outre remplacer 4x par 27 dans la formule (13). 
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La fonction harmonique de P 


1 
qui s’annule sur S et qui a une singularité au point M s’appelle 
fonction de Green relative a D. 

Elle est symétrique en M et P, c’est-a-dire que 


6(M, P) = G(P,M). 


Dans le cas ot M et P se trouvent a la fois sur la surface, 
la fonction G s’annule également, pourvu que les deux points ne 
coincident pas. 

On démontre d’ailleurs qu’une fois connue la fonction de Green la 
résolution du probleme de Dirichlet s’ensuivrait. Mais nous n’utilise- 
rons pas ce fait. 

Au contraire, nous nous servirons d’une autre propriété remar- 
quable exprimée par la formule 


(13) Af #@)G(M.P) dep = ~4nf() 


dans laquelle f(M) est une fonction finie quelconque ayant des dérivées- 
premieres et ot le signe { désigne une intégration triple effectuée dans 
le domaine D dont |’élément de volume est dwp. 

Pour démontrer cette relation, on écrit le premier membre sous la 
forme 


1 
er) dup - A f (M,P). s(P) dor. 


La seconde expression s’évanouit, puisque w est harmonique. 

La premiere intégrale est le potentiel d’une distribution de matiere 
de densité f(P). 

Donc le premier membre de (13) est bien égal a 


—4m f(M) 
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& cause de l’équation de Poisson (2'S), p. 15, a laquelle satisfait la 
premiere intégrale, pourvu que f(M) possede des dérivées premiéres 


Of of of 


Ox’ Oy? Az’ 


Enfin nous établirons aussi la propriété suivante : 
L’intégrale 


(14) J £@)GO4.P) dep 


s’annule lorsque M s’approche de la surface. 

En effet, tous les éléments de cette intégrale s’annulent sauf celui 
qui contient le point M. 

Prenons un point A sur la surface (fig. 2) et autour de A décrivons 
deux spheres cy, cp de rayons r1, r2. Considérons la partie de l’intégrale 
qui s’étend 4 l’espace © entre ces deux sphéres et la surface S (*). 

Soit M un point situé a une distance de A plus petite que 

La fonction harmonique @ n’a ni maximum ni minimum, et elle est 


, « , : 
égale a — sur S; par suite elle reste comprise 
r 


entre le maximum et le minimum de = sur 

la surface. Elle est donc positive. ‘ 
D’autre part G = : —q@ est une fonction 

de P harmonique, sauf au point M, qui 

s’annule sur S, et qui est strement positive 

au voisinage de M. En tracgant autour de M 

une sphére trés petite, on voit qu’entre cette 

sphere et S, gG a un maximum positif sur Fie, 2 

cette sphere et un minimum nul sur S. Par 

conséquent, G est positive partout, et, puisque w > 0, elle est en 


(1) Nous avons substitué dans la figure 2, «contour et aire renfermés» a 
«surface et volume renfermés». 
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module plus petite que a Si K est le maximum de f(P) dans D, on a 


< K [- dwp, 
e étant une partie quelconque de D. 

Il faut démontrer que l’intégrale (14) tend vers zéro quand le 
point M s’approche d’un point A du bord. Choisissons d’abord un 
petit nombre positif «. Nous allons démontrer qu’on peut prendre le 
point M assez prés de A (fig. 2) pour que l’intégrale (14) soit plus 
petite que «. 

Autour de A décrivons une petite sphere c, de rayon p. Soit M a 
Vintérieur de cette sphere. Nous pouvons choisir le rayon p tel que 


donc : 


(15) G(M, P) dwp 


ae. 1 Z & Wiecpet . 
Vintégrale / —dwp étendue a l’intérieur d’une sphere c2 de centre M et 
iF 


de rayon 2p soit plus petite que =) r étant la distance entre le point P 
et le centre. En effet cette derniére intégrale est égale a 


2p 1 
[ = Arr? dr = 87”. 
0 


Remarquons que p étant ainsi choisi et fixé, la sphere fixe c; reste 
tout entiere comprise dans la sphere variable cj, si M est a l’intérieur 
de la sphere fixe c, décrite autour de A. L’intégrale (15) relative a 
la partie « du domaine D a l’intérieur de cette sphére fixe c, sera 
donc plus petite que =. Les éléments de Vintégrale (15) extérieurs a la 
sphére c, sont finis et ils tendent vers zéro lorsque M s’approche du 
bord, c, restant fixe. Donc on peut rendre le reste de l’intégrale (11) 
plus petit que 5 en faisant tendre M vers A. 


Donc on peut rendre l’intégrale (11) plus petite que = + = =e, 
Nous venons de traiter un cas de la fonction de Green, dont 


Vexistence dépend de la résolution du probleme intérieur de Dirichlet. 
Il y a aussi a considérer la fonction de Green relative a l’extérieur 
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d’une surface, et ensuite des fonctions analogues pour les conditions 


sur la surface et 


sur la surface pour les cas extérieur et intérieur. Pour approfondir 
cette question nous avons besoin des résultats du Deuxiéme Chapitre. 
Nous reconnaitrons (n° 14, p. 132) a Vaide de ces résultats que des 
fonctions de Green satisfaisant a ces conditions n’existent pas dans 
tous les cas. Dans les numéros suivants (n° 13, p. 30) nous indiquerons 
Vemploi de certaines de ces fonctions de Green. 


9. La mise en équation. — Nous allons maintenant montrer 
comment a une exception prés (§ 4°, p. 21) les problémes précédents 
se rameénent aisément a des équations de Fredholm, c’est-a-dire de 
la forme (1'), p. 4. Envisageons d’abord le probleme 1° du n° 4, 
c’est-a-dire le probleme de Dirichlet. Nous le remplacerons par un 
probleme plus général : 

Trouver une fonction V harmonique dans tout l’espace (sauf sur la 
surface S$) et qui sur la surface S satisfait & la condition 


(16) VituVe =huf(M), 


ou pw est un parametre arbitraire. 

Pour » = 0, nous avons le probleme intérieur; pour p = oo, le 
probleme extérieur. 

Remarquons maintenant que, étant donnée une distribution de 
matiere, une intégration suffit pour déterminer son potentiel. Donec si 
nous pouvons déterminer une distribution dont le potenticl satisfait a 
notre condition, le probleme sera résolu. Nous chercherons une double 
couche (!)de densité p étalée sur 8 et dont le potentiel W (= V) satisfait 
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a la condition (16). 

Les équations (7) et (8) sont maintenant satisfaites. Et si nous 
substituons dans l’équation (16) les expressions de V; (= W;) et 
de Ve (=W-) tirées de (7) et (8) nous aurons une équation renfermant 
la densité p cherchée. 


(17) 2m p(M) +2 f p(P) P52 as = f(a) 
ou 1+ 
d= a fi(M) = To 1. 


C’est une équation de Fredholm (n°1, p. 1) dont le noyau 
(note (+), p. 1) est est fonction des deux points P, M. Si nous 
pouvons résoudre cette equation, nous aurons une densité p qui nous 
conduira immédiatement a la solution du probleme posé. 

Les valeurs \ = 1 et \ = —1 du paramétre correspondent aux cas 
intérieur et extérieur respectivement. 

Le probleme 2°, celui de Neumann, se traite d’une maniére 
analogue. 

Nous le remplacgons par le probleme plus général ot on a la 
condition 


OV i OV 
P On, One an 


(18) (M). 

Cherchons une simple couche p étalée sur S et dont le potentiel 
satisfasse a cette condition. En nous servant des équations (5) et (6) 
nous trouvons que p doit satisfaire a l’équation de Fredholm suivante 


cos w 


2 


(19) 2m p(M) +2 f p(P) 5" ds = fi(M), 


(*) On pourrait se demander pourquoi nous choisissons une double couche. 
C’est que cette espece de distribution nous permet d’appliquer la méthode de 
Fredholm. En représentant V par un potentiel de simple couche, on serait 
conduit 4 une équation de premiére espéce (voir n° 1, p. 3). 


L’EQUATION DE FREDHOLM 28 


ou 
bb 
fi (M) = —— g(M) 
l~pu 
et A a la méme expression qu’avant. 
Le noyau est maintenant ~ see : et lon voit que l’équation (19) 
Tr 
est associée (p. 1) a l’équation (17). 
Les valeurs 4 = 1 et } = —1 du paramétre correspondent aux cas 


extérieur et intérieur respectivement. 
Nous arrivons maintenant a 3°, le probleme de la chaleur. Cher- 
chons une simple couche dont le potentiel remplit la condition (11). 


En employant les équations (3), (5) et (6) nous avons l’équation 
intégrale 


(19bis) Qn p(M) — df oP) ss ) dS = h1(M), 


oll hy = —Ah et X = +1,—1, selon que nous sommes dans le cas intérieur 
ou le cas extérieur. 
[L’équation (198) correspond a la condition 
1 OV 1 OV 


s(t +9) Fa 34-1) 5 - +A e(M) V(M) = An(M), 


Cette méthode ne s’applique d’ailleurs pas au cas ou la quantité 
appelée plus haut g (n° 4, p. 17) s’annule et par conséquent pas non 
plus au probleme mixte. 


10. Cas de deux dimensions. — Tout ce qui précéde s’applique 


: : ‘ : ee 
au cas de deux dimensions mot pour mot, si l’on substitue log— a -. 


Les problemes qu’on se propose sont identiques, et ils conduisent a 
des équations fonctionnelles analogues. On forme de la méme fagon 
les fonctions de Green pour un contour plan. 


REESE RROOOe 
REESE REO 


NN 
REER—— BOO 
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II. — PROBLEMES RELATIFS A L’EQUATION AV =R(«,y,2)V 
ET A DES EQUATIONS ANALOGUES 


11. Les probleémes envisagés. — Nous nous occuperons dans 
cet article des problemes qui consistent a trouver une solution de 
Véquation (également elliptique) 


(20) AV = R(a, y, z)V 
ou bien de ’équation plus générale 
(21) AV = R(z, y, 2)V + Ri (a, y, 2) 


ou. R(z,y,z), Ri(z,y,z) sont des fonctions finies, ayant des dérivées 
premiéres (R gardant toujours le méme signe), cette solution étant 
assujettie a l’une des trois conditions de Neumann, de Dirichlet, de la 
chaleur, indiquées dans le n° 4. 

La résolution de ces problemes nous permettra (n° 19, Ch. II, 
p. 140) de traiter une autre catégorie de questions, celles qui consistent 
a trouver une solution de 


(22) AV =R(e,y,2) 
ou bien de 

a2v 
(23) AV = Riz, y, z) aD 


satisfaisant a lune des conditions précédentes sur la surface, et 
satisfaisant a des conditions initiales. 
Enfin nous montrerons (n° 24, p. 152) que la méme méthode s’ap- 
plique a l’équation générale elliptique a deux variables indépendantes 
ON OV «OV OV 


Dre + Dy Br aca +cVa=f 


ou a, b, c, f sont des fonctions données de x, y (c’est-a-dire de M). 
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12. Problémes de Physique correspondants. 1° Acoustique ('). 
— Probleme des ondes sonores. La fonction V est le potentiel des 
vitesses : elle satisfait a l’équation (23). La constante R dépend de la 
densité et du coefficient d’élasticité du gaz. Donner la vitesse normale 


autour de la surface revient 4 donner x. Supposer que la vitesse est 
normale revient a se donner V. 

2° Elasticité (t). — Pour le probleme des petites vibrations d’un 
corps élastique, on a les mémes équations que dans le 1°. La fonction V 
est la dilatation. 

3° Chaleur (7). — Pour le probléme du_ refroidissement d’un 
corps, nous avons |’équation (22) ot V représente la température. La 
fonction R est le rapport de la chaleur spécifique au coefficient de 
conductibilité. Si la température a la surface est connue, nous avons 


V = fonction donnée sur la surface. 


Si on donne le flux, nous avons 
OV : : 
aa fonction donnée sur la surface. 
nr 


Sil y a de la radiation, nous avons 


R ~ —V = fonction donnée sur la surface 
nN 


d’apres la loi de Newton, ot R est fonction des points de la surface. 


13. Mise en équation. — Considérons d’abord |’équation (20) et 
soit G(M,P) la fonction de Green correspondant a la condition 


(24) vV=0 sur la surface. 
() Lord RAYLEIGH. — Theory of Sound. Vol. II, n°°241-244 et n° 374; 
O. D. CHWOLSON. — Traité de Physique (traduit par E. Davaux). Acoustique, 


p. 103; A. WINKELMANN. — Handbuch der Physik. Bd. I, p. 206, Bd. I, p. 146; 
M. BOuSSINESQ. — Théorie analytique de la Chaleur, t. I, p. 47. 

(7) Boussinesq. — Loc. cit., t. I, p. 3 et p. 64 (suite); KIRCHOFF. — 
Theorie der Warme, § 6; A. WINKELMANN.—Loc. cit., Bd. III, p. 440. 


i 
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Si Pon a 
1 
(25) V(M) = erie [ G(M, P) R(P)V(P) dwp 


il est évident d’aprés la formule (13), p. 23, que V_ satisfera a 
Véquation (20) et réciproquement. En tenant compte des résultats du 
n° 8, on voit que V satisfera également a la condition (24). Donc la 
question revient a résoudre |’équation fonctionnelle homogéne (25). 

Nous avons une méthode analogue pour les solutions correspondant 
aux conditions 


OV 
(26) A 0 
et 
OV 
(27) Sy t PM) V =0. 
Considérons ensuite la condition 
(9) V = fonction donnée = f(M) sur la surface. 


Soit v(M) une fonction harmonique qui satisfait 4 la condition (9). 
On voit que V est la solution de l’équation fonctionnelle non-homogene 


(28) V(M) = 7% [ G(M,P) R(P) V(P) dep + v(M). 


Et de méme pour les conditions 


(13) oY = o(M) 
(14) oY + p(MV = A(M). 


Les solutions de (21) correspondant aux diverses conditions satisfe- 
ront a l’équation 


(29) V(M) = = Ke €(M, P) [R(P) V(P) + Ri (P)] dp + v(M). 
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14.—Arrivons maintenant aux équations (22) et (23) qui renferment 
la nouvelle variable t. Ce sont — contrairement aux premieres — 
des équations du type hyperbolique ou parabolique et les problemes 
correspondants ne sont pas entierement analogues a ceux du numéro 
précédent : ils en different par Vintervention des données initiales 
(pour t = 0); mais ils s’en rapprochent par les données aux limites 
(c’est-a-dire sur la frontiere du domaine en 2, y, z). 

Nous les résoudrons par un artifice qui fait intervenir au lieu 
des équations hyperboliques ou paraboliques (22), (23), l’équation (20) 
laquelle est du type elliptique. Les solutions s’exprimeront sous forme 
de séries que nous démontrerons plus tard (n° 19, Ch. III, p. 144) étre 
uniformément et absolument convergentes. 

Substituons dans (22) V, =e y(a,y,z); il vient 


(20bis) Ay = -\R(z, y, z)9. 


Donc V, sera une solution de (22), si on prend pour y une solution 
de (20>) qui est de la forme (20). 
Prenons d’abord le cas ot la fonction V doit satisfaire a la condition 


V=0 


sur la surface et cherchons une fonction y, solution de (20>) 
satisfaisant a 

yp=0 
sur la surface. 

Nous venons de voir que ceci revient a résoudre une équation ho- 
mogeéne (25) — ot l’on remplacerait R par —\R. — Nous démontrerons 
plus tard qu’il existe une solution seulement pour une suite infinie de 
valeurs de \ 

M1, AQ, +--+) Any --- 


qui sont toutes réelles et positives. Nous appelons ces quantités, 
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constantes caractéristiques (1) et nous écrivons les solutions correspon- 
dantes, 


Pl, P22, «+5 Pny -es 


que nous appelons fonctions fondamentales (!). Donc nous aurons la 
solution suivante de (22), qui satisfait a (24) 


(30) V= Yo Ane™ nt on (x,y, 2) 


ou les A, sont des constantes arbitraires, qui doivent étre telles 
que la série (30) soit absolument et uniformément convergente, si les 
données sont continues. Enfin il faut que V se réduise 4 une fonction 
donnée a(z,y,z) pour t = 0. C’est-a-dire, il faut que Aj, Ag,...,An,--. 
soient choisies de facon que 


(31) (x YU, -YoMs nl L,Y,% z). 


Nous démontrerons que ce développement est toujours possible, et 
qu’on peut déterminer les coefficients d’une maniére unique analogue 
a la maniére de calculer les coefficients de Fourier. 

Remarquons qu’il est nécessaire que la fonction a(z,y,z) satisfasse 
a la condition 

a(z,y,z) =0 


sur la surface. 

Au point de vue physique, il est important de considérer le 
cas ot. cette condition n’est pas remplie. C’est par exemple le cas 
du refroidissement d’un corps dont la distribution thermique est 
quelconque et qu’on plonge dans de l’eau glacée. On prend dans ce 
cas pour a une fonction qui est égale 4 la température initiale dans 
le domaine ouvert D, mais qui prend la valeur zéro sur la surface. 
La théorie n’est pas encore approfondie dans les cas discontinus, mais 


(1) D’aprés M. Porncare. 


L’EQUATION DE FREDHOLM 34 


nous pouvons prévoir par analogie avec ce qui se passe pour la série 
de Fourier, que la série (31) sera encore valable et qu’elle sera non- 
uniformément convergente sur la surface. Et dans ce cas nous savons 
par un théoréme de Cauchy que V sera analytique (donc continue) 
pour ¢ > 0, de sorte que la série (31) sera toujours uniformément 
convergente. 

En somme tous ces résultats sont admis par les physiciens depuis 
longtemps. 

Nous avons des résultats tout a fait analogues pour les conditions 
(26) et (27). 

Nous considérons maintenant la condition 


(9) V = f(M) 


sur la surface. 

Soit v une fonction harmonique qui satisfait a (9). La fonction 
a(M)—v(M) s’annule sur la surface. Trouvons par la méthode précédente 
une fonction V;, satisfaisant a (22), s’annulant sur la surface, et se 
réduisant & a—v pour t=0. 

On voit tout de suite que 


V=Vitv 


est une fonction satisfaisant a (22), se réduisant a f(M) sur la surface 
et a a(M) pour t=0. 
Il est nécessaire évidemment que a(M) satisfasse a la condition (9), 
mais pour le cas discontinu les remarques précédentes s’appliquent. 
Pour démontrer que la solution est unique, supposons qu’il existe 
deux solutions distinctes V; et Vo. La fonction 


Vi = Vo =U 
1° Satisfait a l’équation 


(32) Au = R(M) 
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2° Se réduit a zéro sur S. 
3° Se réduit a zéro pour t = 0. 
Multiplions l’équation (32) par u et intégrons-la 


puisque R n’est pas fonction de t. 
D’ot, moyennant la formule de Green rappelée au n° 4, p. 119 


L(Y (2) +B) or 4d fate 


En intégrant la derniére équation par rapport a ¢, il vient 


i R(P) u? dwp = —F(t), 
D 
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ou F(t) est une fonction qui ne peut pas étre négative et qui se réduit 
a zéro pour t=0, puisque u se réduit a zéro pour t=0 et que R est 


une fonction positive. 
Donc u est identiquement nul. 


Cette méthode pour établir que la solution est unique s’applique 
a tous les probléemes relatifs 4 l’équation (22). Son application ne 
présente pas de difficulté et nous nous bornerons 4 donner le cas 


précédent (1). 


Nous avons enfin a traiter |’équation (23). Si nous y substituons 


cos 
V = y(z,y,2) X an 


sin 


nous aurons 


(33) Ay = -\ R(z,y, 2) 9. 


(‘) H.-B. Heywoop, Thése, p. 84. 
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Dans ce cas nous pourrons démontrer (n° 23, p. 151) qu’il existe 
un systéme de valeurs de 


Mi NOK eben ats 


qui sont toutes réelles et positives et un systeme de _ solutions 
correspondantes de (33) 


Pl, PQ, «+5 Pry -ees 


La solution pour les différents cas se fait d’une maniére un peu 
différente mais analogue a la précédente (voir le troisieme Chapitre). 

La résolution de l’équation générale elliptique dépend d’une fonc- 
tion de Green; et elle est semblable a celle de l’équation (21) : nous la 
réservons. 


CHAPITRE II 


L’EQUATION DE FREDHOLM 


1. — Nous avons dans ce qui précéde appris a ramener les 
problemes dont nous nous sommes occupés a l’équation de Fredholm 
qui s’écrit sous la forme 


b 
(1) nO ) K(s,t) y(t) dt = f(s) 


ou dans le cas de domaines a plusieurs dimensions sous la forme (1’) 
du n° 3, p. 4 


(V y(M) - i oy SOP) EP) dep = F000. 


Nous ne ferons, au moins jusqu’é nouvel ordre, aucune distinction 
entre ces deux équations et leur appliquerons le méme mode de calcul, 
le meme signe f{ représentant simplement dans la seconde d’entre elles 
une intégrale multiple. 

Dans chacun de ces cas la fonction y est a déterminer tandis 
que f, K sont des fonctions données, toutes ces fonctions étant bien 
entendu supposées intégrables. 

Nous indiquerons plusieurs méthodes pour résoudre l’équation (1). 


I. — METHODE D’ITERATION 


2. — Le premier procédé dont on ait disposé est di a Liouville (*). 
Il n’est autre, lorsqu’on lapplique au probleme de Dirichlet, que celui 
qui est bien connu sous le nom de méthode de Neumann (7). Neumann 


() Journal de Liouville, t. II, p. 24, 1837. 
(7) Voir Untersuchungen tiber das Logarithmische und Newtonsche Potential. 
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luicméme démontre la légitimité des opérations pour les frontiéres 
convexes et non «biétoilées» ('). En écrivant l’équation (1), sous la 
forme 


b 
(2) ols) =f K(s.1) v(t) dt + f(s) 


on est amené a rapprocher cette équation fonctionnelle de l’équation 
ordinaire : 
y = Fly) 
oul y est un nombre inconnu et F(y) une fonction connue de y. 
On sait que sous certaines conditions de continuité(?) on obtient 
la solution de cette équation en «réitérant» indéfiniment la fonction F. 
C’est-a-dire que les nombres 


91 =F(¥0), 92 =F(¢1), ---, Gn =F(Yn-1), --- 


ou yo est arbitraire sont des valeurs approchées qui tendent vers la 
racine de |’équation. 

Une méthode analogue peut s’appliquer a l’équation (2). 

Remarquons d’abord que s’il existe une solution ¢(s) continue et si 
K(s,t) est continue presque partout (n° 5's, p. 6), la fonction f(s) devra 
étre aussi continue. Nous supposerons donc que f(s) est continue. 

Substituons a y(t) dans le second membre de (2) la valeur 
que donnerait léquation si une fonction continue arbitrairement 
choisie yo(s) en était solution; on obtient 


b 
(3) els) =f K(s.t) yolt) dt + £(s) 


Opérons de méme en remplacant yp par vy; et ainsi de suite, on 
aura une suite de fonctions 


yo(s), y1(s), taey Yn(s), tee 


(') Voir la note de la page précédente. 
(7) Tl suffit que F(y) posséde une dérivée F’y et que l’on ait |F’y| < k, nombre 
fixe <1. 
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telles que : 


b 
(4) Oe | K(s,t) @n—1(t) dt + f(2). 


Si moyennant certaines hypothéses sur les données, nous pouvons 
montrer que yp(s) converge absolument et uniformément vers une 
fonction limite ¢y(s), lorsque n croit indéfiniment, on voit que la 
formule (4) deviendra a la limite la formule (2), c’est-a-dire que la 
fonction limite y(s) sera une solution de l’équation de Fredholm dont 
les fonctions y,(s) seront des valeurs aussi approchées que l’on veut. 

Or on a d’aprés (3), (4) : 


b b b 
(4’) onsi(s) = f(s) + f K(s,t) f(t) a+ | K(s.t) | K(t,t1) f(t1) dt, dt +... 


a 


b b b 
+f K(s,t) f Kit)... f K(tn—2; tn—1) f (tn—1) dtn—1 dtn—2 ... dty dt+Rn+1 
a a a 


avec 
b b b 
Ret = i. K(s.t) | K(t, f1) au K(tn—1, tn) vo(tn) din dtn_1...dty dt. 
a a a 


Nous avons supposé que K(s,t) est bornée dans son champ de 
variation. Soit M, la borne supérieure de sa valeur absolue; soient 
aussi P, N les bornes supérieures de |yo(s)| et |f(s)| dans (a,b), on aura : 


[Rn| < [M(b — a)]"P. 


D’autre part, on observe que y,, représente la somme de R, et des 
n premiers termes d’une série dont le terme général est, en valeur 
absolue, inférieur a 

[M(b — a)|"N. 


Dans le cas ott l’on a 
M(b—a) <1, 


(5/) 
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on voit que cette série est uniformément et absolument convergente 
dans (a,b) et que R, tend vers zéro avec =. 

Cette série étant indépendante de yo(s), on voit que la limite ¢(s) 
de y»(s) est indépendante de la fonction yo(s) dont on est parti. On voit 
aussi d’apres les formules (3), (4) que yo(s), étant supposée continue, 
yn(s) sera aussi continue et par suite aussi sa limite uniforme y(s). 
Ayant ainsi obtenu une solution continue y(s) de l’équation de 
Fredholm, il est facile de voir que dans nos hypotheses actuelles il ne 
peut y en avoir d’autres. En effet supposons que y(s) soit une nouvelle 
solution continue de (2). En prenant pour yo(s) la fonction 7(s), il est 
manifeste que tous les y,(s) seront identiques a 7(s). Par suite leur 
limite y(s) (indépendante du choix de yo) sera aussi égale a 7(s). 


Ainsi lorsque la borne supérieure M de |K(s,t)| est inférieure a = 

—a 

l’équation de Fredholm admet une seule solution continue ¢(s) qui est 
donnée par la série absolument et uniformément convergente : 


b 
p(s) = f(s) + f K(s,¢) f(t)dt +... 
(5) b i‘ b b 
+f K(s.t) [ K(tty)... f a ee ae ee ee 


Si au lieu de l’équation telle que nous venons de la traiter, on 
avait pris (comme dans I’Introduction, p. 1) l’équation 


b 
(bis) p(s) — A | K(s,t) p(t) dt = f(s) 


qui contient le paramétre \, ce qui revient a changer K(s,t) en \K(s,t), 
la solution aurait été évidemment 


b 
Omir rf K(s,t) f(t) dt +... 
wat [x60 PK) [Kea astn)atn dln a diy dt+.... 


De la résulte que la présente méthode consiste au fond a développer 
la solution suivant les puissances croissantes de X. 
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3. Les noyaux réitérés. — On écrit souvent cette série en intro- 
duisant les «noyaux réitérés». On appelle ainsi (et nous désignerons 
par les symboles Kj, Ko,...,Kn,.-.), les noyaux successifs définis par les 
relations 


b 
(6) Kj(s,t) = K(s,t),... Kn(s,t) = [ K(s,7) Kn—1(7,t) dr,.... 


Moyennant cette notation, les formules (3), (4), deviennent par un 
simple changement dans l’ordre d’intégration 
b 


b 
(7) onai(s) = f(o)+ [ Ki(syt) f(b) dt +. f Ku(s.t) FO) dt-+ Rost 


a a 


avec : 
Rn+1 = Kn+1(s, t) po(t) dt. 
a 


Dvailleurs, si ’on a encore M(b—a) <1, la série 
(8) — k(s,t) = K1(s,t) + Ko(s,t) +...+ Kn(s,t) +... 


est aussi absolument et uniformément convergente de sorte que la 
solution de l’équation de Fredholm peut s’écrire 


b 
(9) pare | k(s,t) f(t) dt. 


4. — Nous allons démontrer une propriété caractéristique de la 
fonction k(s,t). 
Pour cela, observons d’abord que l’on a évidemment 


b b 
(10) Kp(s, t) = | loans sh K(s, t1) K(t1, ta)... K(tn—1, t) dtn_1... dty. 


D’ot lon tire immédiatement par de simples changements dans 
Vordre d’intégration, la formule remarquable 


b 
(11) Kntn(st) = f Kn(s,7) Kp(z, t) dr. 


a 
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Alors on aura d’aprés (8) 


— k(s,t) — K(s,t) = Ko(s,t) + Ka(s,t) +...+Kn(s,t) +... 


b b 
= Ky(s,7) Ky(7,t)dr +... J Kn-als.t) Kulrt)dr +... 


b b 
=) Ki(s,7) Ki(7,t) dr +... J 16.7) kr artar +... 


d’ot: les deux formes : 


b 
—k(s,t) — K(s,t) = | [Ki(s,7) +... + Kn_1(s,7) +...] K1(7,t) dr 


a 


b 
= | Ki (s,7) [Ki(7,t) +... + Kn—1(r,t) + ...] dr. 


On obtient ainsi la formule caractéristique suivante 


b b 


(12) K(s,t) + h(s,¢) = f k(s,7) K(7,t) ar =f K(s,7) k(7, t) dr. 


a a 


5. Fonctions réciproques. — La formule (12) est démontrée pour 
M(b—a) <1. Mais elle a un sens méme en dehors de ce cas. Alors 
étant donnée la fonction K(s,t), nous appellerons réciproque de K(s,t), 
une fonction k(s,t) bornée et intégrable, telle que l’on ait : 


b b 
(13) K(s,t) + k(s,t) = ‘4 k(s,7) K(7,t) dr = | K(s, 7) k(7,t) dr. 
a a 
M. Volterra a montré que, connaissant une fonction réciproque k(s, t) 
de K(s,t), on peut toujours trouver une solution continue de |’équation 
de Fredholm (1). Soit, en effet, u(s) une telle solution, on aura 


“b 


(14) ult) = F()+ J K(tsta) w(t) at. 
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En multipliant par k(s,t) et intégrant, on aura 


[reo u(t) a= [roo soaes ff moo Kenyutyan dt 


égal d’apres (13) a 


b b 
/ k(s.t) rae | [K(s,t1) + k(s, t1)] w(t1) dt. 


On a donc 


b 


b 
o= | k(s.t) Fae f K(s,t1) u(t) dty 


de sorte que (14), ot l’on a remplacé t par s, devient : 


b 
(15) u(s) = #()— fP W(sst) FC a 


S’il y a une solution continue elle est done unique et donnée par 
cette formule. D’ailleurs pour prouver qu/il y a bien une solution 
continue, il suffit de montrer que le second membre de (15) satisfait 
bien a l’équation de Fredholm. En effet, si nous posons 


b 
(16) sts) = Fle) - i k(s,t) f(t) at 


la fonction donnée f(s) peut étre considérée comme une solution 
continue d’une nouvelle équation de Fredholm, |’équation (16), dont 
le noyau k(s,t) a pour fonction réciproque K(s,t). Alors, d’aprés ce 
qui précéde, f(t) est donnée par la formule (15), ot on substitue v(s), 
K(s,t), f(s) a f(s), k(s,t), u(s) respectivement ; ce qui donne 

b 


f(s) =v(s) - | K(s, ¢) u(t) dt 


a 


c’est-a-dire que v(t) est bien solution de l’équation de Fredholm donnée 


(1). 
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Dailleurs le méme raisonnement appliqué a la formule (13) 
considérée comme équation intégrale en k(s,t) prouve que si une 
fonction réciproque existe, elle est unique. 


5bis, — On peut considérer le second membre de (1) comme le 
résultat d’une opération fonctionnelle effectuée sur y, ceci signifiant 
qu’on en peut calculer la valeur lorsque y est donnée. Représentons 
cette opération par le symbole f = Sy. L’équation (9) exprime alors 
que l’opération qui transforme f(s) en 


b 
f(s) - | k(s,t) f(t) dt 


est linverse de S, c’est-a-dire qu’elle fournit » lorsque f est donnée. 
On peut la désigner par S~! et on a identiquement S~!S, =y 


6. La généralisation de Schmidt. — Nous avons vu que l’on peut 
former une fonction réciproque du noyau et par suite résoudre |’équation 


de Fredholm quand la borne supérieure de |K(s,t)| est inférieure a 5 
M. Schmidt a considérablement étendu ce résultat. 

1° Tl a d’abord montré qu’on peut appliquer la méthode d’itération 
lorsque l’on a 


(17) [ [iso0P dans, 


ww désignant une quantité numérique déterminée. Cette inégalité a lieu 


—a 


1 
quand |K| a une borne supérieure plus petite que — mais aussi dans 
—a 
des cas plus généraux. 
Cette démonstration est fondée sur la formule (11) et sur V’inégalité de 
Schwarz (8), page 8, d’ou l’on tire 


b b 
(18) Kn4i(s,t))? < / [Kn(s,r)]? dr. i [K(r, t)]° dr, 


et par suite 


(19) ae Kn+i(s, t)] * dsdt < [ [te (s,r) Paras. ff pw (r, t)| * dr dt. 
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On déduit facilement de (19) 


(20) ff tne “(| dette Lf ft (s, 2) yPasat] 


Or on tire de (11) 


b pb 
Kn(s.t) = | | K(s,71) Kn—2(11, 72) K(ro, t) dri dra, 


et en se servant de l’inégalité de Schwarz (8), page 8, formée comme (8), 
mais avec des intégrales doubles, on aura 


< [ 4 [Kn—a(r1,72)]° dr, dra. [ i [K(s,71)]°[K(r2, t)]? dra dra. 


En combinant avec (20), on a donc 


n—2 
5) 
K,,(s, t)| <{f [ts (s, t)| ava ” 


En utilisant maintenant I’hypothése (17), on voit que les termes de la 
série (8) sont a partir du troisiéme rang inférieurs en valeurs absolues a 
ceux de la série convergente 


L 
2 


. ie [K(s, r)|? dry ; i [K(r2, t)]? in| 


NIr 


[fistnPar. [isco] ar AWwetGey they hth). 


Alors la série (8) est absolument et uniformément convergente, sa 
somme —k donne la fonction réciproque k, et par suite on peut écrire la 
solution unique (9). 


2° D’autre part, i est facile de résoudre directement l’équation de 
q=p 


Fredholm dans le cas ou son noyau est de la forme S/ ag(s s) By(t). 


q=l 
Pour faciliter la discussion, nous le supposerons multiplié par un facteur 


numérique arbitraire 4. Considérons donc |’équation 


(21) os) —A f “ott bs ag(s) y(t) at = 6). 
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On voit tout de suite quelle est la forme de la solution en écrivant 
Véquation sous la forme : 


P b 

(22) ols) = f(s) + aa(s) | f ante) o(0 at] 
1 a 

c’est-a-dire 

(23) soz )- a Devan) 


ou les ag sont des constantes a déterminer. On les calculera par la méthode 
des coefficients indéterminés en substituant (23) dans (21) ce qui donne 


P PSP b b 
(24) Sole) {aSo f a(t) Bg(t) at — aq — | Ba (t) f(t) a = 0, 
1 r=1 a a 


de sorte qu’on obtiendra une solution du probleme en déterminant les 
constantes a, par les équations linéaires 


(25) Yor [a font a0 a] — 0 = [ato 
G79 um 
Si nous posons : 
Ang = fF alt) Bult) dt 


on voit que le déterminant des coefficients de ces équations sera un 
polynome en X de degré p au plus 


AA —1 AAi2 sa’ AAip 
(26) Ag AAg2 —-1 ... AA2p 
ALG We sate NN com 


Ce polynome n’est certainement pas identiquement nul puisque, pour 
A = 0, il est égal a (—1)?. On pourra donc résoudre les équations linéaires 
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quand A a une valeur quelconque autre que les racines (en nombre au plus 
égal & p) du déterminant. 

Et, si nous supposons que les fonctions a,(s) sont linéairement 
indépendantes, nous voyons méme que les équations (24) ne peuvent avoir 
lieu que si les équations linéaires (25) sont vérifiées. 

En résumé si les fonctions ag(s) sont linéairement indépendantes 
Véquation fonctionnelle (21) admet une solution unique y(s) que nous 
savons déterminer complétement sous la forme (23) et ceci pour toute 
valeur de » sauf peut-étre pour p valeurs de A au plus, racines du 
déterminant. 

Nous remarquerons que ces racines ne dépendent que du noyau, nous 
les retrouverons plus tard dans le cas général (p. 57) sous le nom de 
constantes caractéristiques. 

3° Revenons maintenant 4 la méthode de Schmidt qui résulte de la 
combinaison de 1° et 2°. 

Ici encore nous faciliterons la discussion en supposant que le noyau 
posséde en facteur une constante arbitraire A. Le noyau étant désigné 
par AK(s,t), o& K est une fonction continue quelconque, on voit que la 
condition (17) ne sera vérifiée que pour des valeurs de |A| inférieures 4 une 
certaine quantité finie. Pour résoudre l’équation pour les autres valeurs 
de A, nous montrerons d’abord qu’on peut toujours trouver des fonctions 
continues a,(s) linéairement indépendantes et des fonctions continues /3,(t) 
en nombre égal de facon que l’expression 


2 


ff |x (s,t) — Yeats ) Bg(t)| ds dt, 


soit aussi petite que l’on veut. 


: : . : : <s<b 

Le noyau étant uniformément continu dans le domaine ie s 
a<t<b 

on peut trouver un entier p tel que deux valeurs de K(s,t) prises en 


2(b— a) 


deux points a une distance inférieure a , he puissent différer de 


plus de 6. Divisons maintenant l’intervalle (a,b) en p intervalles égaux 
71,12,.-.,%p et appelons leurs milieux $1, 52,..., 8p. 
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Nous prendrons pour les {,(t) les fonctions 
Bq(t) = K(sq,t). 


Pour définir a,(s), nous construirons une fonction a,(s) égale 4 1 dans 


Vintervalle ig et a 0 en dehors de ig; puis pour la rendre continue, nous 


modifierons légérement ici la construction de Schmidt. Appelons 7, un 


2(b — a) : 


intervalle de longueur et ayant pour milieu l’extrémité commune a 


aes : ; _ r a 
ig et ig41; nous prendrons a4(s) = aq(s), sauf dans les intervalles i,_, et 7% 


ot nous prendrons pour a,(s) une fonction linéaire de la valeur 0 a la 
valeur 1 ou inversement. 


Alors quand s est dans i, sans étre dans 7,_, ni dans 7, 


ge Ol Bi 2 


|K(s,t) — S © ag(s) Bq(t)| = |K(s, t) — K(sq,#)| < 6. 


/ 


Quand s est dans i, on a: 


P 


|K(s,t) 7 eG) Bq(t)| = |K(s, 4) — Qq(8) K(Sq,t) — Og41(8) K(sq41,)], 
1 


et comme on voit facilement que dans 7, on a a,(s) + aq4i(s) = 1, la 
quantité précédente sera égale a 


|ag(s)[K(s, t) — K(sq, t)] + ag41(s)[K(s, t) — K(sq41, 6)] | 
< |K(s,t) — K(sq,t)| + |K(s, t) — K(sq41,t)| < 26. 


Ainsi, on a en tout cas, 


2 


ae K(s,t)— J ‘ae(s) Aolt)| dsdt < 46°(b — a)?, 


et en faisant croitre p nous pourrons rendre le second membre aussi petit 
que l’on veut. 
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4° Ceci étant, soit l’é6quation de Fredholm donnée : 


b 
(27) y(s) — A j. K(s,t) p(t) dt = f(s), 


choisissons les fonctions a, 3, de facon qu’on ait par exemple 


(28) el 


Ecrivons alors (27) sous la forme 


p 2: 
K(s,t) — > ots) dsdt < > 


b 
(29) ols) -A f H(s,4) elt) dt = fils), 


en posant 


H(s,t) = 


K(6,t) — Srag(s) at) , 

(30) ; : ; 

f(s) = f(s) tA- / ovg(s) By(t) y(t) at 
1 a 


Alors d’aprés 1° et (28) il existera une fonction h(s,t, A) réciproque du 
noyau AH(s,t) et l’équation de Fredholm (29) admettra une seule solution 
continue 5 

els) = fils) — [W(s,t,9) Aa 
a 


Si nous y remplacgons f; par son expression (30), nous voyons que 
Véquation (27) est équivalente a l’équation 


P b b 
(8) = £0) + X92 fags) Bult) ot) at — f(s, t,d) #18 at 
“Tl a a 


B b fb 
— x h(s,r, A) ag(r) Bq(t) p(t) dr dt, 
are 
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qui peut s’écrire : 


b P 
(32) y(s)—A | S[fals,) Balt) o(6)] at = F(s), 


1 


en posant 


b 
F(s) = f(s) — | h(s,t, A) f(t) de, 
b 
Fels, A) = Gel 9) = |. Wihsyr Abarat 


Alors on est ramené & une équation de Fredholm (32) de la forme (21) 
déja examinée. La seule différence, c’est que le déterminant (26) n’est plus 
un polynome en 4X, les fonctions fy dépendant déja elles-mémes de A. Mais 
il est encore égal & (—1)? pour \ = 0 et par suite non identiquement nul. 

Ainsi lorsque le noyau AK(s,t) est un noyau continu quelconque et 
lorsque X nest pas égal a certaines valeurs exceptionnelles (racines du 
déterminant mentionné) l’équation de Fredholm a une solution continue 
unique. 

5° Le cas ot: la fonction inconnue dépend de plusieurs variables tel 
qu'il a été indiqué au n° 8, p. 4, se traite par la méme méthode. 


Tt. Equation de Volterra. — On appelle équation de Volterra 
Véquation 


(33) os) - f "K(s,t) o(t) dt = f(s) 


qui dérive de (1) en remplagant la limite supérieure fixe b, par la limite 
variable s. 

Cette équation est un cas particulier de l’équation (1), mais la 
convergence des séries obtenues est plus favorable. Ces séries sont 
uniformément convergentes toutes les fois que K(s,t) est borné et 
intégrable. 
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On peut, en effet, identifier ’équation (33) avec l’équation (1) en 
prenant dans l’équation (1) a la place de K(s,t) une fonction qui est 
égale a K(s,t) tant que t < s, mais qui s’évanouit pour ¢t > s. La solution 
s’exprime par la série — analogue a la série (5), p. 40 — 


p(s) = f(s) + [K(s rears f [* (s, $1) K(s1,t) f(t) ds, dt+. 
(34) +f he 4 hs - #1) K(oiy 80) Klsosiia sc Kya) SO) dar dep db bes 
= s+ [K(s, t) + Ko(s, t) + K3(s,t) +... + Kn4i(s,t) +...] f(t) dt, 
ou 


s S] Sn-1 
Kn(s,t) = | ou, K(s, 81) K(s1, 82)... K(sn_1, t) ds1 dsq...dsy_1. 
a a a 


Soit M la borne supérieure de la valeur absolue de K. On a 
|K (s,t)| <M 


|Ka(s, t)| = a K(s, 81) K(s1, t) ds1 


<[ M? ds; = M*(s — a) 
a 


Ss 
|K3(s, t)| = i K(s, 81) Ko(s1,t) dsy 


(s—a)? 
19 


a M3(s — a) ds = M? 
a 


et en général 


|Kn(s, t)| < ee 


Donc la série 


(35) K(s,t) + Ko(s,t) + K3(s,t) +...+Kn(s,t) +... 
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est comparable avec 


ie ; | im Gag 
ee (n — 1)! 


M+M?(s a) +M3 4 


de sorte que la série (35) est toujours uniformément et absolument 
convergente, et la solution (34) est toujours valable. 

M. Volterra a beaucoup étendu ce résultat, il a traité d’une manicre 
analogue l’équation de premiére espéce 


[ o.petoat =f 


et il en a fait des applications. Nous renvoyons le lecteur a son 
mémoire dans les Annali di Matematica, 1906. 


II. — METHODE DE FREDHOLM 


8. — La méthode de Fredholm pour résoudre |’équation intégrale 


b 
(1) y(s) - | K(s,t) g(t) dt = f(s) 


consiste a étendre au cas des équations fonctionnelles le procédé de 
résolution des équations du premier degré a une infinité d’inconnues. 
Sa méthode sera développée un peu plus loin sous forme 
synthétique, mais elle paraitrait tout a fait artificielle si nous ne 
montrions auparavant comment on y est conduit par un passage a la 
limite. 
Divisons l’intervalle (a,b) en intervalles partiels par les points 


hes 
89 =a, 85 =at+6, s2.=a+20, ..., Ss, =atnd=b, (5- =) 
nm 
et prenons d’abord comme inconnues non pas la fonction ¢y(s), mais 
ses valeurs 


Lp = 9(Sp) (p =0,1,...,n). 
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On aura d’aprés (1) 


b 
(36 ty — f Klsp.t) oll) dt = Hsp) 


ou, d’aprés la définition de l’intégrale 


q=n 
tp — jim f > tqK (Sp, wo} = f (Sp). 


q=1 


La méthode consiste alors & considérer comme valeurs approchées 
des x» les solutions des équations linéaires : 


(36°) Lp — 6 > LqK(Sp, Sq) = f(Sp) 


(p =0,1,2,...,n). 


On voit alors que l’existence de la solution dépend de la valeur du 
déterminant des coefficients des zp 


1—6K(s1,81) —6K(s1,82) ...  —dK(s1, Sn) 
— 06K(s2,81) 1—d6K(s9,s2) ... — dK(s2, $n) 
Dn = 

— dK(sn,81) —O6K(sn,82) ... 1—dK(sp, 8p) 

pFn 

=1-6 > K(Sp, Sp) 
pal 
K(s1,81) ... K(s1,8n) 


K(sp, 8p) K(Sp, 8) 


+...+(-1)" 6" 
K(sq,8p) K(sq, Sq) 


1 
ae Ds 


K(sn,$1) --- K(8n,8n) 
Si on fait croitre n indéfiniment, on voit que chacun des termes de 


cette somme a une limite déterminée. De sorte que D, tend, au moins 
formellement vers la série 


p=1- [Ken wrap f 


K(71,71) K(7,72) 


dt, dTg+.... 
K(72,71) K(72,72) 
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Il y aurait lieu ensuite de voir aussi ce que deviennent les formules 
de Cramer, de rechercher ce qui se passe quand D = 0 et surtout de 
rendre la méthode rigoureuse. C’est ce qui a été fait par M. Hilbert (*) 
qui retrouve ainsi les formules que nous allons établir par la méthode 
synthétique de Fredholm. 

Nous nous contenterons cependant d’avoir montré, au moins sur 
Vexemple de la série D, que certains développements compliqués 
introduits a priori dans la méthode de Fredholm (n° 10, formules 
(40), (41), ...) ne sont pas construits de toutes piéces, mais qu’au 
contraire leur considération s’impose quand on s’inspire du passage a 
la limite. 


9. Théoréme de M. Hadamard (7). — Pour étudier la convergence 
de D, il serait utile de connaitre le maximum de la valeur absolue des 
déterminants qui y figurent. M. Fredholm s’est servi dans ce but d’un 
théoreme de M. Hadamard. Nous nous contenterons de démontrer ce 
théoréme dans le cas d’un déterminant 4 termes réels (3) 


ai ai2 «++ Alin 
A= 


Qnl OQn2 +--+ ann 
Il s’agit de prouver que l’on a : 
(87). AA. < (a? Hate. pat, Me ass baad, cla boty be he?) 
Pour cela, posons 


2 2 2 2 
Oy = Gpy Ogg bene + Gap (p =1,2...n) 


ee Gotting. Nachrichten, 1904. 

eg Bull. des Sc. Math. M. Muir, en réponse a Lord Kelvin (qui possédait 
dés 1885, l’énoncé de la proposition en question), lui avait communiqué, en 
1886, un raisonnement qui démontrait cette proposition ou du moins la rendait 
trés vraisemblable. Ce raisonnement a été publié en 1901 dans les Reprints de 
V Educational Times. (Voir aussi l’article récent de l’auteur dans les Trans. of 
the Royal Soc. of South Africa, t. Tl, 2™° partie) (Note de M. Hadamard). 

(3) La démonstration ne serait d’ailleurs guére plus compliquée dans le cas des 
termes imaginaires. 
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et cherchons la borne supérieure de A? lorsque les ap, varient de fagon 
que les 5? gardent des valeurs arbitraires données. Il est manifeste 
que A? est une fonction continue et dérivable par rapport aux ap¢q qui 
eux-mémes restent bornés quand les 6? sont fixes. Donc A? atteint un 
maximum A? donné par les régles du calcul différentiel. 

On aura ainsi : 


OAo fis 
Oapn He 


= ie ORG Ag 
0 =dAo = Sa Oa 
pour toutes les variations telles que 62 = constante, c’est-a-dire telles 


que 
; 57, = dp dap, + Ap2 dapy +... + Ayn dan = 0. 
On a donc ado 
Dang ~ PAP (p =1,2,...,n) 


ou les \, sont des constantes convenables. 


D’ailleurs wo 
Jang 
dans le développement de Ap. D’ot : 


est, comme on sait, égal au coefficient Apg de apg 


Ang = ApGpq (p = 1,2,...,n) 


et 
Ao = Ap1ap1 aes Apn&pn = Apde, 
ce qui donne 


Ao 


Ang = Fe. Apq: 
Dp 


Alors le déterminant des Ap, (c’est-a-dire le déterminant adjoint 
de Ao, égal comme on sait a Aj ') sera évidemment égal au 
déterminant Ap des ap, multiplié par 


AG 
6253... 62 
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D’ot : 
AS, 
&... 62° 


n 


n-1 _ 
Ay = 


avec |A| < Ao. 

D’ot enfin Pinégalité annoncée A? < 67...62. 

En particulier, si M est un nombre supérieur ou égal & chacune des 
valeurs absolues de tous les termes de A, on aura 62 <nM?, d’ot 


(38) |A| < Vn” (1). 


10. — Nous sommes maintenant en mesure de développer la 
méthode de Fredholm. 

Pour faciliter la discussion nous introduirons comme précédemment 
(n° 2) un paramétre numérique arbitraire \ et nous nous proposerons 
de résoudre l’équation 


b 
(1bis) y(s) — [ K(s, t) p(t) at = f(s), 


ou f(s) est continue dans (a,b) et ott K(s,t) est une fonction continue 
presque partout (n° 5>#5, p. 6). 
Guidés par la méthode de passage a la limite dont nous avons 


(') M. Kneser a fait connaitre récemment une démonstration intéressante de ce 
méme théoreme (Die Integralgleichung, Vieweg, 1911, p. 227). Une démonstration 
analogue a été donnée en méme temps par M. Tommaso Boggio, Bull. des 
Sc. Math., 1911. Elle consiste 4 réduire la proposition au cas (ot elle devient 
évidente) ot! tous les termes d’un c6té de la diagonale principale sont nuls. On 
y arrive en opérant une méme substitution orthogonale sur chacune des lignes 
du déterminant A. 
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donné une idée plus haut (n° 8), nous envisagerons les deux séries 


b 
S $, $8] 
(40) Ne = ok (S ") asi +. 
yn b b 
ap [« & sa) ay dip. don + 
! a t, S51, 82,...,8n 
b 
(41) = | K(s1, 51) ds, +.. 
ae ye b 
La eb ee 
S a a $1, 82,+--,8n 
avec 
ar . K(s1,t1) K(s1,te) ... K(s1,tn) 
(42) Oe ae ee ne et aa eS 
ty,t2,---,tn 
K(sn,t1) K(sn,tg) ... K(Sn,tn) 


[La quantité D(A) est dite le déterminant du noyau K(s,t).] 

Puisque |K(s,t)| a une borne supérieure finie M dans son champ de 
variation, ces deux séries seront convergentes quel que soit A. En effet 
d’apres le théoreme de M. Hadamard, on a 


K $1, 52,---,8 
ty, ¢9,... ‘ty 


quand les variables restent dans (a,b). Donc les termes de D (>) et 


< vn" M", 


de D(\) restent en valeurs absolues respectivement inférieurs & ceux de 


AI" 


n! 


1+ pim@o—a)+...4+ 0 vir Moa)” 


M+ 2\A|(6—a)M?4+...4 (n+ 1)"+1 (b— a)" M™t1 + 


Le rapport d’un terme au précédent est pour la premiere série de 
la forme 


Ia (1, 1) ee a)M 


n 
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et pour la seconde 


nN (+2) ae a) M; 


dans les deux cas il tend vers zéro quel que soit \. Ainsi les deux 
séries représentent deux fonctions entiéres de 4. 

Les zéros de la fonction D(A) sont déterminés par la fonction K(s, t) 
seule. On les appelle les constantes caractéristiques du noyau K(s,t). 
On sait (n°5, p. 6) qu’elles sont isolées les unes des autres. Lorsque 
A West pas égal & une constante caractéristique, le quotient 


BUA) 
4 K = 
(43) (61) = i 
représente une fonction de s, ¢ bien déterminée. Quand )_ varie, 
c’est une fonction méromorphe de \ qu’on appelle la résolvante du 
noyau K(s,t). 


11. — Nous allons maintenant chercher & montrer que la fonction 
k(s,t) = —AK(s,t,A) est la fonction réciproque (n° 5, p. 42) de AK(s,t) 
lorsque X nest pas une constante caractéristique. Ce sera suffisant 
pour prouver que dans ce cas, l’équation de Fredholm a une solution 
continue unique donnée par la formule (9), p. 40. 

Etudions d’abord le genre de discontinuité de D () ) par rapport 
aux s et ¢. On a en développant (42) par rapport a la premiére ligne 


eee Keon) +An 


t, $1, $2,---,8n §1,52,.--,5n 
avec 
0 K(s, 81) K(s, sg) ... K(s, sn) 
8,81,---,8n\ _ |K(s1,¢) K(si,s1) K(s1,82) --- K(s1, $n) 
An=A = 
t,81,-+-,8n 


K(sn,t) K(sn,s1) K(sn,s2) .-. K(Sn, Sn) 


L’EQUATION DE FREDHOLM 59 


En portant cette expression dans la formule (40), on a 
(44) D G ) = D(A) K(s, t) + R(s,t, \) 


avec 


b (yy? b b 
R(s,t,4) = -a f Ards +...4 5° f | An ds, ...d8y + 
a © a a 


En développant A, par rapport a la premiére ligne, on voit que 
le terme général de R est une somme de termes de la forme (6) 
considérée dans l’Introduction, p. 6, puisque K(s,t) est convergente 
presque partout; donc R est une somme de termes continus. Et 
comme la série R est aussi uniformément convergente on voit qu’en 
définitive R(s,t,) est une fonction continue quand s, ¢ varient dans 
a, b. 

Des lors, la formule (44) montre que D c A) est une fonction qui a 
les mémes points de discontinuité que K(s,t). 

D’autre part on a: 


ie [andor ane fi [tc 1)? K(s, sp) x 


K(s1,t) K(s1,81) ... K(s1,8p-1) K(s1,8p41) ... K(81,8n) 
K(sp,t) K(sp,81) -.. K(Sp,Sp-1) K(Sp,8p41)  «-. K(sp, 8n) | ds1-.. dsp 
K(sn,t) K(sn,s1) ... K(S8n,Sp-1) K(sn,8p41) ... K(Sn, Sn) 

Pp=N Ab b 

--y ff K(s, Sp) x 

p=1 70 a 
K(sp,t) K(sp,81) ... K(8p,Sp-1) K(S8p,5p41) ... K(8p, Sn) 
K(s1,t) K(s1, 81) nent K(s1, 8p—1) K(s1, 8p41) gine K(s1, 8n) ds ,...d8n 


KiGSAst) ae oe Re ck A Se ee oe Boe URS an) 
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en remontant la p’@”¢ ligne jusqu’A la premiére place. Si ensuite on 
remplace sp, 8p41,---,5n, Pal 7,Sp,---,8n—1 Ce qui ne change rien dans 
la valeur de lintégrale, on voit que le signe > porte sur n intégrales 
égales a 


K(7, ft) K(7, 81) ... K(7, 8n-1) 
b b 
K t K .. K ok 
if a, K(s,7) Bint) (Sty 81) (S1, 82-1) dr ds,...d8y—1 
a a a 
K(S,_1, ¢) K(Sn—1, $1) see K(Sn—1, $n—1) 
b b 
=} ah K(sir) Km) K (S08) de dst dina 
a a 


Bij xx Sami 
b b T,S1,+-+55n—-1 
+f soe K(s,7)A drds,...d8py—4 
a a t, $1,--+;8n—1 
En sorte que 


R(s,t, A) = 
b N=0CO 7 y\\n b b 

-| K(s,7) K(r,t) dr x asa ey i (8 doy dona 
a a a 


Sigicg Sl 
b = 22 —))n b b 

-{ K(s,7) ey n a a (Pe) day dona dt 
a , n) a a t,$1,---,Sn—-1 


K(s, 7) K(z, t) én x AD(A) + K(s,7) «A. R(z,#, A) dr. 


a 


R(s, t, 4) = sf Ks) {K(r, t) D(A) + R(z, t, A) } dr. 


En tenant compte de (44) on a enfin 


D(FA) — D(A) K(s,t) af 


a 


K(s,r)D @ A) dr. 
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En procédant avec les colonnes des déterminants comme nous 
avons fait avec les lignes, on obtiendrait de la méme maniére pour 
valeur du second membre 


ao C ) K(r, 8) dr. 


On a ainsi Videntité fondamentale 


(45) p(*a) - DENK) =A fK(s,7)D (7 X) dr 


= afr e A) K(r,t) dr. 


En particulier si \ n’est pas une constante caractéristique, D(A) n’est 
pas nul, la formule (43) a un sens et on a 


b 
(46) RG) = koa = af KG KGa 
= af Ko X) K(7, t) dr 


qu’on peut écrire sous la forme 


b 
k(s,t) + [AK(s,#)] = i) [AK(s, 7)] k(r,t) dr 


=| k(s, 7) [AK(r, t)| dr 


ce qui démontre la proposition. 

En résumé nous voyons que (d’aprés un raisonnement déja em- 
ployé au n°5, p. 42) ’équation de Fredholm (1) ott K(s,t) est 
une fonction continue presque partout (n° 5, p. 6), ott f(s) est 
une fonction continue et ott \ n’est pas une constante ca- 
ractéristique (n° 10, p. 57), admet une solution continue unique 
donnée par 


b 
i (47) (s) = f(s) +a K(s,t,A) f(t) dt, 
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ot. K(s,t,A) est une fonction méromorphe de \ donnée par les 
formules (43), (40), (41). 


12. — On voit l’avantage de la méthode précédente sur la méthode 
d’itération : c’est que celle-ci ne réussissait que si, M étant la borne 
supérieure de |K(s,t)|, on avait (n° 2, p. 39) 

|A|M(b—a) <1 
c’est-a-dire : 
M(b—a) 

Il est vrai qu’on peut traiter les autres cas au moyen de la 
généralisation de Schmidt que nous avons exposée plus haut, p. 44. 
Mais celle-ci ne fournit pas une formule générale unique comme la 
précédente et ne méne pas a une discussion aussi simple. Elle ne 
fournit pas une expression explicite du résultat, dans laquelle chaque 
terme soit directement donné par une expression analytique connue, 
ce qui est en général trés important au point de vue des calculs 
théoriques. 

Au point de vue du calcul numérique, au contraire, ceci constitue 
plutot un désavantage de la méthode de Fredholm, puisque le calcul 
d’un terme ne sert en rien pour celui des suivants. 

Mais ce désavantage est compensé par un avantage essentiel, celui 
d’une convergence trés rapide, lequel est précisément en rapport avec 
ce fait que nos séries actuelles ont en A, un rayon de convergence 
infini et les premieres (5), (5’) un rayon de convergence fini. Les 
séries (40), (41) convergent a la facon des développements de fonctions 
entieres, c’est-a-dire comme des progressions géométriques a raison 
infiniment petite a partir d’un rang assez élevé. Dans beaucoup de 
cas, les deux ou trois premiers termes suffiront, dans chacune d’elles, 
dés que A ne sera pas trop grand. 

La série (5’) représente, comme nous l’avons remarqué, le dével- 
oppement de ¢y(s) suivant les puissances de , c’est-a-dire le quotient 
des deux séries (40), (41). Son rayon de convergence est égal au module 


|Al < 
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du premier pdle (en \) de ¢(s), c’est-a-dire (voir plus loin, n° 19) 
du premier zéro A, de D(\). Elle converge comme une progression 


géométrique de raison +]; par conséquent, comme une progression 
1 
de raison a s'il s’agit de la série (5). 
1 

Il faut enfin remarquer que, si l’on doit calculer les constantes 
caractéristiques (1) (n° 10)—ce qui est indispensable dans les appli- 
cations mentionnées au Chapitre I, — on ne peut employer notre 
premiere méthode qu’avec la généralisation de Schmidt, et méme en 
prenant l’entier p du n°6 d’autant plus grand qu’on veut calculer des 
constantes \,, d’indice plus élevé. 


13. Noyaux infinis. — Nous n’avons pas encore considéré le cas 
ou le noyau AK(s,t) tout en restant intégrable deviendrait infini en 
certains points. Fredholm a aussi indiqué pour ce cas une méthode 
qui revient a substituer a l’équation donnée une équation équivalente 
ot le noyau est fini. 

Revenons en effet a la formule (7), p. 40, en y remplacgant yo(t) par 
une solution y(s) de l’équation (1); alors y, = v2 =...¢n41 = ct 
nous aurons en remplagant n par n—1 et K(s,t) par \K(s, t) 


b 
(48) y(s) — a” i} Kn(s,t) o(@) dt = fale) 


ou fr(s) est une fonction continue connue 


b b 
(49) als) = (5) +2 f Kils,t) fat +... bart f Ot On 


On voit que y(s) est une solution de l’équation de Fredholm (48) ot 
AK, f sont remplacés par \"Ky, fn. Inversement, soit y(s) une solution 
de l’équation (48). Si dans les relations (3), (4), on remplace yo par , 
(K étant toujours changé en AK), yp sera aussi égal a y. Il suffit alors 


(') Les constantes en question sont évidemment (voir plus loin, n° 19) les pédles 


de K(s,t, A). 
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d’ajouter membre 4 membre les n premiéres de ces équations (3), (4) 
pour constater que la fonction 


_ Yor Pr +---+Pn-1 
n 


wv 


est une solution de (18). C’est d’ailleurs aussi évidemment une 
solution de (48) comme y, donc VW = y et gy est bien une solution 
de (bis). Les équations (1's) et (48) sont donc équivalentes. 

[En vue des applications numériques, il y a lieu d’observer que 
Von obtient la transformée en ” de l’équation D(A) = 0 relative au 
noyau K(s,t) en égalant a zéro le «déterminant» du noyau n fois 
réitéré Kn(s,t).| 

Or, il arrive dans un grand nombre d’applications que le 
noyau K(s,t) étant infini en certains points, on peut trouver une 
valeur de n pour laquelle le noyau réitéré K, est borné. On sera donc 
ramené au cas précédent. 


14. — La discussion des cas — qui se rencontrent fréquemment en 
Physique—ot cette condition se trouve réalisée sera, contrairement a 
ce qui avait lieu dans tout ce qui précéde, différente suivant le nombre 
des variables. 

Soit d’abord une intégrale simple. Nous supposerons que K(s,t) n’est 


infint que pour s=t et comme pode’ pour que K soit intégrable, il 


faut d’ailleurs que a qui est > 0, soit <1. Pour préciser, supposons 


K(s,t) = al 


ou H(s,t) est une fonction bornée intégrable. On a 
b 
Ko(s,t) = | K(s,7) K(z,t) dr. 
a 


Si M est la borne supérieure de |H(s,t)|, on aura : 


a aes | 1 
249, S : 
|Ko(s,t)| <M dr 
o PoP FH 
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Or posons 


¢=sy+t(l—y), 


Vintégrale deviendra 


Joie ——_ 
Gee alt Salt 


Donc : 


dy = aa 
|Ka(s, t)| < M? i wet — yl x |s t|! 2a = Pals — tl! 2a 


ou Py est un nombre positif fixe. 
De méme, comme on a 


b 
K3(s, t) = Ko(s, 7) K(z, t) dr 
a 
on aura 
|K3(s,7)| < Pals — t/?- 3% 


et en général 
|[Kn(s,t)| < Pnls— fie tae, 


, F ‘ pe ‘ 1 
Soit alors n le premier entier supérieur a j (> 1); on aura 
Qt 


n—-l—-na> 0. 


Done P,|s —¢t|"-!-" sera borné (et < Pp(b-a)"-1-"2). 
De sorte que le noyau de l’équation (48) sera fini. 
En particulier, si a < 5 il sufhra de prendre n = 2. 


K 1 . , 
15. — Dans ce meme cas a < es M. Hilbert a donné un moyen 


élégant de traiter directement léquation (1>) en modifiant la méthode 
de Fredholm pour le cas d’un noyau fini. 
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Il s’agit de tourner la difficulté provenant de ce que les termes des 
diagonales principales des déterminants 


8).3] pe ee 
K( Fl, ? ") 
t, 81,..-,8n 
sont infinis sauf le premier comme étant égaux a 
K(s1, 81),---,;K(Sn, $n). 


Pour cela, M. Hilbert remplace la fonction donnée K(s,t) par une 
fonction Ko(s,t) égale partout a K(s,t) sauf pour s =¢ ot l’on prend 
Ko(s,s) = 0. Il est manifeste que la nouvelle équation intégrale obtenue 
sera entiérement équivalente a la premiere, car une intégrale simple ne 
change pas de valeur si l’on change la valeur de la fonction a intégrer 
en un seul point. 

En opérant ainsi, les nouvelles expressions de D (F>) et de D(\) 
ne contiendront que des termes finis. I] reste 4 montrer que quoique 
Ko(s,t) ne soit pas bornée, ces deux séries sont convergentes. 

En effet, considérons par exemple ce que devient 


ei), 
81,---38n 


En supposant d’abord s1 > sg > 83... > sn, nous multiplierons la 
premiére ligne de ce déterminant (voir formule 42), par (s; — s2)%, la 
seconde par [(s, — s2)~% + (s2 — s3)~°]7!..., la derniére par (sn—1 — sn). 
Les éléments du déterminant obtenu restent tous inférieurs, en valeur 
absolue, 4 un nombre fixe P. D’aprés le théoreme de M. Hadamard 
(n° 9, p. 55), on a donc 


i Ga a) < (81 — 82)” [(s1 — 52)~% + (s2 —83)"%] >... 


81, 82,--+,8n 
(Sn—1 — 8n)~»Vn™. P™ 


dans le domaine D défini par les inégalités b > s1 > sg... > 8, ><a. 
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Or, les intégrales de |K| dans chacun des n! domaines analogues 
a D obtenus en permutant les s, sont égales. Le terme général u, de 
la série (41) qui représente le «déterminant» D(A) relatif a Ko est donc 
au plus égal a 


; |A”| vn P” +++ f (81 — 8) [(s1 — 82) * + (s2 — 83) *] 
(49Pis) D 


fer —a7—l 
[(Sn-2 — 8n—1) ~~ + (Sn—1 — $n) °] (Sn—1 — $n)% ds1...d8n. 


Un changement de variables tel que s; = a+ (b—a)s; ramene cette 
expression a la forme 


||" (/n)” P"(b _ aja)" i 


ot I, représente Vintégrale qui figure dans (49s), mais étendue au 
domaine 1 > 81 > sz... > 8» > 0. Le calcul de I, montre (!) que l’on peut 
n 


eerie Lx Te ou A est un nombre fixe. De sorte qu’on a pour le 
7nd 


terme général up de D(A) 


1 
Vlun| < |A|P(b—a)!-° A x = 


n2 > 


1 F si 
Comme a < 5 %/|un| tend vers zéro et la série D(A) est absolument 


convergente. M. Poincaré(?) est parvenu a étendre le procédé de 
M. Hilbert au cas ot a est un nombre quelconque compris entre 
0 et 1. Il arrive au résultat suivant : Soit n un entier quelconque 


On peut encore appliquer les formules de Fredholm 


supérieur a i 2 = 
(40), (41), (43), (47) comme si le noyau K(s,t) était fini pourvu qu’on 
ait soin de supprimer dans le développement de chaque déterminant 


(*) Voir Hilbert, Grundztige ..., Géttinger Nachrichten (Erste Mitteilung), 
1904, Heft 1, p. 84. 

(7) Remarques diverses sur Véquation de Fredholm. Acta Mathematica, t. III, 
1910, p. 73. 
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$1,52,..-,5 c 
K( a il tous ceux des termes qui sont de la forme 
S1,52,---,5n 


+K(Sq, 8g) . K(s,, Sy) Sigice K(s), Su) . K(sy, Sa) 
ot les indices a, 3,7,...,A,u,a forment un cycle de n +1 lettres. 


16. — S’il s’agit d’équations intégrales de la forme (1’) (p. 5), 
on admettra — pour se borner toujours aux seuls exemples vraiment 
usuels — que le noyau peut devenir infini lorsque les points M, P 
dont il dépend coincident, et seulement dans ce cas; que, de plus, 
il est alors de l’ordre de a, a étant encore un exposant constant 
et r désignant la distance MP. Cela revient dire — si, par exemple, 
l’équation est a deux variables, de sorte que M est défini par deux 
coordonnées s1, s2 et P par deux coordonnées t;, t2, — que le noyau 
est de l’ordre de : os 

[(s1 — t1)? + (sg — t2)?] 2 

Ce cas sera étudié dans une note 4 la fin du volume (note A, p. 155), 
ou nous démontrerons que les méthodes exposées au n° précédent 
réussissent pour a < 2 si, comme nous venons de le supposer, le 
nombre des variables est de deux. 


17. Etablissement de quelques formules. — Nous ne nous 
sommes pas occupés encore du cas exceptionnel ot \ est égal a 
une constante caractéristique du noyau K(s,t). Pour y arriver, nous 
commencerons par le cas de l’équation intégrale homogeéne, celle qu’on 
obtient en supposant dans (1>i), p. 56, que f(s) = 0. Il nous sera 
nécessaire auparavant d’établir quelques formules. 

Tout d’abord, on tire évidemment des développements (40), (41) 
(valables quel que soit 2) 


b b b pb 
[o(arje- | K(s,s)ds—a | [K (22) asasi 4 
Ja s Ja Ja Ja 8, 5] 


et par suite 


(50) £pe)=- f'D(2a) de 
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formule qu’on peut écrire aussi d’aprés (43) 


: b 
(50bis) Rae) = -| K(s, s, \) ds. 
dd , 


18. — La formule (50>) permet d’établir le développement suivant 
les puissances de \ du logarithme du déterminant D(,). Si, en effet, 
on utilise pour K(s,s,\) le développement obtenu par notre premiere 
méthode, — c’est-a-dire celui que l’on déduit de la formule (8) (n° 3, 
p. 40) en y changeant K en \K et k(s,t) en —AK(s,t,) d’apres le n° 11 
— il vient (|\| étant supposé plus petit que le module du premier zéro 


de D(\)) 


d b 


(50°25) a [log D(A)] = 


a 


b 
K(s,s)ds+A | Ko(s,s)ds+... 
a 
b 
+9n f Kn(s,s)ds+.... 
a 
D’oti (puisque D(0) = 1) 


Nr b 
log D(A) = -S- — | Kn(s, s) ds. 
n a 


19. THEOREME. — Si \ est une constante caractéristique de 
AK(s,t), c’est, pour quelque point (s,t), un pole de la résolvante 


K(s,t, \) = ee) 


Car pour quil en fit autrement, il faudrait que (quels que 
soient s, t) \’ soit racine de D @ ) et méme racine d’un ordre de 
multiplicité p,; au moins égal a l’ordre de multiplicité p de 4’ considéré 
comme racine de D(\). Or on tire de (50) 


qh b qh-i 8 
pe bO) =~ f ayer P (39) 
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Si on fait dans cette identité \ =’, on aura 


ie =0 pour h=1,2,...,p1 
et par suite 
qh 
apa DA) = 0 pour \=2' et h=1,...,pr. 
Donec 
(51) p2pi tl. 


II]. — RESOLUTION DE L’EQUATION HOMOGENE 


20. — Considérons maintenant l’équation intégrale homogéne 


b 
(ate) p(s) —d | K(s,t) y(t) dt = 0 


c’est-a-dire l’équation (1>8, p. 40) ot Von suppose le terme connu 
f(s) identiquement nul. Elle admet évidemment la solution ¢(s) = 0. 
Si A n’est pas une constante caractéristique, le théoreme fondamental 
du n° 11, nous apprend que cette solution est la seule. I] suffit donc 
d’examiner le cas ou \ est une constante caractéristique. 

Il est d’abord évident que s’il existe une solution non nulle de 
Véquation, il en existera une infinité d’autres qu’on obtiendrait en 
multipliant l’une d’elles par un facteur numérique arbitraire. On voit 
meme que s’il existe plusieurs solutions linéairement indépendantes 
de (1%), toute combinaison linéaire et homogéne de celles-ci sera aussi 
une solution non identiquement nulle. Reste le probleme d’existence 
qui est résolu par le théoréme suivant. 


21. THEOREME. — Si \ est une constante caractéristique du 
noyau K(s,t), ’équation intégrale homogeéne 


b 
(52) Bay i! K(s,t) y(t) dt = 0 
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admet au moins une solution ¢(s) non identiquement nulle. 

En effet, d’aprés le n° 19, la constante caractéristique \’ sera un 
pole de la résolvante K(s,t,A) d’un ordre r au moins égal a l’unité. 
Nous pourrons donc écrire cette résolvante sous la forme 


PI (s, t) 
(\’— A) 


yr(s,t) | Yr—1(s,t) 
(= A)r : (7 — Ajyr-t aie 


(53) K(s,t,A) = + yo(s,t, A) 
ot ¢,(s,t) est certainement non identiquement nulle et ott yo(s,t, A) est 
une fonction de qui est holomorphe prés de \ =’. 

Nous utiliserons maintenant la relation (46) entre le noyau et sa 
résolvante sous la forme 


b b 
— K(s,t) + K(s,t,A) = (A-’) [ K(s,7) K(7,t, A) dt + x | K(s,7) K(7,t, A) dr 


a 


=(A-2’) i K(s,7, \) K(7, t) dr + a K(s,7, A) K(z, t) dr. 


Si nous substituons ensuite l’expression de K(s,t,A) tirée de (53) 
dans (54) nous pourrons identifier les coefficients de 


OCS Ol SA es Oo 


La derniére des relations ainsi obtenues sera 
b b 

(53) erls,)=N f Kis,rertrt)dr =X [K(G,1) poloyr) ar 
a a 


Cette formule nous montre qu’en donnant at une valeur quelconque 
comprise entre a et b, la fonction ¢,(s,t) vérifie ’équation (52). Or 
yr(s,t) n’étant pas une fonction identiquement nulle en s et t, il existe 
au moins une valeur numérique @ comprise entre a et b telle que y,(s, #) 
soit une fonction de s non identiquement nulle. Cette fonction étant 
solution de (52), le théoreme est démontré. 

Ce théoreme est fondamental. I] peut en effet s’énoncer ainsi : 

Si le déterminant d’une équation de Fredholm est nul, l’équation 
homogeéene correspondante a au moins une solution non identiquement 
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nulle. Si au contraire, le déterminant est différent de zéro, la for- 
mule (47) s’applique en particulier si f(s) est identiquement nulle, 
quand on y prend encore pour ¢(s) une solution de l’équation (1) qui 
devient homogéne puisque f(s) =0. Elle montre que si D(A) 40, toute 
solution de l’équation intégrale homogeéne est identiquement nulle. En 
rapprochant ces énoncés du théoreme fondamental du n° 11, on a la 
proposition suivante : 

Pour que l’équation de Fredholm (15%) admette une solution 
(fournie d’ailleurs par la formule 47), il faut et il suffit que 
V’équation homogéne correspondante (1*") n’ait pas d’autre 
solution que zéro. 

On a ainsi pour établir l’existence ou la non existence d’une solution 
de l’équation avec second membre, un procédé souvent commode en 
ce quil évite le calcul du déterminant. 

On voit que pour s’assurer que le déterminant D(A) est différent 
de zéro [que par conséquent la méthode des numéros précédents est 
valable et fournit une solution unique de l’équation (1')|, il n’est pas 
nécessaire de calculer ce déterminant. I] suffit de constater que pour 
la méme valeur de \, l’équation homogene correspondante (1'") n’a 
pas d’autre solution que zéro. Cette remarque sera d’une application 
constante dans la suite. 

Remarque. — Nous avons obtenu plus haut dans le cas ot \ est égal 
a une constante caractéristique \’, une solution y,(s,t) de l’équation 
homogeéne (52) qui dépend d’un parametre arbitraire ¢. Il y aurait 
done une infinité de solutions, ce qui ne saurait nous étonner (n° 20). 
Mais nous montrerons plus loin (n° 24) qu’on les obtient toutes de la 
maniére indiquée précédemment (n° 20) ou d’une maniere plus précise, 
que toute solution de l’équation (52) est une combinaison linéaire et 
homogéne de solutions linéairement indépendantes en nombre au plus 
égal (n° 25) au degré de multiplicité de la racine \’ de D()). 


22. — Le théoréme que nous venons de démontrer s’applique 
encore aux cas oti le noyau devient infini, tels que nous les avons visés 
au n° 14; c’est-a-dire que dans ces cas, ou bien l’équation donnée 
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a une solution fournie par la méthode du n° 13 ou bien |’équation 
homogéne correspondante a une solution non nulle. C’est en effet ce 
qui a lieu pour l’équation (48) et nous avons vu que cette derniére est 
équivalente a la proposée. 

On étendra également ce résultat a la méthode de M. Hilbert 
(n° 15). 


23. — Pour trouver le nombre exact des solutions de l’équation 
homogéne (52) pour chaque constante caractéristique ’, on a recours aux 
mineurs du déterminant D(A). 

On appelle ainsi avec M. Fredholm les séries 


(56) D Cire bone COR ea vmera. 
We Daa eachn i, 19)..25tp . ty, to,..-,tp,T1 
Sy? ‘ S19 sixty Sata. 5 
ie nok Pf fe ie Td odin to. 

ni a 4 t1,---,€p,T1,---5Tn 


On voit comme pour les séries (40), (41), que si K(s,t) est borné, 
la série (56) est uniformément convergente, quand les s, t varient de 
facon quelconque dans (a,b) et que \ a une valeur fixe quelconque. C’est 
donc une fonction entiére de A qui est continue par rapport a (sp,tp), 
(h = 1,2,...,p) en tout point ot K(sp,t;,) est continue. 

On tire de la formule précédente (ot lon remplace 71,...,7 par 


Sptly+++5 Spin) 


b b 
pe) Dp (32°72) day dsp.» dsp 
a a S1,+++5,5p 


a Sy] Sy, Sp+1 
=f. fx een cael 4 ) 51... dsp dspaa 
a a $1,+++,Sp, Sp4+1 


—r)r f° g eee Pee 
Pn PKB | do. donap to 
n) zs A Biya end 
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D’ott en comparant avec (41) 


dP 8p 
(67) pa)= -y f.. fp y (78a) dn dp 


Développons le déterminant 


(58) K oo 


ty, t2,...,lp,T1,+++5Tn 
K(s1, t1) Sas K(s1, tp) K(s1, 71) eae K(s1, 7%) 
_ |K(sp,t1) ... K(sp,tp) K(sp,71) ... K(Sp, Tm) 
— K(7, t1) ane K(r1, tp) K(i. 71) ie K(11, 7) 
Kiet) ~ sas. Kltede) IG), «en leet) 


par rapport a sa premiére ligne et intégrons par rapport a 71,72,-.-7n- 
Nous aurons 


b b 
peat sD daca 
a a ty, t2,.--,tp,T1,-++5Tnr 


h=p 
= 31)" Kin 
h=1 
lee pl Led eecettee oe SO tated 
a - th, -ytp—1, thai,» «+5 bp, T1,+++57n 
h=n b 
+ nett [KGa ny) x 
h=1 


b 
de [x G eee rR ee. Nee A aan 
159224 lpg bps Tis ae 29 THY TRL) = 09 Te 


dt dt2...dT,-1 dtn41-. any}, 
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Mais en faisant passer la (p +h — Tjere ligne au premier rang dans 
le déterminant qui figure dans l’accolade, on voit qu’en remplacant les 
variables d’intégration Tp, Th+41,7h+2;---)Tn Pal T,Th,---,Tr—1, Sa valeur 
deviendra la méme quel que soit h, de sorte que la deuxiéme partie du 
second membre de l’égalité deviendra 


b b rb 
89,...5S8yn-1; 8 Bese ze 
nf Konr){ ff cate mE A cae ee drm hd 
a a a th, eats ytp—1, bp, T1y ++ +5 Tn—-1 
(=A)? 


n! 
de la série obtenue en portant dans (56); nous aurons 


Db, oS a) 
ty, t2,...,tp 


Multiplions maintenant l’égalité obtenue par et faisons la somme 


h=p 
SD ig tis weve ord aod one e hag ook sys Ss 
59 = aif kK 81, tp, Dp ( , : a) 
( ) 2 ) ( ) > t1,-++yth—1)th4i,.++ stp 
b 
+f K(s1,7)Dp (°°?""?"? )) dr. 
a ti, t2,...,tp 


En opérant avec une des p premieres lignes autres que la premiére 
dans le développement de (58), on aurait de méme obtenu 


Dp Gree 
t1, t2,...,ty 
h=p 


(60) = (anes K(s;,th) Dp-1 be ae sea r) 
Fz ty, t2,.--,tn—1,th4i,--- bp 


b 
815 825 e065 S$L15 Ty SELL 5 is 2 8 
+ | K(sis7)Dp ( ne eres ra) dr. 
a ty, ta,..-, 0-1, bi, biga,.--, tp 
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En opérant de méme avec les colonnes, on voit qu’on aura aussi 
D BLA aa sae py 
BN Metis t 
1, 42,+++5 lp 


S1,.++5Sh—158 nag 
(61) = (—1)*+4 K(sp, ti) Dies ( 1 9 Sh—-1) Sh+1> ’ Pa) 
ea t1,...,bj-1, tigi, .-- 5b 


b 
8450445 S715 845 $4 sb gee 
+ / K (rte) Pp ( 1; 9 84-1) 4) F741; ) ra) Ar 
a by, . 6-5 bya, T, bys 5 bp 


24. — Reprenons maintenant l’équation intégrale homogéne (52), 
A étant égal & une constante caractéristique \’. Soit donc D(X’) = 0. Nous 
avons déja observé que si l’équation 


b 
(62) ols) -¥' | K(,t) o(t) dt =0 
a q solutions ®;(s), ®o(s),...,®(s), la fonction 
(63) C1 ®1(s) + c92 ®o(s) a[ee goe Cg ®,(s) 


sera aussi une solution quelles que soicnt les constantes ¢1,C¢2,...,Cq- 

Nous allons donner une nouvelle démonstration de l’existence d’une 
solution de l’équation (62), démonstration beaucoup plus longue que celle 
du n° 21, mais qui aura l’avantage de nous donner |’expression générale 
de toutes les solutions possibles. 

Pour cela, utilisons les formules (53), (54), (55). Comme D(A) est une 
fonction entiére égale 4 1 pour \ = 0, D(A) n’est pas identiquement nulle 
et par suite les quantités 
dP 


= D(a)... 


d 


ne sont pas toutes nulles pour la valeur \ = X. Il résulte alors de la 
formule (57) que les fonctions 


S1\, $1,52 \, $1,52,-+-+,5py7 
4 D r eee DB) a eee LD r aes 
vo aby ) ; Apacs ) : a asia ) 
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ne sont pas toutes identiquement nulles par rapport aux s et aux t. 
Soit gq, Pindice de la premiére de ces fonctions qui n’est pas identiquement 
nulle. Nous allons montrer que toute solution continue de l’équation 
homogeéne (62) est de la forme (63) ot. ®;,69,...,@, sont q solutions 
linéairement indépendantes — q étant le nombre que nous venons de fixer. 
— En effet, d’aprés (60), on aura 


Dg Gare +++, Si-1,8, Si41,---5 Sq ’) 
ty, t2,...,t—-1, ti, bi41,---5tg 


b ; 
= xf K(s,7) Dg Ge eee x) dr 
a 


th, HHS »bi-1, ti, tigi, vee ytq 


$1, 52,+++,8q \ 0 
lea )# 
on peut supposer a partir de maintenant qu’on prenne pour les s et 


les t affectés d’indices, non pas des quantités variables, mais des valeurs 
numériques déterminées s’, t’ telles que l’on ait l’inégalité numérique, 


SV aavesne “4 
(66) Dg yy! a ee, 
ia epeenealy 
Dans ces conditions, on voit qu’en posant 


/ / s/ / / 
ely 0; (one ye cane) 
a / / 


be pl: / 
Se are arenes 


(65) 


Puisque 


la fonction y;(s) sera une fonction continue non identiquement nulle qui 
d’aprés (65) vérifiera l’équation homogéne (62). Il en sera de méme de 
toute fonction obtenue en multipliant y;(s) par une constante non nulle. 
Nous utiliserons les fonctions 


D S15 8-158; Si4d) »5q / 
Bh: GrpaaGal cvbecl ang aok ae 
(67)  ;(s) = ae Z (§=1,2,... 


/ 

$1) $3) 83)-+ +1 8y 1 
q tH . tht! t! 
12 Gd eG 
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comme solutions de l’équation homogéne. 
Il est facile de prouver, avec M. Plemelj, que ces qg fonctions sont 
linéairement indépendantes, c’est-a-dire que si l’on a Videntité 


(68) 1, ®i(s)+...+1, 8,(s) =0 


ou 1;,l9,...,1, sont des constantes, ces constantes sont nécessairement 
toutes nulles. En effet, d’aprés la définition de ®;(s), on a 


66.) =. Os) = 0" pour 7A kere = (12... 


En faisant maintenant s = s1,55,...,5, dans (68), on obtient successi- 
vement 


l=0, b =0, ..., 4 =0. 


Il reste A prouver que toute solution yo(s) de ’équation (62) est de la 
forme (63). Pour cela remarquons que la fonction de s 


b 
(69) Bie = i K(s,) wo(t) dt 


étant identiquement nulle, il en sera de méme de |’expression 


(70) oN [oo t) {volt —N i. K(t,7) yo(T) ar} dt 


a 


ou g(s,t) est une fonction continue quelconque. D’ot aussi, en retranchant 
les deux expressions (69), (70) 


b b 
0= ols) -¥ f K(s,t) o() at f ols,fe0lt) tt 
a : a" a 
+ eh | g(s,t) K(t, 7) yo(r) dr dt; 
a a 
ou en permutant les notations t, 7 dans la derniére intégrale 


b 
(71) o(s) —%' [H(s,t) po(t) dt = 0 
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b 
(72) fe) =k. + {atest =a | g(s;7) Kind) ar} 


Servons-nous maintenant de l’identité (61) pour 7 = 1. En y remplagant 
p par q+1,h par h+1, d par X’ et 51, t1, 5n41, thoi par s, t, 5), t}, elle 
devient : 


! , b , / 
8, 81,+++589 \7 , 5, 515+++55q y7 = 
D TN) —X K(r7,t) D IN )dr= 
aes [ a) —— 


K(s.4)D Bigs es $q K( 31, 6) Pho 8 Sh pays 8q y7 
(s, ) qd t! t -5 (s,,¢) ie / , yt t! : 
At sess ; 


er en eee 


Or, prenons en particulier pour g(s,t) la fonction 
Deut $, 845 825-++ 4 8q yy 
SN Eth excel 
9(3,1) = . 


/ / / 

$11 82)+++98q yy 
qd t! t! t! 

iaboneeraty 


L’identité précédente deviendra 


b cam! 
g(s,t) — xf g(s, 7) K(r, t) dr = —K(s,t) + SK (sh,t) p(s). 
7 h=1 
Et l’équation (72) pourra s’écrire 
h=q h=q 
His.) =Klat)— KG. — YRC t al) = Yo BN) 
h=1 h=1 


En portant dans (71), on voit que 


h=q b 
cols) =X On(s) f K(sh.t) eolt) at 


h=1 wo. 
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De sorte qu’en définitive ~o(s) s’exprime bien sous la forme (63) en 
prenant pour constantes cp, les quantités 


b 
(74) ae | KG ede. 


En résumé : Si » n’est pas une constante caractéristique, l’équation ho- 
mogéene de Fredholm (1°") n’a aucune solution continue non identiquement 
nulle. 

Si est égal ad une constante caractéristique , léquation ho- 
mogéne (62), a une infinité de solutions continues non identiquement 
nulles, et si q est l’indice de la premiére des fonctions (64) qui ne s’annule 
pas identiquement, toute solution continue est une combinaison linéaire 
de q fonctions continues linéairement indépendantes données(') par les 
formules (67). 

D’aprés ce qui précéde, il existe une infinité de systémes de fonctions 
continues tels que toute solution de l’équation intégrale homogéne (62) 
soit une combinaison linéaire et homogéne des fonctions 61,62,...,q 
appartenant a l’un de ces systémes. (Il suffit en effet de prendre pour 
61,...,9g, q combinaisons linéaires, homogénes et indépendantes des 
fonctions ®;). Nous dirons que 61,62,...,0q forment un systéme de 
fonctions fondamentales correspondant a la constante caractéristique \' et 
au noyau K(s,t). 


25. Remarque. — Le nombre q des solutions indépendantes varie 
avec \’; on l’appelle lindex de la constante caractéristique ’. D’autre 
part, \’ étant racine de D(A), on peut définir son ordre de multiplicité 
p2i. 

Or, nous savons (n° 19) que \’ est (pour quelque valeur de s, t) un 
pole de la résolvante K(s,t,) : soit r son ordre de multiplicité. On peut 
déterminer exactement q au moyen des D,;,Do,.... Nous allons voir qu’il 
suffit de connaitre r et p pour avoir une limite supérieure de gq (?). 


Cy Remarquons que la condition (66) est la seule imposée aux quantités s’, t’; 
de sorte que les fonctions ®; dépendront en général de 2q paramétres (du moins 
en apparence). 

(?) H.-B. Heywoop, Thése, p. 16. 
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En effet, appelons p,,p2,... les ordres de multiplicité de X’ pour les 


fonctions D; ic A) , D2 aes s) ,-... On a d’aprés (56) 
1, 42 


d S1,-++58k ° Siscivs SE 
dX aaa ) ‘i wisn t1,..., tk, T 7 
Pk+i < pr — 1. 


p2 <pi—1,...,pq-1 < pq-2 — 1, 1 < pg-t. 


En ajoutant les inégalités, on a 
b<pr=@=2) 08 pi -g=—: 


D’autre part, d’aprés la définition de K(s,t, A), on a 


r=p— Pi. 
D’ot: 
ce p—@ 1) 
et enfin 
(APs) q<p-rtl. 


En particulier, puisque r > 1, q < p; on a d’ailleurs par définition de gq, 
q > 1; donc & une racine simple de D(A) ne correspond qu’une fonction 
fondamentale. 


26. Résolution de l’équation intégrale singuliére avec se- 
cond membre. — Revenons a |’équation intégrale avec second 
membre (1>8), mais ol. nous supposons que \ est égal & une constante 
caractéristique \’. Non seulement le raisonnement qui nous a fourni 
expression (47) de la solution de l’équation (15'S) n’est plus valable, 
mais cette expression n’a plus de sens dans le cas actuel. Nous allons 
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montrer qu'il n’y a pas alors en général de solutions continues de 
Véquation 


b 
(75) p(s) — | K(s,t) p(t) dt = f(s) 


si f(s) est une fonction continue quelconque. Mais il est évident que 
s'il y a une solution y(s) de cette équation et si 41,62,...,@, sont les 
solutions indépendantes de |’équation homogéne correspondante (62), 
la solution la plus générale de (75) est de la forme 


p(s) + .e1 01(s) +... + eg Aq(s) 


ou c1,c2,..., sont des constantes arbitraires. 

Pour étudier les conditions d’existence d’une solution de (75), nous 
allons introduire I’équation associée a Véquation homogene (62); ce 
sera par définition l’équation 


b 
(76) w(s) — xf K(t, s) w(t) dt = 0. 
a 
A chaque solution x(s) de cette équation correspondra une condition 
de possibilité de l’équation donnée (75). 


Supposons en effet que l’équation (75) ait une solution continue ¢(s). 
On aura en multipliant par y(s) et intégrant 


[sox sjas= [ols xoyas—¥ PK (s,t) p(t) x(s) ds dt 
falas K(t,s) ialen 


Vaccolade étant nulle par définition de y(s). Il faut donc qu’on ait 


[ “FOx@aex 


ce qui n’aura pas lieu pour une fonction f(s) quelconque. 
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27. — Pour l’équation (76), le déterminant D(A) est le méme que 
pour l’équation (62); les fonctions K C ) joie Dg & ae - Sa 
Lye bp 


pour (76), s’obtiennent en permutant les s avec les t dans les fonctions 
correspondantes pour (62). I] en résulte que si \ est égal A une constante 
caractéristique 2’ de (62) correspondant 4 q solutions fondamentales, 
’ sera aussi une constante caractéristique de (76) qui correspondra a 
q solutions analogues aux ®; 


/ / / / / / 
D $1, 895+ ++) 844) 85) Sigs +++) 8g \! 
ale! t! t! t 
Anevenee ta, Stig) 


(ot l’expression D, est celle qui correspond & (75) comme jusqu’ici). 
Soient X1,-..,Xq wn systéme de fonctions fondamentales de (76) 
correspondant 4 ’. On aura d’aprés le n° 26 les q conditions de possibilité 


b 
(78) [ fera)ds=0 @=1,2,...,0). 


Supposons maintenant que la fonction f(s) satisfasse A ces q conditions. 
Les W; étant des combinaisons linéaires des y;, les conditions (78) peuvent 
s’écrire en remplacgant les x; par les Wj. 

Connaissant déja la forme (9), p. 40, de la solution unique, dans le cas 
ot. A n’est pas une constante caractéristique, il est naturel de chercher si 
Yon pourrait trouver une solution de l’équation (75) qui serait de la forme 


b 
(79) ls) = #05) -X f Gls,t) F(t) at 


ot G(s,t) est une fonction & déterminer. 
Il faudrait pour cela que V’on eut [en substituant cette expression 
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de y(s) dans (75)| 
b b 
f(s) / G(s,t) Fld) dt — / K(s,t) f(t) dt 
b b 
+22 f ii Kener nrodaare 


ou 


b b 
(80) i, f(t) {Gb t) + K(s,t) — xf K(s, 7) G(r, t) ar} dt = 0. 

a a 

Si \’ n’était pas une constante caractéristique, on satisferait 4 cette 

équation en choisissant pour G une fonction telle que l’accolade fut 
identiquement nulle, il suffirait de prendre pour X’G la fonction réciproque 
de X’K et on retomberait sur la formule (47) déja obtenue. Mais, X’ étant 
une constante caractéristique, et puisque nous admettons les conditions (78) 
réalisées, il suffit qu’on puisse choisir G de facon que l’accolade soit une 
combinaison linéaire des W(t) (a coefficients indépendants de t, sinon 
de s) pour que Videntité (80) ait lieu. Or la formule (59) peut s’écrire 
pour p=q-+1 sous une forme analogue a (73) 


I / b / / if 
S,S1,.-.-5,8 T,81,89,.-.-,8 
Dg+1 is : nN} — x | K(s,7) Dg+1 a . ; N | dr 
ees a bt Gypnccgth 


roo ! 
S15 Bbyser gS 

= K(s, t) Di (? u - ") 
toby ++ bg 


= IGE contrat tech as 2 cplusiaatee s! 
+ 5(-1)*K(s, 44) Dy (°17?” oe 


(81) 


L’EQUATION DE FREDHOLM 85 


ou en tenant compte de l’expression (77) des WU; 


Donc en prenant 


(82) G(s,t) = 


on voit que le premier membre de (80) se réduira a 


h=q 5 
DDK (sth) [FO Yalt) dt =0 
h=1 Ja 


d’apreés (81) et (82). 

En résumé, si ’ est une constante caractéristique, lV’équation avec 
second membre (75) n’a pas en général de solution. Pour qu’elle ait une 
solution continue, il faut que f(s) satisfasse aux q conditions 


b b 
(83) J teraw)as=o, oe / f(s)xXq(s) ds = 0 


ou les x; forment un systeme de solutions indépendantes de l’équation 
homogeéne associée. Et alors il y en a une infinité donnée par la formule 


By Sassen Se 
t,s saute =e 
usd 


qi ; Ch On(s) 


Dai 
p(s) = f(s) —X i fo 
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ou les 8, sont les q solutions indépendantes de l’équation donnée privée de 
second membre. 


28. Remarque. — Nous avons vu que toute constante ca- 
ractéristique \’ de l’équation (62) était constante caractéristique de 
léquation associée (76) et que le nombre des solutions linéairement 
indépendantes de ces deux équations intégrales homogenes est le 
méme. 

On peut montrer en outre que si deux solutions respectives de 
deux équations intégrales associées ne correspondent pas a la méme 
constante caractéristique, elle sont orthogonales (voir n° 7, p. 8). Car 
silona: 


6(s) =r [8.000 dt 


U(s) =” [ K(t, s) U(t) dt 


on aura 
(N — 2X") i 6(s) U(s)ds = Xr" i i 0(s) K(t, s) U(t) dt ds 
a i. [vowe.nao dt ds. 


Or la derniére intégrale devient égale a la précédente quand on 
y permute les variables d’intégration. Le premier membre étant nul, 
avec \’—” £0, les fonctions 6, Y sont bien orthogonales. 


28bis, Cas exceptionnel de M. Picard. — Nous avons vu (n° 14, 
p. 64) que sous des conditions trés larges, on peut étendre la méthode de 
Fredholm au cas oti le noyau devient infini dans le domaine d’intégration. 
On pourrait étre amené a penser qu’il est de méme facile d’étendre la 
méthode au cas ot le domaine d’intégration n’est pas borné. M. Picard 
a montré qu’il n’en est rien : non seulement la méthode de Fredholm 
suppose essentiellement que le domaine est borné, mais les propriétés 
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qu’elle permet de démontrer peuvent devenir complétement inexactes dans 
le cas contraire. Sans préciser toutes les particularités nouvelles qui peuvent 
se présenter dans ce cas exceptionnel et pour lesquelles nous renvoyons a 
la note de M. Picard (‘), nous nous contenterons des indications suivantes 
qui suffisent pour mettre en garde contre une généralisation hative. 
Considérons l’équation différentielle linéaire du second ordre 


(84) uy — pu(t) = f(t) 


ot u(t) est la fonction inconnue et p une constante donnée, et servons-nous 
de Videntité 


u; v(t) — u(t)o, = 5 [u,v(t) — u(t)u,]. 


Nous y remplacerons u(t) par une solution de l’équation (84) et w(t) 
par une solution de l’équation 


— p(t) = 90, 


ou p est une constante réelle positive. Prenons méme en particulier la 


solution 
u(t) = VE (Rs) (eg = 41) 


qui dépend d’une constante arbitraire s. L’identité deviendra 
(0 — wh etVEO-9) u(t) + ef VF) F(t) = © [af — eyault)) VFO) 


En intégrant par rapport a t de —oo a s pour ¢ = +1 et de s a +00 


pour € = —1, on aura si u(s) et u’(s) sont bornées et > 0 
(p — ») | e VE(S—)) u(t) dt +f e VES) F(t) dt = ul, — /uu(s), 


oo +00 
(p — n) [ e VEE) a(t) dt +f e VES) F(t) dt = —ul, — /uuls), 


et en ajoutant : 


+00 +00 
(85) (p—p) i, eo VEIS—t u(t) dt + ‘ e VE Ist F(t) dt = —2,/pu(s). 


—oo —oo 


(1) E. Prcarp, Comptes Rendus du 13 octobre 1910. 
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Prenons en particulier 44 = 1, on voit que toute solution u(s) de (84) 
qui est bornée ainsi que sa dérivée est solution de l’équation 


(a6), a= (=) |. ie elt u(t) dt — 5 ‘a ie elt f(¢) at, 


laquelle est du type de Fredholm (en posant A = +) mais avec des 
limites d’intégration infinies. 

Or, pour cette équation, les propriétés essentielles de l’équation de 
Fredholm cessent d’étres vraies. 

Supposons en effet p > 0, c’est-a-dire \ < alors on peut faire u = p 


dans la formule (85) de sorte qu’on a 
+co 
(87) u(s) =-——= eV Ist £(t) dt. 


Si done f(t) est une fonction bornée non nulle, le second membre est 
une fonction bornée (ainsi que sa dérivée premiére), qui est solution de (84) 
et par conséquent de (86). On a donc pour solution d’une équation de 
Fredholm, une fonction qui n’est évidemment pas en général méromorphe 
par rapport au parametre X. Il y a plus : les singularités de cette solution ne 
dépendent plus exclusivement du noyau comme dans le cas général (n° 11, 
p. 59), mais aussi du second membre. 

Prenons en effet par exemple f(t) = —cosnt, ot n est une constante 
arbitraire; on voit alors par un calcul facile que le terme indépendant de u 
dans |’équation (86) est ici : 


‘ge cos ns 
PEST ape 
et que la solution bornée de l’équation intégrale 

-+-0o 
88 = d —|s—t| t dt = COs NS 
(88) u(s)—> fo ertHu(t) at = 


sera, d’apres (87), ou plus simplement d’aprés (84) 
cos ns cos ns 


us) = = 


~ ptn2  14+n2—2n 
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14+7n? 


On voit bien que le pdle de cette fonction de A, soit A = dépend 


essentiellement de la fonction y(s) puisque le premier membre de (88) ne 
contient pas la constante n. 

Supposons enfin f(s) = 0; autrement dit, occupons-nous de 1|’équation 
intégrale homogéne 


(89) Neca / 7 It u(t) dt = 0. 


D’aprés ce qui précéde toute solution de l’équation 
iW I 
u,—pu=0 ou u,—(1—-2A)u=0 


qui est une fonction bornée (ainsi que sa dérivée ui) sera une solution 


de (89). Des lors, pour p < 0 ou A> 3 la fonction cos(V2A — 1.8) est 


1 
une solution non nulle de (89). De sorte que toute valeur de \ > 3 peut 


étre considérée (n° 21, p. 71) comme une constante caractéristique de (89). 
Ainsi, contrairement encore 4 la théorie générale (n° 10, p. 57), on a une 
équation intégrale homogéne (mais & limites infinies) dont les constantes 
caractéristiques ne sont pas isolées. 


i L’EQUATION DE FREDHOLM A NOYAU SYMETRIQUE 
ET LES DEVELOPPEMENTS EN SERIES DE FONCTIONS 
FONDAMENTALES 


29. — Etant donnée Véquation de Fredholm 
b 
(ibis) el) — af K(s,1) o(t) dt = F(6) 


on dit que son noyau AK(s,t) est symétrique, si l’on a 


K(s,t) = K(E, s). 
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On voit alors que cette équation ne differe pas de son équation 
associée (2), p. 1, et Von concoit que ce fait puisse amener des 
simplifications. Cette prévision est confirmée pleinement par les 
travaux de Hilbert et de Schmidt principalement, qui ont examiné 
en détail ce cas particulier. Ils ont été guidés par le rapport étroit 
entre la question actuelle et la théorie des formes quadratiques et de 
«l’équation en s» (1). En suivant encore la méthode qui a été exposée 
au n°8, p. 42, on voit en effet que si le noyau est symétrique, les 
équations (36>), p. 42, ot lon remplace K(s,t) par AK(s,t) peuvent 
étre écrites sous la forme : 


a) 
Lp — Az— K(H1, @2,..., Ln) = f(Sp) (p =1,2,...,n) 


Oxy 
ol. A(z1,...,%n) est une forme quadratique en 2z}1,...,2n. Et le 
déterminant des coefficients de 2,...,2, fournit quand on légale 


ee 1 E ‘ : 
a zero en y remplacgant \ par —, une équation de la forme dite 
Ss 


«équation en s» (1). M. Hilbert a procédé en suivant la méthode de 
passage a la limite, fondant en méme temps la théorie des formes 
quadratiques infinies. M. Schmidt a ensuite attaqué directement la 
méme question par la méthode des fonctions orthogonales (déf. au 
n° 7, p. 8). L’importance de leurs résultats consiste surtout en ce que 
le cas du noyau symétrique est tres fréquent dans les applications. 


30. THEOREME. — Deux solutions y(s), 0(s) d’une équation 
intégrale homogéne a noyau symétrique, qui correspondent a deux 
constantes caractéristiques différentes ’, X”, sont orthogonales (2). 


31. THEOREME. — Un noyau symétrique réel ne peut avoir de 
constantes caractéristiques imaginatres. 

Supposons en effet que le déterminant D(\) du noyau symétrique 
K(s,t) ait une racine imaginaire \’ =a +73; D(A) étant a coefficients 


(') Voir par exemple, Encyclopédie des sc. math., I, 11, p. 390. 
(7) La proposition actuelle peut étre considérée comme un cas particulier de 
celle du n° 28. 
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réels admettra aussi la racine conjuguée ” = a—i3. Nous avons 
vu que l’équation homogéne a au moins une solution continue non 
identiquement nulle 


(s) = P(s) + 1Q(s) 


quand \ est égal a une racine 2’ de D()). 
Si dans (62), p. 75, on remplace i par —i, on voit que la fonction 


A(s) = P(s) — iQ(s) 


est une solution non identiquement nulle de l’équation homogeéne (1°, 
p. 69), correspondant 4 une constante caractéristique \” = a — ip 
différente de d’. 

D’aprés le théoréme précédent, on a donc 


b b 
o- | o(s) (8) ds = f (P? + Q?) ds 


et la fonction continue P? + Q? = y(s) @(s) devrait étre identiquement 
nulle sans qu’aucun de ses facteurs ne le soit. 

La propriété actuelle généralise le fait qu’une «équation en s» a 
coefficients réels, n’a pas de racines imaginaires. 

De ce que le déterminant D(\) n’a pas de racines imaginaires, il 
ne s’ensuit pas nécessairement qu’il ait des racines réelles. On peut 
donner un exemple trés simple d’un noyau continu réel qui n’a aucune 
constante caractéristique. Il suffit(') de prendre K(s,¢) = sinscost en 
réduisant l’intervalle (a,b) a (0,27); car alors D(\) =1, tous les termes 
de son développement (41) s’annulant, sauf le premier. 

L’équation de Volterra fournit un autre exemple important. On 
voit en effet d’apres le n° 7, p. 51, que cette équation a une solution 
continue quelle que soit la fonction continue f(s) et pour toute valeur 
du paramétre A. I] n’en peut étre ainsi que si son déterminant D(A) 
n’a pas de racines réelles. En fait, tous les déterminants tels que 


(') Heywoon, Loc. cit., p. 32. Journ. de Math. (6° série), t. IV, p. 310. 
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S1,52,---,5n 
éléments principaux, ON a 


arate) (voir formule (42), p. 56) tant réduits a leurs 


D(A) <e A Sa K(s,8) ds. 


32. — On voit alors importance de cette proposition. 

THEOREME. — Tout noyau symétrique posséde au moins une 
constante caractéristique (+). 

Je suppose bien entendu que le noyau K(s,¢) n’est pas identiquement 
nul. Nous allons d’abord prouver qu’aucun des noyaux réitérés n’est 
identiquement nul. La formule (10), p. 40, montre immédiatement 
qwils sont symétriques et la formule (11) montre que si l’un est nul 
identiquement, il en est de méme des suivants. Soit dans ce cas, K2,, le 
premier des noyaux de rang pair qui est identiquement nul. On a 
d’apres (11) 


Kam(s, 8) = i: Km(s,7) Km(7, 8) dr 


et a cause de la symétrie 


b 
Kam(s.s) =f [Km(s,7)]? dr. 


Il faudrait donc que K,,(s,r) soit identiquement nul comme Ko, 
et comme m ne peut étre égal a 1, Kom ne serait pas le premier 
noyau réitéré identiquement nul et de rang pair. De plus la formule 
précédente montre que non seulement aucune des fonctions K,(s,t) 
nest identiquement nulle en s et ¢, mais encore aucune des fonctions 
Kn(s, 8). 


(') Ce théoréme a été énoncé par M. Schmidt sous la forme suivante : toute 
équation intégrale homogéne & noyau symétrique continu posséde au moins une 
solution continue non identiquement nulle. Grace a la théorie de Fredholm, 
cet énoncé équivaut a celui du texte et, comme l’a montré M. Kneser, la 
démonstration de M. Schmidt s’en trouve un peu simplifiée. 
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D’autre part, d’aprés la formule (50°), on a 


1 od 


co b oo 
(90) “DO D D(A) = y oy Kn41(s, 8) ds = » aka 


en posant : 
b 
tas f Kn(s, s) ds. 
a 


La série du troisicme membre de (90) est convergente pour \ assez 
petit (n° 12, p. 61). Je vais prouver qu’elle n’est pas constamment 
convergente, c’est-a-dire que le premier membre n’est pas une fonction 
réguliere en tout point 4 distance finie. Comme D(A) et *D(r) sont 
des fonctions entiéres, on en déduira que D(A) s’annule pour au moins 
une valeur finie de \, autrement dit que le noyau possede au moins 
une constante caractéristique. 

On a d’apres l’inégalité de Schwarz (8’), p. 8 


b pb ' 
| [ Kn—1(8,t) Kn41(s, t) wnat 


<f [Kr aeoP asd x [| knerlso) aca 


Or la symétrie des noyaux nous permet de déduire de la for- 
mule (11), p. 40 


b 
Knin(s,s) = f Kn(s,t) Kp(s, t) dt. 


a 


L’inégalité précédente devient donc 


b 


2 
b b 
[Psntos < | Kana(8,s)ds x | Kon+2(s, 8) ds 
a a 


a 


ou 
[Uan]? < Usn—2 x Uan4e. 
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D’ailleurs, aucun des noyaux réitérés n’étant identiquement nul 
Vexpression 


b 
Von = | [Kn(s, 8)]* ds 


est positive et non nulle pour toute valeur de n; on peut donc écrire 


quel que soit n 
Vante 5 _U2n 


> 0 
Weg. ” “Usa6 


et par suite : 
Uanta . Us 
Urn ~ Ug’ 
Or dans la série (90) le rapport de deux termes, d’exposants impairs 
successifs, est égal en module a 


)2| Uon+2 
Van 


Il serait donc plus grand que 1 (et par suite le terme général 
d’exposant impair de la série ne tendrait pas vers zéro), pour 


Ug 
nN — 
| I> Vg, 


Nous voyons bien que la série (90) est divergente pour de telles 
valeurs de et par conséquent que le noyau admet au moins une 
U2 
U4" 

33. — La recherche des constantes caractéristiques, autrement dit, 
la résolution de l’équation D(A) = 0 se fait par les méthodes connues. 
On voit toutefois que lapplication de l'une d’elles, la méthode de 
Graffe, présente une simplification notable dans le cas actuel. 

Le procédé de Graffe consiste(t) 4 calculer une puissance des 
racines assez grande pour que leurs rapports deviennent considérables. 


Poste: . /U 
constante caractéristique comprise entre Th et 4 
4 


(1) Voir CARVALLO, Ann. Fac. Sciences Toulouse, t. II, 1890. 
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Or, ici nous savons former directement l’équation qui a pour racines 
les valeurs de \. Il suffit (n° 13, p. 63) de remplacer le noyau K par 
le noyau réitéré K,, et d’égaler a zéro le déterminant correspondant. 
L’application de la méthode consiste donc tout simplement a écrire 
cette équation en prenant m suffisamment grand. 

On peut remarquer aussi que la démonstration précédente donne un 
moyen de calculer approximativement la constante caractéristique 1 
la plus petite en valeur absolue. En effet, son carré sera d’apreés le 
n° 13, p. 63, la constante caractéristique la plus petite du noyau 
réitéré Ko(s,t). Or, le déterminant Do() de celui-ci vérifie évidemment 
Pégalité analogue a (90) 


1 d = 
Do(u) =— > Bw” Uanse. 
Do(H) du we d “~ 
Dvailleurs, d’aprés les inégalités du n° précédent, le rapport U2n 
2n+2 


tend vers une limite déterminée pio > 0, donc la série du second membre 
a pour rayon de convergence ji9. De sorte qu’on a 


J. Schur(') a généralisé ce résultat de fagon a calculer les autres 
constantes caractéristiques de proche en proche. 


34. — Nous avons vu que chaque constante caractéristique d’un 
noyau quelconque est un pole de la fonction réciproque. Ce pole peut 
étre simple ou multiple (?). Dans le cas ot le noyau est symétrique, 


(1) Math. Annalen, t. 67, 1909, p. 327. 
(7) Prenons, avec Korn, K(s,t) = as + Gt, a et 6 étant des constantes. En 


portant dans les expressions (40), (41) de D ( ) et de D(A), on voit que les 


déterminants qui y figurent sont nuls dés qu’ils ont plus de deux colonnes. De 
sorte que K(s,t,) se réduit a une fraction rationnelle en A dont les termes 
sont de degrés respectifs 1 et 2. Si l’on prend par exemple a = 0, b= 1, on voit 
facilement que K(s,t,) a, en général, deux péles distincts qui se réduisent & un 
pole double en A pour § = —3a. 
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ce pole est toujours simple. 
En effet, servons-nous des notations du n° 21, p. 71. La méme 
méthode qui nous a fourni |’équation (55) nous donnera aussi, s2 r > 1 


b b 
a ee | Kea ide’ i Kiar eer Dar. 


Multiplions (55) et cette identité par y,_1(s,t) et y,(s,t) respec- 
tivement, intégrons par rapport a s de a a b et retranchons; on 
aura : 


b b 
pi ‘ / aCe Ce a 
a a 


b rb 
(91) -¥ f i K(s, 7) Yr—1(7, t) pr(s, t) dr ds 
b 


=) or(s,t) 


Quand on permute les lettres s et + dans la seconde intégrale 
double, on n’altére pas la valeur de cette intégrale et d’autre part, 
elle prend la forme qu’aurait la premiére intégrale double si on y 
remplagait K(s,r) par K(r,s). Par suite de la symétrie du noyau, ces 
deux intégrales sont égales; si on a soin de simplifier le crochet du 
dernier terme au moyen de (55), l’égalité (91) devient donc 


b 
| K(s,T) Yr (7, #) tn ds. 


b 
/ [¢r(s,t)]? ds = 0. 


De sorte que, ayant supposé r > 1, on obtient y,(s,t) = 0 contraire- 
ment a l’hypothese. 

Remarquons qu’il ne faudrait nullement conclure de ce qui précéde 
que le dénominateur D(A) de 
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n’a que des racines simples. Au contraire l’ordre de multiplicité p 
dune racine \’ de D(A) est d’aprés la formule (74's), p. 80, au moins 
égal au nombre q des fonctions fondamentales. En fait, M. Hilbert a 
démontré Végalité de ces deux nombres : p= q('). Ce fait fournit une 
nouvelle analogie avec les propriétés classiques de l’équation en s a 
racines multiples. 


35. Développements en séries de fonctions fondamen- 
tales. Normalisation des fonctions fondamentales. — A chaque 
constante caractéristique \’ correspond un nombre fini (p. 75) de solu- 
tions linéairement indépendantes de l’équation intégrale homogene (62). 
Nous appellerons ici ces solutions : fonctions fondamentales. Si on rem- 
place celles-ci par un systeme normalisé équivalent (p. 9), l’ensemble 
de fonctions obtenu en opérant ainsi successivement sur toutes les 
constantes caractéristiques formera luicméme un systeme normalisé, 
puisque deux fonctions fondamentales correspondant a deux constantes 
caractéristiques différentes sont déja orthogonales (p. 90). 

Si on range maintenant les constantes caractéristiques suivant 
Vordre de grandeur de leurs valeurs absolues, chacune étant répétée 
autant de fois qu’elle a de fonctions fondamentales, on voit qu’on 
obtiendra une suite énumérable de constantes caractéristiques réelles 


(92) AL, A2, A3; aie 


suite a laquelle correspond un systéme principal, c’est-a-dire un 
ensemble normalisé réel 


(S) Pils), pals), y3(s), --- 


de fonctions fondamentales, c’est-a-dire telles que l’on ait 


b 
(93) Yn(s) = rn ff K(s, t) Gn (t) dt. 


a 


(*) Voir par exemple, Goursat, Ann. Toulouse 2 (10), p. 49. 


L’EQUATION DE FREDHOLM 98 


D’apres ce qui précéde, toute fonction fondamentale (s) corres- 
pondant au noyau K(s,t) est une combinaison linéaire d’un nombre 
fini de fonctions de la suite (S), 4 savoir de celles qui correspondent a 
la meme constante caractéristique que (s). 

Mais d’autres fonctions que les solutions fondamentales sont aussi 
exprimables linéairement en fonction des »y,(s) — ne correspondant 
pas, bien entendu, a la méme valeur de 4 —. Nous nous proposons 
maintenant d’étudier ce mode de représentation dont nous savons déja 
(page 11) déterminer les coefficients. 


36. Développement du noyau. — Considéré comme fonction 
de s, le noyau K(s,t) peut-il étre considéré comme une combinaison 
linéaire des yp(s) 


(94) K(s,t) = c1 y1(s) +... +¢n Yn(s) +... ? 


D’apres la méthode de la page 11, on voit qu’on aurait 


b 
on = f K(s, t) gn(s) ds 


ou en tenant compte de la symétrie de K(s,t) et de la formule de 
définition (93) des fonctions yp, 


en = SO. 
Ainsi, on aurait 
_ grils) vit) , v2(s) yet) |, Ynls) ent) , 
(95) (aye) = SREP Pa a OO ee 


Réciproquement si la suite 


(S) yi(s), vals), .-. 


désigne un systeme principal de solutions fondamentales correspondant 
aux constantes caractéristiques \1,A2,... du noyau symétrique K(s,t), et 
si la série 


i (96) vri(s)erlt) , vals)e2) |, Gls) nit) 


ay T rs Tee eT 
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est limitée ou uniformément convergente, sa somme est égale au 
noyau K(s,t). 
Soit H(s,t), la somme de la série, et posons 


Q(s,t) = K(s,t) — H(s,t). 


La fonction Q(s,t) est symétrique; si elle n’était pas identiquement 
nulle, en la considérant comme noyau d’une équation, on pourrait 
trouver (page 79), une constante c et une fonction continue 7(t) non 
identiquement nulle telle que l’on ait 


b 
= of Q(t, 7) w(7) dr 


b 
i K(t,7) on(t) al eae 


B T 
fF ent onte sa PAO) (7) ar, 
a P 


Si l’on tient compte de la formule de définition des y, et de ce 
que tout systéme principal est normé, on voit que les deux termes du 
second membre sont égaux, comme se réduisant a 


Cc 


b 
Xf only vr) ar 
D’ot enfin 
b 
(97) | w(t) gn (t) dt = 0. 


Or, ona 


nae faervenar=ef Ker u@)dr-e° x a 
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donc d’aprés (97) 
b 
v(t) = | K(t,7) (7) dr. 


Des lors y(t) serait une solution fondamentale relative au noyau 
K(t,7) et serait par suite une combinaison linéaire d’un nombre fini 
des Yn 

w(t) = a1 Yn, (t) + a2 Yng(t) +... + ag Yn, (t)- 


En multipliant par yn,(¢) et intégrant, on aurait donc d’apreés (97) 


b 
aa W(t) Yn, (t) dt = an, (p =1,2,...,q) 


d’ou : 
w(t) =0 
ce qui est impossible. 
37. — Il résulte en particulier de ce qui précede que le noyau 


pourra toujours étre mis sous la forme (95) lorsque le nombre des 
constantes caractéristiques (et par suite aussi le nombre, en général 
supérieur, des y,) sera fini. Nous avons déja vu la réciproque de cette 
proposition (p. 46). De sorte que la condition nécessaire et suffisante 
pour qu’un noyau symétrique K(s,t) puisse étre écrit sous la forme 


K(s,t) = S1(s) S1(¢) +... + Sn(s) Sn(t) 


oun est un entier fini, est que le nombre des constantes caractéristiques 
soit fini. I est d’ailleurs bien évident que les noyaux de cette forme 
sont exceptionnels. Nous complétons ainsi le théoreme du n° 32; 
non seulement un noyau symétrique posséde au moins une constante 
caractéristique, mais en général il en posséde une infinité. 


38. Développement des noyaux réitérés. — Soit (s) une 
solution fondamentale relative a la constante caractéristique \’, on a 


B(s) = [0.090 di 
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d’otl 


b b b 
Hs) =X f K(s,s1) ®(s1) ds = a? [ [ Ss K(s.has ®(t) dt 


ou d’aprés (11) 
b 
®(s) = X? i. Ko(s, t) ®(t) dt 


et en général 


b 
(98) &(s) =? i. K,,(s, t) ®(¢) dt. 


Ainsi toute fonction fondamentale d’un noyau est fonction fonda- 
mentale de chacun des noyaux réitérés correspondants. Il résulte de 
plus du théoreme précédent sur le développement de K(s,t) et de la 
formule (98) que si la série 


gils)pi(t) vols) p2(t) _ , Yn(s) Yn(t) , 


(99) ll Tee ee T Tie sud 
ne pid 


est uniformément convergente, elle représente K,(s,t). Je dis qwil y a 
convergence absolue et uniforme au moins pour p > 3(!). En effet, la 
question ne se pose que si le nombre des constantes caractéristiques 
distinctes est infini; je remarque qu/alors leurs valeurs absolues 
croissent au dela de toute limite (car ce sont les zéros d’une fonction 
entiére D(A)). On pourra donc trouver un nombre gq tel que pour n > q, 
on ait |A,| >1. Alors pour n>q et p>3 


h=n+m 


S pn(s) ee (t) eS ealsh ent) (t) 
h=n = 
sa pes yee R(t) 
> SS 
© 2)Anl bg 


(1) Ceci est encore vrai pour p = 2, mais la démonstration est différente. 
Voir par exemple Kowalewski, Einfiihrung in die Determinantentheorie, Teubner, 


1900, p. 533. 
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a? + b2 


d’apres linégalité générale |ab| < et Vhypothése que la suite 
des |A;,| ne va jamais en décroissant (p. 96). Or d’aprés l’inégalité de 


Bessel (11), p. 12, ona 


h=oo 


P 


h=1 


b ‘ b 
[socnontoat < | [K(s, t)]? ae; 


en transformant le premier membre d’apres (93), on a 


(101) 


ce qui donne pour (100) 


h=n+m 


Ly 


Ls 


Ph (s) Pni(t) 
a 


1 b 2 
< ma. [K(s, t)]~ dt 


lorsque p > 3, n> q. Quand n croit indéfiniment, le second membre tend 


vers zéro avec ae D’apreés le théoréme de Cauchy sur la convergence 


rT 
des séries(!), cela suffit pour démontrer ce que nous avions en vue. 
Ainsi, pour p >3 (voir la note précédente), on a toujours 


gil(s) yilt) _ vals) yalt) Yn(s) Yn(t) 


102 K,(s,t) = ae 
( ) p(s ) Bu iy ae a 2 


ty sie 


ou le second membre est nécessairement une série absolument et 
uniformément convergente. 


h=oo 
@ Si une série S- un(s,t) est telle que, A tout nombre ¢ > 0, on peut faire 
RA1 h=n+m 


correspondre un entier n pour lequel x |un(s,t)| <e€ quel que soit m, la 


h=n 
série est absolument convergente et si l’inégalité a lieu quand s, t varient dans 
un domaine D, pour une valeur de n indépendante de s et de t, la convergence 
est uniforme. 
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39. — A cette occasion, on peut remarquer que la série 
1 1 1 
(103) Pat..tat... 
Mp oe rN, 


est absolument convergente pour toute valeur de p > 2. Il suffit 
évidemment de le prouver pour p= 2. Or d’aprés (101) on a 


h=m 2 b 
PA) PKs #)]2 at 
Ye < f 166.0) 


quel que soit Ventier m; d’ot, en intégrant encore une fois, en tenant 
compte du fait que les fonctions y, ont été normalisées (n° 8, p. 9) 


ham 4 b pb ' 
y= ze | [K(s, t)]” ds det. 
et Oh wee 


Ceci ayant lieu quel que soit m, prouve que la série (103) est 
convergente et meme, de plus, que pour p= 2 sa valeur est inférieure 


ou égale a Vintégrale 
b pb 
/ / [K(s, t)]? ds dt. 
a a 


Ces deux propriétés ont été étendues par J. Schur(!) au cas 
dissymétrique. Comme dans ce cas les racines peuvent étre imaginaires, 
il y a lieu de remarquer que la série (103) reste absolument convergente, 
mais que dans le cas de p= 2, c’est pour la série des modules qu’on a 


1 1 
wea K t)| 24 dt. 
Dae * Paez t Oh, pt: roi [i ee 


En outre, il est entendu que dans la série (103), chaque constante 
caractéristique \’ doit étre répétée autant de fois que l’indique son 
ordre de multiplicité comme racine de D(A). Dans le cas symétrique, 
notre démonstration remplace ce nombre par le nombre des fonctions 


(1) Math. Annalen, t. LXVI, 1909, p. 501. 
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fondamentales correspondant 4 4’. Nous avons fait remarquer (note 1, 
p. 96) que dans le méme cas symétrique, ces deux nombres sont égaux. 


40. Développements généraux. — Démontrons d’abord un 
lemme : la condition nécessaire et suffisante pour qu’une fonction 
continue I(s) vérifie Videntité 


b 

(104) | K(s,t) I(t) dt = 0 
a 

est que l’on ait quel que soit h 


b 
(105) / pn(s) l(s) ds = 0. 


En effet, supposons vraie (104), on a d’apreés (93) 


) iOIOn= 5s | 


=n i. en(t) 


b 
[ K(s,t) yp (t) “| I(s) ds 


[0.0 I(s) a dt =0 


quel que soit h. 

Inversement, partons des relations (105). La série (102) étant 
uniformément convergente pour p = 4, et égale a Ka(s,t), on peut 
Vintégrer terme a terme apres avoir multiplié par I(s) l(t); on obtient 
ainsi 


b pb jis Ga b 
J [ Kss.9u(s)1) asae = 2. a en(s)t(s) as [ pnt) l(t) dt = 0. 


En remplagant K4 par son expression, d’apres (11), on a 


of f 
fh [senna x 
-[ [satan] ae 


i Ko(s,7) Ko(r,t) in i(s) U(t) ds dt 


i Ka(z, €) l(t) “| dr 
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Pour que lintégrale de ce carré soit nulle, il faut que la quantité 
élevée au carré soit nulle, ou encore en la multipliant par I(7) et 
intégrant : 


A [ve T)I(s) U(r) dsdr =0. 


En recommengant le méme raisonnement sur Kp et K, comme nous 
venons de le faire sur Ky et Kz, on obtient bien enfin (104). 

Maintenant, nous nous proposons de développer en série de 
fonctions fondamentales, non pas une fonction absolument arbitraire, 
mais toute fonction g(s) qu’on peut représenter par la formule 


b 
(106) g(s) = f K(s,t)n(t) at 


ot p(t) est une fonction continue quelconque convenablement choisie. 
Si ’on admet pour un instant la possibilité du développement 
de g(s), on voit immédiatement (page 11) qu’il serait de la forme 


b 
Y vals) | g(t) p(t) de. 


Or cette expression peut se remplacer terme a terme par 


[ea arx f 
> [Ken ea cay a [un 
b 


b 
=) i. K(s,7) pp (r) dr x | piel ae 


b 
| K(é-m)otn) pn(t) dt 


b 
i Kin) nO) a i 


et sous cette forme le lemme de la page 12 nous apprend qu’elle 
représente une série uniformément convergente. 
Nous allons donc considérer a priori la fonction 


h=oo b 
(5) =9(3)- Yo ents) f oenttat 
h=1 Ja 
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et montrer qu’elle est identiquement nulle. 
On a quel que soit n 


b 
| pn(s) I(s) ds 
=f als)ents)ds— >  ga(s)en(s)ds x f g(t) vale) at =o 
a h=1 a a 


puisque le systeéme des y,, est normé. Donc d’aprés le lemme 


b 
(108) o= | Kesnide 
D’autre part 


b 


b b h=co ph 
I(s) "ds = I(s) g(s) ds — U(s) pp(s) ds x (t) yp (t) dt. 
[ [U(s)| : 9 d [ Ph i: g(t) eh 


Le deuxiéme terme du second membre est nul d’aprés (107), le 
premier est égal d’apres (106) a 


b 
l 
d’apres (108). Done I(s) est bien identiquement nulle. Ainsi toute 
fonction g(s) de la forme 


[0.0 I(s) ts p(t) dt = 0 


b 
(106) g(s) = i. K(s, t) p(t) dt 


peut étre représentée par une série absolument et uniformément 
convergente de fonctions fondamentales de la forme 


h=oo 


b 
(109) ae) = 5 | faorentoat} ents 


h=1 
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ou 
b 
= i p(t) pn(t) dt 
(110) g(s) = » a yn(s). 
h=1 


En remplacant g(s) par son expression (106), en multipliant par 
une fonction continue arbitraire q(s) et intégrant, on obtient une 
formule remarquable que M. Hilbert a obtenue par |’étude des formes 
quadratiques infinies 


(111) [ [x (s, t) p(t) q(s) ds dt = DE L fn t) yp (t) ya x fats) enls)as 


41. — Appelons avec M. Hilbert, noyau fermé, un noyau K(s,t) tel 
quw il n’existe aucune fonction I(s) vérifiant l’identité : 


[x0 I(t) dt =0 


Pour simplifier les résultats, nous n’imposerons pas a la fonction I(s) 
la condition d’étre continue mais seulement d’étre de carré intégrable. 

D’apres le lemme du n° 40, on voit que si un noyau symétrique 
n’est pas fermé, il existe au moins une fonction I(s) qui est orthogonale 
a toutes les fonctions fondamentales et réciproquement. C’est ce qu’on 
exprime en disant que le systeme orthogonal des fonctions fondamen- 
tales n’est pas complet. Autrement dit, un noyau symétrique fermé 
est un noyau dont le systeme principal des fonctions fondamentales 
forme un systéme orthogonal complet. 

Il en résulte en particulier qu’un noyau fermé posséde toujours une 
infinité de constantes caractéristiques distinctes. Sans quoi il n’aurait 
qu’un nombre fini de fonctions fondamentales qui devrait former un 
systeme complet. Or on peut toujours évidemment former une fonction 
qui soit orthogonale 4 un nombre fini de fonctions données. 
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En outre, un noyau fermé est tel qu’on peut, autant que l’on 
veut, approcher de toute fonction g(s) de carré intégrable au moyen 
d’une combinaison linéaire convenable d’un nombre fini de fonctions 
fondamentales ('). L’approximation doit étre entendue ici au sens de 
la théorie des erreurs. C’est-a-dire qu’étant donnés la fonction g(s) et 
le nombre positif «, on peut choisir des nombres cy, c2,...,cn, tels que 
Von ait : 


b 
J [9(s) - ex eu) - e2.921s)...— en Pn(s)]?ds < 


42. Développement de la solution de l’équation de Fred- 
holm 4 noyau symétrique. — Soit ]’équation 


b 
(1bis) ls) = (8) + J K(,4) (tat 


ou K(s,t) est une fonction symétrique. S’il y a une solution continue, 
on peut l’écrire 

9(s) = f(s) + g(s) 
en posant 


b 
ata af K(Gtoede. 


D’apres le paragraphe précédent, on peut développer g(s) sous la 
forme d’une série uniformément convergente de fonctions fondamen- 
tales 


h=oo 
9(s) = >> cn vals) 
h=l 


h=oo 


o(s) = f(s) + > en vals). 


h=1 


(2) is. FISCHER, Comptes Rendus, t. 144, p. 1148. 
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Il ne reste plus qu’a déterminer les c,. Pour cela substituons dans 
Véquation intégrale (158), ’expression trouvée 


h=oco h=o00 
+ > cavals) = f(s) + i. K(s,t) s(ae+ SO ps i K(s,t) vat) a| 
h=1 h=1 io 


ou, d’aprés (93) 


dN b 
(112) a. Ch c — 2 pn(s) = af K(s,t) f(t) dt. 


Or, d’aprés le paragraphe précédent, on a encore 


b 
; nacor | f(t) en(t)at 
[ kootoa= > | pals). 
@ h=1 


De sorte que l’égalité (112) devient 


(113) x {s (a=), (foes 08} 0) 20 


Les fonctions y, formant un systeme normé, nous avons vu (p. 9) 
qu’elles étaient indépendantes, c’est-a-dire que l’accolade précédente 
est nulle. 


1° Si done A n’est pas une constante caractéristique, on a 


® F(t) or (t) 
=f SERS at 


(1bis), elle s’obtient sous la forme 


h=o0 
i (114) v(s) = f(s)-A D> | 
p=1 


b 
py < i f(t) en(t) a pp(s) 


De sorte que, s’il y a une solution continue de |’équation intégrale 
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et le calcul précédent luicméme nous montre qu’une telle expres- 
sion vérifie bien l’équation. Dés lors, si \ n’est pas une constante 
caractéristique, l’équation intégrale (1°") a noyau symétrique a une 
solution continue unique et donnée par la formule (114). On voit de 
plus que si, pour une valeur déterminée de 4, la série 


vals) en(t) 


— Ny, A 


est uniformément convergente en ¢, on pourra écrire la solution 
correspondante sous la forme déja obtenue 


j 
(47) (8) = f(s) +A J K(a,t,d) Fb ae 
mais ici avec 
Yn(s) prt) 
K(s,t,A) = 2 ed 


L’expression (114) de y(s) présente sur celle de Fredholm l’avantage 
de mettre en évidence le caractere méromorphe de la solution par 
rapport a \, en précisant la partie principale relative a chaque pole \. 

2° Si \ est égal a une constante caractéristique ’, elle correspondra 
a un nombre fini de fonctions fondamentales 


Yn+1(s), sey Yn+q(s), 
et on aura 
NS Anti =---=Ante: 


Alors, on voit que Videntité (113) n’aura pas lieu en général, car les 
coefficients de gnii,---;Pn+q Se réduisent a des quantités indépendantes 
de Ch 


b 
(115) -f f(t) yp (t) dt (h=n+1,...,n+q) 
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en général différentes de zéro. 
Il ne peut donc y avoir de solutions continues que si ces gq quantités 
sont nulles (voir p. 81) et alors on peut prendre 


Cnt1s +++> Cntq 


arbitrairement, les autres c, étant déterminés comme précédemment, 
de sorte qu’on aura pour solution 


(116) o(s) = f(s) +en41 Pn4i(s) +... + entqPntq(s) 


h=n h=oo \ b 
4 (Ss > face f oentyath 


h=1 h=n+qt+1 


Ainsi lorsque \ est égal a une constante caractéristique \’, il n’y 
a pas, en général, de solution continue de l’équation intégrale (1's). 
Il ne peut y en avoir que si (A étant égal & Andi, Ant2,---,Antg); les 
q conditions 


b 
[ fen ar=0 (h=n+1,...,n+q) 


sont satisfaites. Nous retrouvons ici les résultats de la page 84. Mais 
de plus, nous voyons que dans le cas symétrique, sil y a une solution 
continue de (1"%), elle est donnée par la formule (116) ou les c sont 
q parametres arbitraires. 


43. Noyau dissymétrique. — M. Schmidt a étendu au cas dis- 
symétrique la notion de fonctions fondamentales. Etant donné le noyau dis- 
symétrique K(s,t), il appelle couple de fonctions fondamentales conjuguées 
un couple de deux fonctions y(s), v(s) non identiquement nulles et 
satisfaisant aux deux équations 


b 
(117) Bax / K(s,t) w(t) dt 


b 
(118) w(s) = | K(t, s) y(t) dt. 
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On voit qu’en éliminant w, puis y entre les deux équations (117), (118), 
on aura : 


b 
eyes? / K(s,t) v(t) dt 


b 
(s) =? | K(s,t) a(t) dt, 


a 


en posant 


b 

Kis.t) = | K(s,7) K(t,7) dr 
b 

Kiset) = ih K(t, s) K(z, t) dr. 


Donc ¢ et w sont les solutions de deux équations intégrales homogénes 
dont les noyaux K et K sont évidemment chacun symétrique. Il en 
résulte que les équations (117), (118) ne peuvent étre vérifiées que pour 
certaines valeurs de » (Aq, A2,..-,An;---) telles que A? soit une constante 
caractéristique de K et de K. Ces valeurs de \ que nous pourrons 
appeler constantes fondamentales du noyau K sont en général distinctes 
des constantes caractéristiques; et les fonctions fondamentales ne sont pas 
en général des solutions de l’équation de Fredholm homoggéne (1°), p. 69. 

D’aprés ce qui précéde, A? est réel. M. Schmidt démontre méme que 
\? est positif de sorte que A est aussi réel. Il prouve ensuite qu’une 


fonction de la forme j 
g(s) = / K(s, t) h(t) dt 
a 


se développe en série uniformément et absolument convergente de fonc- 
tions y(s) comme au n° 40. Et de méme une fonction de la forme 


b 
n(s)= | K(t, s) hi (t) dt 


se développe en série de fonctions w(s). 
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Il montre aussi que si la série 


Pn($) Yn(t) 
a ane 


est uniformément convergente, elle représente K(s, t). 

M. Schmidt est parvenu aussi 4 établir plusieurs des résultats que 
fournit la méthode de Fredholm dans le cas dissymétrique en ramenant ce 
cas au cas symétrique. 


Il montre que les solutions y(s) de l’équation dissymétrique 


b 
wis) =x | K(s,t) p(t) dt 


sont aussi des solutions de l’équation 


b 
p(s) = / Q(s,t, X’) ple) dt 


et réciproquement, en désignant par Q(s,t, A), la fonction symétrique : 
b 
G6 NS Ks 4 KG) i K(r, s) K(r,t) dr. 
a 


Mais, il faut remarquer la difficulté introduite par le fait que le noyau 
symétrique AQ(s,t, A) contient le paramétre \ autrement qu’en facteur. 


44. — On peut aussi appliquer directement au cas dissymétrique la 
méthode des fonctions orthogonales. Elle permet de retrouver le résultat 
fondamental de Fredholm : l’équation intégrale (155) a une solution continue 
unique, sauf pour des valeurs exceptionnelles de A; et c’est seulement pour 
ces valeurs exceptionnelles de 4 que la méme équation (15), privée de 
second membre, admet une solution continue non identiquement nulle. 

L’intérét de cette méthode consiste en ce qu’elle s’applique aussi 
facilement en remplacant partout la condition de continuité pour les 
fonctions considérées par la condition infiniment plus générale que leurs 
carrés soient sommables (intégrables au sens de M. Lebesgue (*)). Il faudra, 


(') LeBesGur, Legons sur Vintégration. Paris, 1905. 
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par contre, considérer comme identiques deux fonctions f1(x), fo(x), telles 
que 


b 
/ Lfa(a) — fale) |? de = 0 


(égalité qui, en fait, n’empéche pas deux fonctions discontinues fi, fo, 
d’étre différentes en une infinité de points). 

On voit alors directement sur l’équation (1°!) et au moyen de l’inégalité 
de Schwarz (8), p. 8 qu’en tout point so ot f(s) et K(s,t) sont a la fois 
continues, la solution y(s) est aussi continue (*) 


45. Noyau de la forme K(s,t) p(t). — M. Erhard Schmidt a 
montré que le cas du noyau K(s,t) p(t), ol. K(s,t) est symétrique et 
p(t) est une fonction qui est toujours positive ou nulle, se réduit au 
cas symétrique. Considérons |’équation 


b 
ols) —A i K(s, t) p(t) v(t) at = f(s). 


Multiplions les deux membres de cette équation par la fonction 
+,/p(s) qui est réelle et uniforme. On a 


b 
9(s) Vp(s) — af K(s,t) Vp(s) Vt) « o(t) Vp) dt = f(s) v/s). 


Si nous prenons pour fonction inconnue la fonction ¢(s) \/p(s), cette 
équation devient une équation de noyau symétrique K(s, t) \/p(s) \/p(t). 

Cette généralisation a une certaine importance dans les applica- 
tions, comme nous le verrons plus tard. 


(') F. Riesz, Comptes-Rendus du 8 avril 1907. 


CHAPITRE III 


RESOLUTION DES PROBLEMES POSES AU PREMIER 
CHAPITRE 


I. — INTRODUCTION 


1. — Nous commengons par un résumé des résultats du_se- 
cond Chapitre dont nous nous servirons apres les avoir transformés 
conformément au n° 3, p. 4. 

Les équations associées de Fredholm 


(1) (Mt) —2 f o(P) K(M.P) dop = FM), 
(2) w(M) — d [ v@) K(P, M) dop = g(M), 


ot. lintégrale s’étend a une surface S fermée dont dop est l’élément 
d’aire au point P, et ot M est un point fixe de la surface, ont chacune 
une solution unique en général. (Dans ces équations, lintégration 
peut également s’étendre 4 un domaine D limité par une surface S, 
dop étant un élément de volume). 

Si nous faisons f(M) = 0, la solution correspondante y(M) sera nulle 
identiquement en général. 

Il y a cependant un nombre fini ou infini de constantes ca- 
ractéristiques 


(x) Adie Dag eek Ange oes 


(qu’on peut ranger par ordre de modules non décroissants), telles que 
pour A =Aj4,A2,... les équations homogenes 


(2!) e(M) —» [ o(P) K(M,P) dap = 0 


(3') )-a f ve) K(P, M) dop = 0 
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ont chacune une solution, ou un méme nombre fini de solutions 
linéairement indépendantes différentes de zéro (n° 21, p. 70). 

Nous supposerons que dans la suite (©) chaque constante ca- 
ractéristique est répétée autant de fois qu’il lui correspond de solu- 
tions indépendantes. A la suite (©) correspondront donc deux suites 
de solutions associées de (2’), (3’) 


Pl, P2, PB, +++, Yns --- 
1, W2, U3, ---, Un, ---- 


On peut multiplier chacune de ces fonctions par une constante 
arbitraire sans qu’elles cessent d’étre des solutions de (2’), (3’). 

Si les constantes caractéristiques sont toutes différentes, ces fonc- 
tions sont biorthogonales (n° 28, p. 85), c’est-a-dire qu’on a 


(3) [erm P) an (P)dop =0, pour mm’. 
s 
On choisit la constante arbitraire de fagon que 
(4) [ emlP) ¥m(P) dew =1. 


Dans le cas ot il y a plusieurs (soient q) paires de solutions 
indépendantes correspondant 4 une constante caractéristique An, on 
8 hn = DW =H ow See = Aegis Cb les equations (3), mont: plus én 
nécessairement. Mais on peut s’arranger de sorte que les équations 
précédentes soient toujours valables. I] suffit d’effectuer sur les q paires 
de solutions qui correspondent 4 \, une transformation linéaire dont 
les coefficients se détermineront par une méthode analogue a celle du 
n° 8, p. 9. 

Les solutions des équations non homogeénes (1), (2) pour les valeurs 
ordinaires de \ sont déterminées par une fonction K(M,P,.) appelée 
résolvante du noyau K(M,P). Cette fonction est méromorphe en \ et 
a pour poles les constantes caractéristiques (n° 10, p. 57). Lorsque 
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nest pas une constante caractéristique, les équations (1), (2) ont 
chacune une seule solution (n° 11, p. 61), donnée par les formules 


y(M) = f(M) +2 [ K(M, P, ) f(P) dop 


(M) = f(M) +2 | K(P,M, d) f(P) dop. 


Quand nous sommes dans le cas singulier, les equations (1) et (2) 
n’ont pas de solutions en général. Pour que |’équation (1), par exemple, 
ait une solution, il faut et il suffit qu’on ait (n° 26, p. 80) 


(5) i. f(P) ¥(P) dop = 0, 


pour chacune des solutions indépendantes w de (3’) qui correspondent 
a la méme valeur de 4. Dans ce cas, la solution ne sera pas unique : 
elle renfermera g constantes arbitraires (n° 26, p. 80). 

Réciproquement, si l’équation homogene a une solution nous 
sommes dans le cas singulier. 


2. — Dans le cas symétrique ot K(M,P) = K(P,M), les deux 
équations associées coincident. Leurs solutions yp, w, coincident donc 
aussi. Elles ne forment plus qu’une suite de fonctions fondamentales 
qui sont orthogonales (n° 30, p. 90) 


[em(P) %mu(P)dop =0, pour m Am, 
5 


La solution de (1) s’exprime sous la forme 


(6) y(M) = f(M) + Ai g1(M) + Az yo(M) +... 
ou 
(7) An = sy [ FP) enlP) dor. 


lorsque \ n’est pas une constante caractéristique (n° 42, p. 108). Dans 
le cas singulier, il n’y a pas en général de solution. Il n’y en a 
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que si f(M) est orthogonale a toutes les fonctions fondamentales qui 
correspondent a la valeur de 4. L’expression de la solution devient un 
peu plus compliquée (n° 42, p. 107). 

Si la série (6) est infinie, elle est uniformément et absolument 
convergente. Une fonction quelconque g(M) qui a la forme 


(8) g(M) = d. K(M, P) p(P) dop, 


ou p(M) est une fonction quelconque, se développe suivant une série 
absolument et uniformément convergente, 


(9) g(M) = Ai gi(M) + Ao yo(M) +... 
ou 
(10) oe i yn(P) g(P) dop. 


Quand le noyau K(M,P) est fermé (n° 41, p. 107), la fonction p(P) 
est unique pour chaque fonction g de carré intégrable. 

Le nombre des constantes caractéristiques est, dans ce cas, infini 
et le systeme des fonctions fondamentales est complet (méme n°). 


II. — LA SOLUTION DES PROBLEMES DE POTENTIEL 


3. — Nous avons vu (n°9, p. 26) que les problémes de Dirichlet 
et de Neumann se ramenent respectivement aux deux équations 


(11) om p(M) +2 [ p(P) 5? as = prc, 
(12) 2m (Ml) +A f o(P) Sas = fC), 


ou nous avons 


pour \ = +1 cas intérieur Dirichlet, extérieur Neumann, 


» A=-l » extérieur » _, intérieur » 
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Ce sont deux équations de Fredholm associées dont les noyaux 
sont respectivement 


cos 


Qnr2 


cos y 
Qrr2’ 


K(M, P) = — K(P, M) = — 


Le noyau K n’est pas continu partout. Mais s’il devient infini 
quand M et P coincident, c’est dans les conditions examinées au n° 16, 
p. 67 et aussi note A, p. 155. Pour ne pas introduire de complications, 
dans les démonstrations, nous supposerons, d’ailleurs que la surface S 
est analytique et qu’elle a en chaque point un plan tangent unique. 
En outre, nous supposerons que S$ est fermée et d’un seul tenant. 

La solution du probleme physique s’obtient en prenant le poten- 
tiel V d’une couche [double pour (11), simple pour (12)] étendue sur S 
et dont la densité p est solution de (11) ou (12). 


4. Il nous sera utile pour la suite d’écrire les formules de 
Green (') pour le domaine D intérieur 4 S et le domaine D’ extérieur 


aS. 
OV\?  (av\? | (avy? 
Ox ' \ ay ' \ dz 
=- [vi~as- fff V AV dz dy dz 
s On; D 
av\? (avy? 
dy) ° \ az 
= [ves V aS at V AV dz dy dz. 
© One D’ 


Quand on prend pour V une fonction harmonique (en particulier, 
un potentiel de simple ou double couche), AV est nul dans D et 


dx dy dz 


dx dy dz 


——~ 
Q|_& 
al< 

NY 

i) 


(*) Voir par exemple : APPELL, Traité de Mécanique, tome III, p. 7. 
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dans D’ de sorte qu’il ne reste plus dans chaque second membre que 


la premiére intégrale. On voit alors que si V; ou wv est nul en chaque 


On; 
point de S, la premiére intégrale est aussi nulle et par suite, V est 


a ee i : OV 
constant a lintérieur de S. De méme si Ve ou est nul en chaque 


One 
point de S, V est constant a l’extéricur de S. 

5. — Montrons maintenant que \ = —1 est une constante ca- 
ractéristique de (11) et (12). Il suffit d’aprés le n°1 de le prouver 
pour (11). Conformément au n° 21, p. 71, nous montrerons dans ce 
but que l’équation homogene, 


cos 


(15) 2n (at) — f o(P) “7 dS =0 


a une solution différente de zéro. 
Or si da est l’angle solide déterminé par le point fixe M et |’élément 
de surface dS, nous aurons 


cos 


da = ds. 


r2 


Il faut donc satisfaire a l’équation 
(16) 2m p(M) — [ p(P) da = 0. 
8 


Or on voit de suite qu’on a la solution non nulle 


(17) p(M) = constante. 
Puisque \ = —1 est une constante caractéristique, l’équation 
(18) 2n p(M) — ‘| o(P) ay dS =0 
S r 


aura aussi au moins une solution non nulle p;(M). 

Interprétons ces deux solutions du point de vue physique. En 
appelant V le potentiel correspondant a la couche de densité p 
répandue sur S, on sait que : 
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1° S’il s’agit d’une double couche, le premier membre de (15) est 
égal a —V,-, comme on le voit en retranchant les équations (7), (8) 
(p. 15). 

La proposition physique correspondant a l’existence de la solu- 
tion (17) de (15) est la suivante : une double couche (magnétique) 
fermée dont le potentiel extérieur est nul a une densité constante. 
L’équation (7) du Chapitre Premier montre que le potentiel intérieur 
est constant et égal & 4a p('). 

2° Si au contraire V est un potentiel de simple couche, on voit 
en retranchant les équations (5), (6) (p. 15) que le premier membre 
de (18) est égal a = La solution p,(M) de (18) représente alors la 


densité d’une simple couche pour laquelle on a: 


dv 


re 


Le probleme physique ainsi résolu est celui qui consiste a trouver le 
potentiel intérieur d’une couche électrique en équilibre sur une surface 
conductrice. On ne connait ce potentiel, comme p,, qu’a un facteur 
constant pres qui sera déterminé par l’équation 


i (NAS 
S 


si on se donne la masse totale m de la charge. D’ailleurs ce potentiel 
est constant, d’apres le n° 4, a l’intérieur de S. 


(')Tl nest pas nécessaire dans cet alinéa (ni dans l'autre alinéa ot nous 
renvoyons 4 la note actuelle) de préciser pour l’expression «potentiel intérieur 
(extérieur)» sil s’agit du potentiel en un point quelconque intérieur (extéricur) 
a S ou de la limite de ce potentiel sur le bord intérieur (extérieur) de S. En 
effet : 1° Si V; est constant, les valeurs de V a J’intérieur de S n’ayant ni 
maximum ni minimum seront aussi égales 4 cette constante; 2° Si V. =0, on 
peut employer la formule de Green (14) et on en conclut comme au n° 4 que 
V est constante a l’extérieur de S. Comme V, = 0, cette constante sera d’ailleurs 
nulle. Donec si Ve = 0, V est nul partout a lextérieur. 
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L’existence des solutions de ces deux problemes physiques étant 
établie, il y a lieu d’en connaitre le nombre. Pour le second probleme, 
on voit qu’au point de vue physique il ne doit y avoir qu’une 
seule solution indépendante. Car une charge électrique m (constante 
arbitraire) s’étale d’une seule maniére sur une surface. On démontre ce 
fait analytiquement en prouvant que |’équation (18) n’a qu’une solution. 
Supposons qu'il y ait deux solutions linéairement indépendantes 
pi(M) et po(M). L’expression p!(M) = api(M) + 8p2(M) sera aussi une 
solution, et l’on peut choisir les constantes a et @ de facon que 


(19) [eo dS =am, + Bm2 = 0. 
S 


Si V’ est le potentiel de la simple couche de densité p’, on voit 
comme plus haut que le résultat de la substitution de p’ a p dans (18) 
peut s’écrire 


dv' 
dn; 


(20) =0. 


Alors en appliquant a V’ la formule de Green (13), on en conclut 
comme au n°4, que V’ est constant, a l’intérieur de S. Le potentiel 
d’une simple couche étant partout continu, on a donc V4 = constante. 

Or intégrons sur S la formule obtenue en remplacant p, V, par 
p’, V’ dans Videntité (5) de la p. 15. On obtiendra en tenant compte 
de (19), (20) 

av! 


dS = 0. 
on” u 


Appliquons maintenant la formule de Green (14) ot ’on remplace V 
par V’. On aura encore dans le second membre A?V!/ = 0 et puisque 
vi est constant, le premier terme du second membre deviendra 


Ov’ ov' 
/ = 
[ ¥ wan g One asi 
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On en conclut que le premier membre est nul et par suite que 
v’ est aussi constant a l’extérieur de S. De sorte qu’on a aussi 


Ov’ 


One = 


(20bis) 


et en utilisant la formule (5) de la p. 15, ainsi que (20) et (20S), on 
voit que la densité p’(M) est nulle identiquement, c’est-a-dire que 


a pi(M) + 8 p2(M) = 0 


et les densités p; et pp ne peuvent etre indépendantes. 

L’équation (18) n’ayant qu’une solution, l’équation associée (15) 
n’en aura aussi qu’une seule a un facteur constant prés. 

Remarquons que les équations de Fredholm expriment tout aussi 
bien ces problemes pour le cas de plusieurs surfaces distinctes 
extérieures les unes aux autres, en particulier le probleme qui consiste 
a étaler sur plusieurs surfaces conductrices en présence les unes des 
autres, des charges indépendantes, et on peut démontrer que, comme 
il est évident au point de vue physique, le nombre de solutions 
indépendantes correspondantes a \ = —1 est égal au nombre de 
surfaces. Par contre, les résultats qui viennent d’étre établis en 2° 
subsistent sans modification si D est le domaine compris entre une 
surface extérieure et une ou plusieurs surfaces intérieures a la premiere, 
le potentiel V étant encore visiblement constant dans les parties de 
Vespace intérieures a ces surfaces. 


6. — La valeur \ = +1 n’est pas une constante caractéristique des 
équations (11), (12). Pour le prouver, nous montrerons comme plus 
haut que l’équation homogene (18) correspondant a (12), par exemple, 
n’a pas de solutions non nulles pour A= +1. Cette valeur donnerait a 
Véquation homogeéne la forme 


(21) =0, 
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qu’on obtient comme plus haut en appelant V le potentiel de la simple 
couche de densité p sur S et en combinant les équations (5), (6) de la 
p. 15. 

On en conclut comme au n°4, au moyen de la formule de 
Green (14), que V est constante a l’extérieur de S, donc nulle, puisque 
V s’évanouit 4 Vinfini. Mais un potentiel de simple couche est continu 
a travers la surface : donc V; = 0. Et par conséquent — en appliquant 
cette fois la formule de Green (13) — on a V =0 partout. Donec 


dV dV ions 7 
ee ae 0 et d’aprés (5), p. 15 p(M) =0. 
7. — Nous pouvons maintenant formuler les solutions de nos 


problemes. Indiquons d’abord les résultats suivants qui se déduisent 
tout de suite des formules du n° 4 comme dans la note (*), p. 120. 

Il ne peut y avoir qu’une seule fonction V harmonique a l’intérieur 
de S, qui prenne des valeurs données sur le bord, — ou harmonique 
a Vextérieur de 8, qui prenne des valeurs données sur le bord et qui 
s’évanouisse a linfini, — ou encore harmonique a l’extérieur de §, telle 
que 7 prenne des valeurs données sur le bord, et que V s’évanouisse 
a Vinfini. 

Une fonction V harmonique a l’intérieur de S, telle que oy prenne 
des valeurs données, n’est pas unique, mais elle est définie 4 une 
constante arbitraire additive pres. Elle n’existe que si la condition 


| ~ dS = 0 est remplie, puisqu’on a pour toute fonction V 
g ON; 


(22) // AVardydz+ [ as =o. 
D gs Onj 


Proposons-nous de trouver une fonction harmonique a |’intérieur 
ou l’extérieur de S et vérifiant l'une des quatre conditions : 
go WV _ dV 


=0, 4° 
dn; , dne 


1° V;=0, 2° Ve=0, =0, 
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relatives aux problemes intérieur et extérieur respectivement. 

Pour les conditions 1° et 4° nous avons la solution évidente V = 0. 
C’est la seule d’aprés ce qui précéde. 

Pour 2°, l’équation intégrale correspondante (16) donne une solution 
a savoir le potentiel d’une double couche constante. Mais nous venons 
de voir que ce potentiel s’évanouit identiquement a l’extérieur de S. 

Pour 3°, nous avons le potentiel d’une couche d’électricité en 
équilibre, c’est-a-dire une constante a l’intérieur de S comme nous 
Vavons observé. 

Considérons maintenant les quatre conditions 


[oVese 22 Ves... 33° ch = Ae ON fonction donnée = f(P). 


On, One 

Les équations intégrales correspondantes sont (11) pour 1°, 2°, 
(12) pour 3°, 4°; avec A = +1 pour 1°, 4°, avec \} = —-1 pour 2°, 3°, 
comme nous l’avons établi au n° 9, p. 26. Puisque 4 = +1 n’est pas une 
constante caractéristique, l’équation de Fredholm nous donne pour 
1° et 4°, des solutions qui sont uniques (n° 1, p. 116). 

Elle ne donne des solutions de 2° et 3° que si la fonction donnée f(M) 
satisfait a certaines conditions de la forme (5), n° 1. Pour pouvoir 
les écrire, il faut former l’équation homogeéne associée correspondante 
et déterminer ses solutions. On a ici a prendre \ = —1. Or nous 
savons (n° 5) que pour \ = —1, les équations homogénes correspondant 
a (11), (12) n’ont qu’un seul couple de solutions : p = constante et 
p = pi(M), densité d’une couche d’électricité en équilibre, qu’on peut 
multiplier par un facteur constant de facon a donner a la couche une 
masse déterminée. 

Donc pour que nous ayons une solution de 2°, il faut que 


(23) [ f(P) p(P) dS = 0 


ou p = y(M) représente la densité d’une couche d’électricité en équilibre 
de masse égale a l’unité, par exemple. 
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Pour que dans le cas 3° notre méthode nous fournisse une solution, 
il faut de méme que 


(24) [ f(P) dS =0. 


Nous aurions pu écrire de suite cette derniére condition grace a la 
formule (22) qui peut s’écrire sous la forme de Gauss 


[ au dS = 4m x masse intérieure = 0. 
g dnj 

Ceci nous montre que la condition (24) est nécessaire pour que le 
probleme 3° ait une solution. Au contraire quand la condition (23) 
n’est pas remplie, nous savons seulement que la méthode de Fredholm 
ne donne pas de solution dans le cas 2°. Cela ne veut pas dire qu’une 
solution n’existe pas. Tout ce que nous savons c’est que cette solution 
ne peut étre représentée par un potentiel de double couche. 


8. — Les solutions précédentes sont exprimables sous forme de 
séries d’intégrales comme nous l’avons démontré dans le Deuxieme 
Chapitre |formules (47), (40), (41), (43)|. On a démontré qu’elles 
s’expriment aussi en séries de fonctions fondamentales, ainsi que les 
solutions des cas généraux (\ quelconque) des équations (11) et (12) ('). 
Nous nous bornons a établir les résultats suivants : 

1° Les constantes caractéristiques des équations associées (11), (12) 
sont réelles ; 

2° Elles sont des péles simples de la résolvante (7) ; 

3° La premiere est \=-—1; les autres ont des modules plus grands 
que l’unité. 

9. — Supposons qu’il y ait un pole complexe Ay = a +i. Alors 
léquation homogeéne correspondante a (12) a pour cette valeur de \ une 


(') Poincar#, Acta Mathematica, t. XX, 1897; PLEMELJ, STEKLOFF, voir la 
bibliographie, p. 175. 

(7) On remarquera que ces deux propriétés ont lieu bien que le noyau ne soit 
pas symétrique (comparer avec les n°*31, p. 90 et 34, p. 95). 
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solution non nulle p. Il y a un potentiel de simple couche correspondant 
V=V1+iVe2 qui, d’aprés le n° 9, p. 27, satisfait a l’équation 
oy aw w) | 19 | | W) a 


dne dn; dne dn; 


Car cette équation correspond a l’équation (12), lorsque nous y 
faisons 
gi(M) = 0, A= o- 
Séparons les parties réelles et imaginaires. 
Nous avons 


dV, dvi dV,  dvy dV9 dV9 
(26) dne dn, ‘ ES | me ES =| : 

dV2 dV2 Sag oe ] s oe i ad 

dne dn; dne dn; dne dn; 


Utilisons maintenant les formules de Green 


(27) i Vi ve Vo oN) igs ay, (V, AV2 — V2 AV) da dy dz =0 
S On; On; D 


(28) [ Vi em Vo ek ds II (V1, AV2 — V2 AV) dx dy dz = 0 
JS One One D 


en remarquant qu’ici les deux intégrales triples sont nulles puisque 
V1 et V2 sont harmoniques. 

Si nous multiplions les équations (26) par V2 et V; respectivement, 
si nous retranchons la seconde de la premiére, et si nous intégrons le 
résultat sur la surface S, nous aurons a cause de la formule de Green 
Véquation 


B [wie [wwe [ug [vite =. 
S dn; S dine S dn; S dne 


De méme en multipliant par V; et V2, et en ajoutant apres avoir 
intégré sur S$, nous avons 


dV, dV2 
(1 + a) if Vi Ais dS [v2 de as| 


dV, dVo2 
=(1 — d dS}. 
oe) i VA te [ dns s 
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Done pourvu que nous n’ayons pas 3 =0, puisque 


l+a ame 
l-a 
il faut que 
dV, dVo 
(29) [» Hig dS [ve iia dS = 0, 
dV, dV2 - 
(30) fu a ds + [Ve is dS = 


Or chacune des deux intégrales de (29), est supérieure ou au moins 
égale a zéro, comme nous voyons en faisant AV = 0 dans la formule 
de Green (14). Il faut donc que chacune s’évanouisse. On raisonnerait 
de meme avec |’équation (30). D’ot il résulte, en se servant encore des 
mémes formules de Green, que V est constante dans tout l’espace. Et 
ce résultat exigerait que p=0 d’aprés l’équation (5) (Chap. I). 

Donec 6 =0, et il ne peut y avoir une constante caractéristique 
complexe. 


10. — Supposons maintenant que A, soit un pole de degré m 
de la résolvante. La valeur \; sera (n° 1) un pole de degré m de la 
solution p de (12) (sauf pour des formes particuliéres de f;(M); mais 
cette fonction est arbitraire). Nous aurons donc si m> 1 


— _ pi(M) p2(M) R(M) 
p(M) = CBRL +7 re 


ot p1(M), p2(M) sont indépendantes de , et R(M) n’a pas un pole pour 
A=A\1. 
Substituons cette expression dans |’équation (12) mise sous la forme 


cos wy 


2 


cos w 


re 


(31) 2 p(M) + (A—) o(P) 8" as +4 i p(P) 8" as = f,(M) 


et égalons a zéro les coefficients de (A— \4)~™, et (A—A,)7™*1. 
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Nous obtiendrons 


cos 
2m pi(M) + At f pi(P) a dS = 0, 


2m p2(M) + Ax [ p2(P) we dS = — [ pi(P) oe ds. 


Considérons les simples couches p;(M), p2(M). Les potentiels V, et Vo 
de ces couches satisferont aux équations 


dV, dVi_ dVi | dVil _ 
2) dng dry” a Ee =| =0, 

dVo  dVa dVo . dVa dV; dVy 
33 r = : 
oe dre dnj a ES =| EB ?: =| 


Multiplions (32) par V2 et (33) par Vj ,, retranchons les deux 
équations et intégrons le résultat. On aura d’aprés (27), (28) 


[ue YY ag 4 [» ~ dS = 0. 


Multiplions ensuite (32) par V; et intégrons 


(+d) [vi a a0 [v1 avi dS = 0. 
Sg Maine dnj 


Ces deux derniéres équations sont distinctes : il s’ensuit donc 


[ugee- [vu Gees=o 
dne 


ce qui conduit a p;(M) =0 comme dans le dernier numéro. 

Il en résulte qu’on ne peut pas avoir un pole dont le degré est plus 
grand que l’unité. On voit facilement que ce qui précede ne s’applique 
pas au cas ol m=1. 


11. — Démontrons maintenant que les modules des constantes 
caractéristiques autres que \ = —1 sont plus grands que lunité. On 
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le voit immédiatement en multipliant l’équation (25) par V et en 
Vintégrant sur la surface. I] vient 


pul [Neos S| 
OT vs. Lys. = 3) 


Rappelons que d’apres (13), (14) 


[v® aca [v%as0, 
gs dn gs dne 


nous voyons bien d’aprés (34), que 


(34) 


|Ao| = 1. 


12. — M. Fredholm a montré qu’on peut appliquer sa méthode 
au cas d’une ou de plusieurs surfaces distinctes (1). 


M 
S2 


Si 
Fig. 3. 
Prenons comme exemple le cas de deux surfaces S$), So. (fig. 3). 


Etant donnée une série de valeurs sur ces deux surfaces, il s’agit 
de déterminer une fonction harmonique a l’extérieur de S$; et de So, 


(1)M. Picard a étudié ce probleme dans son cours de 1906. Voir H. B. 
Heywood. These, loc. cit., p. 63-69. 
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dont la dérivée normale prend sur ces surfaces les valeurs données. 
On considére que les points M et P prennent toutes les positions sur 
S; et Sp. La suite des valeurs données sera simplement une fonction 
de M. Le probleme est exprimé par l’équation (12), ot r est la distance 
entre les deux points M et P qui peuvent étre soit tous deux a la fois 
sur l’une des surfaces S$; et So, soit l'un sur $;, ’autre sur So. Toutefois 
il y a lieu de compléter la théorie générale du «déterminant» D(,) 
dans le cas des surfaces multiconnexes. 


13. Probleme de la chaleur. — Considérons maintenant le 
probleme 3° du n°4, p. 17, qui consiste 4 trouver une fonction 
harmonique satisfaisant a la condition 


~ + p(M)V =hA(M) sur la surface. 


Il est exprimé par |’équation 


|" ! | as = hi(M) 


E 


(35) om p(M) — X )-A foe 


(voir le n° 9, Chap. I) ot nous avons les cas intérieur et extérieur pour 
A= +1 et —1 respectivement. 
Pour le cas intérieur, ’équation de Green (13) nous donne 


fh [CY (BY + 
p|\ dz Oy Oz 
ot. V est un potentiel correspondant a l’équation intégrale (25) rendue 
homogeéne. 

Done si p est toujours négatif, cette équation homogene n’a pas de 
solution non nulle : nous ne sommes pas dans le cas singulier, il y a 
une solution unique de l’équation avec second membre (35). C’est ce 
qui a lieu pour le probleme d’un corps en équilibre de température 
avec rayonnement, ou p est une constante négative (voir 5°, n°6, 
Chap. I). Si p est positif ou change de signe, nous ne pouvons rien 
dire, et il faut examiner chaque cas séparément. 


dx dy dz = Lyz. as= | paveas 
S 
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Pour le cas extérieur, nous pouvons dire que nous ne sommes pas 
dans le cas singulier (autrement dit qu’il y a une solution unique), 
si p est essentiellement positif. 


14. Fonctions de Green. — Les paragraphes précédents vont 
nous permettre de construire des fonctions de Green G(M,P) (n° 8, 
p. 21) pour les conditions suivantes : 


1° G=0 sur la surface, intérieur et extérieur; 


xe 


— =0 sur la surface, extérieur seulement; 
0 
8 +p(M). -G =0 
intérieur si p(M) est essentiellement négatif; 


extérieur » » » positif. 


Pour les autres cas de 3° (p quelconque) nous ne pourrons rien 
dire (') sauf qwil y aura en général une fonction de Green. 

En effet, en posant G(M,P) = , —qw, tout revient a trouver une 
fonction @ harmonique dans un domaine D ou D’ limité par S$, et 
satisfaisant a l'une des conditions 


O O : , 
=, We =, = ei saat + p(M)w = fonction donnée; 
One On 


problemes qui ont été traités au n° 7, p. 14. 


WwW , - f 
5 sont données, nous savons que la solution existe et est 
Ne 


unique. Si we est donnée, on obtient par inversion (?), en ramenant 


Si w;, OU 


(1) On n’a pas encore approfondi cette question. 
(7) Il suffit pour cela d’utiliser la remarque de Lord Kelvin suivant laquelle si 
V(x,y,2) est une fonction harmonique il en est de méme de la fonction 


1 Vv 4a Yy 4 
[p24 y2 + 22" xe? + y2 + 22” xe? + y2 + 22? x* + y2 4+ 22 . 
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au cas intérieur, une solution a qui ne s’annule pas a l’infini. Enfin si 


6) 2 : : ; 
a. +p(M)a est donné, nous savons que la solution existe et est unique 
n 


pour les deux cas mentionnés plus haut au 3°. Si p est d’un signe 
quelconque, nous savons seulement que |’équation (35), n’est singuliére 
que pour des valeurs particulieres de \ qui en général sont différentes 
de +1. De sorte qu’en général, le probleme a une solution unique. 

En ce qui concerne la condition 2° pour le cas intérieur, si nous 
écrivions 


(36) G=--2 


r 


il faudrait que nous eussions une fonction harmonique a, qui satisfasse 
a la condition 


Om, 064 sur la surface 
On; 7 On; r : 


Or c’est impossible puisque, contrairement a la condition (24) du 
n° 7, ’expression 


i @ (=) dS=47r nest pas nulle. 
S 


ony \r 
Cependant si nous remplacons la condition 2°par la condition 


aG An 
40 2 = — 
st) On; 5 


ot. S désigne l’aire de la surface, nous aurons a trouver une fonction 


N 


@ harmonique a Vintérieur de S, et satisfaisant a 


Ow O1 4n 
On; On. ) 


sur la surface, 


ce qui est possible, puisque 


amas = [ p (5) as [ysqn te Zoe 
s On; g On; r 3 8 
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15. Les problémes généralisés. — I] s’agit ici de chercher (n° 5, 
p. 18) des solutions de 


(38) AV = y(z, Yy, z) 


qui satisfont a diverses conditions sur la frontiere S d’un domaine D. 
Nous avons déja traité au meme n°5 le cas intérieur pour la 
condition 
V=f(M)_ sur la surface S. 


Le probleme étant ramené au probleme de Dirichlet, a une solution 
unique. 

Pour le cas extérieur, en ramenant par une inversion au cas 
précédent, nous aurons également une solution, mais elle ne s’évanouit 
pas a l’infini. 

Cherchons ensuite la solution intérieure pour la condition 
av 


(39) on 


f(M). 
Nous utiliserons la fonction de Green G, telle que 


@ 
3 Me (ta, 40) 


Ecrivons l 
Vi=-— G M,P) y(P) dup. 
1 4 i ( ) ) ( ) P 


On a, d’aprés (13), p. 23 
AV, = p(x, y, Zz) 


et sur la surface 
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Ecrivons ensuite 


VSN: 
Nous aurons 
Mie S0: 
et sur la surface Vo , 
ae iM g 70: 


Nous aurons donc |d’aprés la condition (24) du n° 7] une solution 
de notre probleme pourvu que 


[ sones +0 =0 


ou 
(41) i f(M) dS + iE y(P) dup = 0, 


condition que fournirait directement la formule (22), p. 124. 

Le cas extérieur ne présente pas de difficulté, puisque la fonction 
de Green correspondante existe (n° 14, 2°). 

Nous arrivons a la condition mixte 


(42) Ad +p(M)V = f(M) sur la surface. 

Pour ce probleme, il y aura une solution pour le cas intérieur 
si p(M) est essentiellement négatif, et pour le cas extérieur si p(M) 
est essentiellement positif. En effet, dans ces deux cas, nous pouvons 
opérer comme précédemment, G étant soumis non a la condition (39), 
mais a 3° du n° 14 et V2 étant une fonction harmonique vérifiant (42). 
Les fonctions G, V2 existent et sont uniques d’aprés les n® respectifs 
14 et 13. Nous pouvons dire aussi que si p(M) n’a pas un signe 
constant, il y aura en général une solution, et nous pouvons prédire 
que pour certains cas isolés, il n’y aura de solution que si f(M) satisfait 
a une certaine condition. 
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III. — PROBLEMES RELATIFS A L’EQUATION AV =R(«,y,2)V 
ET A DES EQUATIONS ANALOGUES 


16. — Nous avons posé ces problemes au n° 11, p. 28 du Chap. I. 
Pour commencer, envisageons |’équation 


(43) AV = R(az, y, z)V 

Soit d’abord a trouver une solution analytique de (43) a l’intérieur 
de S, satisfaisant a la condition 
(44) V=0_ sur la surface. 


Nous avons vu (n° 13, Chap. I, p. 30) que cette question revient a 
la solution de ’équation intégrale homogéne 


(45) V(M) =-2 i 6(M, P) R(P) V(P) dup 


pour la valeur \=1('). 

Il n’y aura de solution non nulle 4 notre probleme que si \ = 1 
est une des constantes caractéristiques, autrement dit il n’y aura 
de solution que pour certaines formes particuliéres de la fonction 
R(z,y,z). Nous allons démontrer qu’il n’y a pas de solution si R(z, y, z) 
est toujours positif dans D (). 

Remplacgons V dans la formule de Green (13) par la solution 
cherchée de (43). Il vient 


wo fff |e) + (ae) +(e) [eae 
— f [f Rene V? aeay de 


(*) Le noyau G(M, P) R(P) devient infini quand M, P coincident, mais seulement 


1 
comme G et par conséquent comme ;) nous sommes donc dans les conditions 


étudiées au n° 16, p. 67, et note A, BA 155. 
‘6 ) Dans ce cas le noyau sera de la forme étudiée au n° 45, p. 114, puisque 
G(M, P) est symétrique. 
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d’ot il est évident qu’une solution non nulle V n’existe pas lorsque la 
fonction R(z,y,z) reste positive. 
Remplagons maintenant la condition (44) par 


(47) V = fonction donnée = f(M) sur S. 


Nous avons vu (Chap. I, n° 13) qu’il faut : 

1° Trouver une fonction harmonique v(M) qui satisfait a la 
condition (47). 

2° Trouver une solution, V; de l’équation non-homogéne 


(48) V(M) = = [ G(M,P) R(P) V(P) dep + o(M) 


(voir n° 13) ot dwp est un élément de volume. Cette fonction V; sera 
alors la solution désirée. 

Donc il y aura une solution unique pourvu que la solution pour la 
condition (44) n’existe pas. 

En particulier, nous aurons certainement une solution unique si 
R(z,y,z) est essentiellement positif. 

Considérons maintenant le cas extérieur pour la condition (44). 
Il est évident que cette solution est donnée de la méme facon par 
Véquation (45) pour un domaine extéricur D’. D’aprés la formule de 
Green correspondante (14), on voit comme plus haut qu’elle n’existe 
pas pour R constamment positif. 

Pour la condition (47), il y a une nouvelle difficulté. Si nous 
cherchons a résoudre le probleme par |’équation (48), il faut avoir une 
fonction v(M), harmonique a l’extérieur de la surface S$ et qui satisfait 
a la condition (47); ce n’est possible avec une fonction v nulle a Vinfini 
que si (n° 7, p. 125, cas 2°) 


(49) [ve seas =o. 
Posons 


(50) F= [ s(n sonyas. 
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En général, la quantité F est différente de zéro. 
Cherchons alors une fonction v, harmonique a l’extérieur de S et 
qui satisfait a la condition 


v(M) = f(M) -F 


sur 8, ce qui est toujours possible (remarquer que [ y(M)dS = 1, 
S 


p. 125). 
La solution sera donnée par l’équation 


(51) v(M) = = i, G(M,P) R(P) V(P) dwp + o(M) + F. 


17. — Envisageons maintenant la condition de Neumann 


Ov 


(52) — = 9(M) 


sur S, en commengcant par le cas particulier 
(52>i#) = =0 


sur S, ot l’on cherche une fonction extérieure. 

La fonction de Green relative a cette méme condition (52>!) 
existe d’aprés le n° 14 (cas 2°, extérieur). Dés lors la solution vérifie 
l’équation intégrale homogéne (45), ot G(M,P) est la fonction de Green 
propre a ce cas et ot l’on prend \ =1. La solution n’existe donc que 
si \=1 est une constante caractéristique. 

La formule de Green nous montre encore ici que cette fonction 
n’existe pas si R est constamment positif. 

La solution extérieure pour la condition (52) vérifiera l’équation 
intégrale (48), ou v(M) est maintenant une fonction harmonique 
satisfaisant 4 (52) et obtenue sans difficulté (voir n° 7, cas 4°). 

C’est pour le cas intérieur quwil y a des difficultés. Commencons 
par la condition (52%), 
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On peut former (voir n° 14, 4°) une fonction de Green G(M,P) telle 
que 


0G _ An 


On; S 


(53) 


sur la surface. 
Nous pouvons ajouter a G(M,P) une fonction de P, C(P) telle que 


I [G(M,P) + C(P)] R(M) duyy = 0. 
Cela revient a supposer que 
(54) He 6(M, P) R(M) dum = 0, 


G ayant les propriétés d’une fonction de Green, sauf d’étre symétrique. 
Considérons maintenant |’équation 


(55) V(M) = -2 i. G(M, P) R(P) V(P) dup. 


Nous avons 


OV nN 4a 
On; 7 Ar D S 


(56) RPI VPs I R(P) V(P) dwp, 


et d’aprés la formule (13), p. 23 
AV = AR(M) V(M). 
Or multiplions l’équation (55) par R(M) dw, et effectuons l’intégration 
dans le domaine extérieur D’. Nous avons 


, 
iL R(M) V(M) dist = i | i G(M, P) R(M) dwxy| R(P) V(P) dwp = 0. 
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Donec 
wv _, 
On; =? 
et Véquation (55) nous donne bien la solution de notre probleme 


pour A= 1. Cette solution n’existe que si \ = 1 est une constante 
caractéristique du noyau 


= G(M.P)R(P). 


Passons maintenant a la solution pour la condition (52) (intérieure). 
Posons 


: g(M) dS = G. 
S 
Construisons maintenant la fonction 


G 1 


= — — th ‘ 
ee Poe era 


D étant le volume de la surface. 


Nous avons d’aprés la formule de Gauss (n° 7, p. 125) 


Ova 


Loa ds =G. 


Déterminons enfin une fonction v,;, harmonique a l’intérieur de S 
et telle que 


Ov} Ove 
— M — 


ce qui est possible (n° 7, cas 3°), puisque 


Ov, 
S On; 


ds =G-G=0. 


La solution de notre probleme sera alors la solution de l’équation 


(57) V(M) = -< E G(M, P) R(P) V(P) dep + 01 (M) + v2(M) 
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pour \= 1, G(M,P) désignant la fonction de Green utilisée dans (55). 

La solution existera lorsque la solution pour g(M) = 0 n’existe pas. 
On voit comme plus haut qu’il suffit que R soit essentiellement positif 
dans D. 


18. — Nous aborderons enfin le probleme de la chaleur ot la 
condition sur la surface S est la suivante 
Ov 
Aa tT POMV = ACM). 


Nous avons vu au n° 14 que nous pouvons en général construire 
une fonction de Green G, satisfaisant a la condition 


aa p(M)G =0 
sur la surface. 
Cette fonction nous permet de résoudre de suite le cas ot 


h(M) =0. 


Nous avons en effet l’équation homogene de Fredholm (45), ou G est 
remplacé par G, et il y aura une solution dans le cas ot \ = 1 est une 
constante caractéristique. 

Et de méme nous avons une équation analogue a (48) pour traiter 
le cas ott h(M) £ 0. 


19. — Nous avons déja indiqué (voir n° 14, p. 31) comment on 
traite les équations 


(58) AV = R(a, y, z) oN, 
Ot 
av 


Il s’agit de trouver une fonction V qui vérifie une de ces équations 
a Vintérieur d’une surface réguliere S, qui se réduit a une fonction 
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donnée a(z,y,z) pour t = 0, et qui satisfait a certaines conditions sur 
la surface, conditions de Dirichlet, de Neumann, de la chaleur. 

Nous commengons avec (58) par le probleme intérieur avec la 
condition 


(60) V=0 


sur la surface. 
Substituons dans (58) 


(61) Vee Y(x, y, 2) [= et y(M)]. 
Il vient 
(62) Ay = —AR(M) py. 


Nous avons déja vu que la solution de cette équation avec la 
condition (60) revient 4 l’équation homogéne (voir n° 16), 


(63) eM) = 7 f GM.P) RIP) o(P) den, 


ou G est la fonction de Green nulle sur S. 
Donc il y aura des déterminations de y pour certaines valeurs de A 
seulement : les constantes caractéristiques du noyau 


= G(M,P) R(P). 


Ce noyau a la forme d’une fonction symétrique G(M,P) multipliée 
par une fonction R de P seul. Si R est essentiellement positif, ce qui est 
vrai pour les applications physiques, ce noyau entre dans la catégorie 
considérée par M. E. Schmidt (voir Ch. II, n° 45, p. 114). Il se raméne 
a un noyau symétrique : ses constantes caractéristiques sont réelles et 
ce sont des poles simples de la résolvante (+). 


(") Pour le cas ot R est quelconque, voir T. Marry, Comptes Rendus, 
fév. 1910. 
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Or au moyen de la formule (13) obtenue au n° 8, p. 23, on peut 
écrire 


(64) Af G(M,P) £(P) dep = —4n FM); 


nous allons montrer qu’a toute fonction g(M) continue ainsi que ses 
dérivées premiéres et secondes et nulle sur S, on peut faire correspondre 
une fonction f(M), telle que 


(65) [ G0?) 1) dep = 9000, 


En comparant (65) avec (64) on voit en effet que si la fonction f(M) 
existe, elle vérifie 
Ag(M) = —47 f(M). 


Il faut montrer que si l’on substitue cette expression de f(M) dans 
le premier membre de (65), la fonction obtenue 


aM) =—7- [| GOP) Apg(P) dep 


Ar 
n’est autre que g(M). 
Or on a d’apres la formule générale (64) 


Aq(M) = Ag(M) ou A[q(M) — g(M)] =0. 


Donc la fonction q(M) —g(M) est une fonction harmonique a 
intérieur de S, qui s’annule sur S$ (n° 8, p. 24), et qui est bien, par 
conséquent (n° 7, p. 124), identiquement nulle. 

Il en résulte que toute fonction g(M) nulle sur S$ et continue avec 
ses dérivées secondes est de la forme exigée par la condition (8) du 
n° 2, p. 118. Elle peut donc étre développée en une série absolument 
et uniformément convergente de fonctions fondamentales 


y1(M), p2(M), a) ~yn(M), ain 
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relatives au noyau G(M,P).R(P). D’ailleurs ceci exige manifestement 
que le nombre de ces fonctions fondamentales soit infini. 

Rappelons-nous maintenant qu’outre la condition aux limites, on 
assujettit la solution V de (58) a une condition initiale, laquelle consiste 
a se donner l’expression a(z,y,z) de V pour t= 0. La fonction a est 
continue avec ses dérivées secondes et d’aprés la condition (60) écrite 
pour t = 0, elle s’annule sur S. Elle satisfait donc aux conditions 
exigées pour le développement en série de g(M). 

Développons de cette maniére la fonction initiale a(z, y, z) 


(66) a(M) = Ay y1(M) + Ag yo(M) +... + An gn(M) +... 


ou les A sont donnés par la méthode de Fourier (n° 9, p. 11). 
Une solution de notre probleme sera alors 


(67) V(M, t) = Aj 1 (M) eit ahs 2(M) er2t oes ay. n(M) eT rnt Sheree, 


série qui est convergente puisque (66) l’est absolument. 

Cette solution satisfait a toutes les conditions prescrites. 

1° Elle satisfait a l’équation (58), puisque chaque terme satisfait a 
cette équation. 

2° Elle a pour limite zéro sur la surface. 

3° Elle se réduit & a(M) pour t= 0. 


20. — Au point de vue physique (n° 13, § 3°, p. 29), nous 
savons que la solution est unique. En effet, la loi de variation de la 
température V d’un corps dont la surface est a une température nulle, 
ne peut étre indéterminée quand on connait la température initiale a. 
Nous pouvons prévoir aussi que toutes les constantes caractéristiques 
doivent étre positives : un corps dont la surface a constamment la 
température zéro aura une distribution de température qui tendra 
vers zéro. Or en se reportant a l’expression (67) de la température V, 
on voit qu’il ne peut en étre ainsi si les A; ne sont pas tous positifs. 

Mais ces deux faits sont susceptibles d’une démonstration mathé- 
matique. Nous avons supposé que la fonction R est positive. 
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La formule de Green (13), p. 119, nous donne 


dv\?  (dv\? (dy? a 
LA ae) + (35) + (22) for=— foams fered 
(68) D v y z S Ny D 


a rf R(P) y? dwp. 
D 


D’ot il est évident que \ est positif. 
Nous avons déja démontré que la solution est unique au n° 14, 
p. 34. 


21. — Nous avons donc traité le probleme intérieur relatif a 
Véquation (58) pour la condition V = 0 sur la surface. Les problemes 
intérieurs pour d’autres conditions, et les problemes extérieurs se 
traitent d’une maniére analogue. Cependant on rencontre certaines 
difficultés que nous allons indiquer. 

Considérons la condition 


(70) V = f(M) 


sur la surface S, la fonction V vérifiant (58) a Pintérieur de S. 
Il est toujours possible de trouver une fonction harmonique indé- 
pendante de t, w(M), qui prenne la valeur f(M) sur la surface. 
Posons 
V=vt+u. 


La fonction v satisfera aux équations suivantes 


Ov 
Av = R(M)— 
(71) v=0 sur la surface 


v=a(M)-—w(M) pour t=0. 
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Nous supposons pour que le probleme soit possible que la valeur 
de V pour t =0, soit a(M), satisfait 4 la condition (70), de sorte que 
a—w Ssatisfait a (71). On est ainsi ramené au probleme du n° 20. 

Donc v sera déterminée par la méthode précédente. 

Le probleme extérieur se traite d’une maniére identique : dans 
ce cas la fonction V s’évanouit a Vinfini; on suppose que sa valeur 
initiale, la fonction a, s’évanouit aussi a l’infini. 

Nous ne nous arréterons pas aux conditions 


Nous avons en effet démontré (n° 14, p. 132) existence en général 
dune fonction de Green pour la premiere de ces conditions. Elle nous 
permettrait de traiter ces cas d’une maniére tout a fait analogue a la 
précédente, pour les problemes intérieurs et extérieurs. 

Pour le cas extérieur les conditions 


av av 
= = g(M) 


2 — = eek 
) On : On 
sur la surface, se traitent par une analyse identique, mais il faut 
supposer que la fonction V et la fonction initiale a sont nulles a 
Vinfini. 

22. — Nous arrivons maintenant au cas intérieur pour la condition 


Ov 
an 0 
sur la surface [toujours pour l’équation (58)]. 
Nous pouvons former une fonction de Green telle que 
6) An 


— @(M,P) = — 
Ba oO? S 


sur la surface (voir n° 14, cas 4°). 
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En refaisant la discussion du n° 17 apres avoir effectué la substi- 
tution (61), nous obtenons l’équation intégrale analogue a (55) 


Xr 
(73) eM) = f GULP) RIP) oP) dup, 
dont les solutions satisfont aux équations 
Ag = —rAR¢, 


i, R(M) v(M) duyy =0 dans le domaine D 
D 


et 


0 
°? 0 sur la surface. 
On 


En raisonnant comme au n°19, on voit que nous pouvons 
développer une fonction g(M) suivant une série de fonctions fon- 
damentales relatives a l’équation (73), pourvu que 


= sur la surface 


et 
i: g(M) R(M) dung = 0. 
D 


Nous supposons, bien entendu, que la valeur initiale de V, a(M), 
satisfait & la premiere de ces conditions. Elle ne satisfait pas en général 
a la seconde. Posons 


| a(M) R(M) dwyy = A, | R(M) dwyy = B. 
D D 


Alors la fonction a(M) -= est développable en série de fonctions 
fondamentales et l’on a 


A 
a(M) =prA yi(M) + Ag yo(M) +...+AnYn(M)+.... 
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La solution sera 
A 
(74) V(M,t) = 5 +A ¢1(M) e184 Ao yo(M)e7™?4+...4+AnGn(M)e t+... 
Car cette expression satisfait a toutes les conditions. 
Nous arrivons enfin a la condition 


Ov 


Ecrivons 


: g(M) dS = C. 
S 
Et considérons la fonction 
Vit) = f ROP) dup, 
D r 
ou r est la distance MP. 
On a d’aprés (13), p. 23 


AV, = —4r R(M), 


[F ds =[ if R(P) i (=) top| dS = in [ RE) dwp = C'. 


On; 


Donc la fonction a 
Up = Ci (—4rt+ V1) 


est une solution de 


av 
(75) AV=R=- 


telle que 
ov dS =C. 


g On; 
Construisons une fonction harmonique U;, indépendante de ¢ et 


telle qu’on ait 


OU: OU 


On; = of ) On; 
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Cela est possible par la méthode déja exposée au n° 7, p. 125, 


cas 3°, puisqu’on a 
OU, 


=C-C=0. 


La fonction U; ne sera déterminée qu’a une constante arbitraire 
additive prés; mais ce fait n’a pas d’importance. 
Ecrivons, enfin, la solution sous la forme 


V=U0+ Ui +u. 


La fonction u sera une solution de (58) qui pour t=0 se réduit a 
la fonction connue 


C 
ug = a(M) — G Vi - U1 
et qui satisfait a la condition 
du | OUp OU, | 
ae re oe 


La fonction initiale up satisfait aussi & la condition 


Oug 
On; per’ 


pourvu qu’on ait bien entendu = = g(M). 


Donc il est possible de déterminer u par la méthode précédente 
et V est complétement connu ('). 


23. — On verra que |’équation des ondes 
av 


est susceptible d’un traitement presque identique. I] est nécessaire 
seulement d’indiquer en peu de mots quelles sont les différences. 


(*) On peut démontrer qu’é la premiére constante caractéristique correspond 
une seule fonction fondamentale, laquelle est essentiellement positive. Voir 
Boussinesq. Traité analytique de la chaleur. Heywood, Thése, p. 90. 
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Pour fixer les idées, proposons-nous ce probleme : un gaz parfait est 
renfermé dans une surface fixe : on provoque une petite vibration 
quelconque du gaz, et on demande le mouvement du gaz a un instant 
ultérieur. A l’équation (59), il faut ajouter les conditions 


ov ine dare 
(76) Pn 2 Su la surface S 
et 
(77) V=a(M), ~ =6(M), pour t=0, 


la fonction V étant le potentiel de vitesses. 


Un type de solution sera y(M) a At, ou 
sin 


(78) Ay = —d*R(M) y(M). 
Construisons une fonction de Green G(M,P) telle que 


OG _ An 


ony S 


sur S (n° 14, 4°) ; 4 cette fonction nous pouvons ajouter une constante 
arbitraire (c’est-a-dire une fonction arbitraire de P). 
La fonction y sera donnée par 


(79) eM) —# [ GIMP) RP) oP) dup = 0 


pour une valeur convenable de \?. En effet, on a bien (78). Pour 
satisfaire a (76), remarquons que nous avons 
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Or 
[ R(P) y(P) dwp = 2? i I: G(M,P) R(M) R(P)g(P) divp deers 
=f | [gon P) ROM) aon R(P) y(P) dup. 


, 6) 
Par suite, © sera nul pourvu que 
nm 


i G(M, P) R(M) dunt = 0. 


Nous réaliserons cette derniére condition en choisissant convena- 
blement la constante arbitraire mentionnée ci-dessus. 

Nous pourrions aussi montrer comme au n° 20 que les constantes 
caractéristiques de l’équation intégrale (79), qui sont des valeurs de )?, 
sont positives, de sorte que A est réelle; mais a chaque constante 
caractéristique correspondent deux valeurs de \, donc deux solutions 
de l’équation (59). Nous aurons de méme qu’au n° 22, comme solution, 


la série 
CO 


V= S (an cos Ant + bn sin Ant) yn(x, y, 2). 
1 


Nous satisferons aux conditions initiales en choisissant les a et les b 
de sorte que 


a(M) = Sean Yn) B(M) = sbi An Yn; 
1 1 


ce qui est possible pourvu que 


dag 28 
On” On 


sur S. 
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24. Equation générale elliptique. — Considérons maintenant 
Véquation 


(80) Bat t Gye to ag to ay TeV ah 

ou a,b,c, f sont des fonctions de xz, y (que nous supposons avoir des 
dérivées premieres et secondes), c’est-a-dire des fonctions d’un point M 
de l’aire renfermée par un contour plan. 

L’application de la méthode de Fredholm @ cette équation (80) a 
été faite par M. Emile Picard(!); nous allons tirer notre discussion 
de son mémoire (7). 

Nous cherchons une solution de (80) qui s’annule sur le contour C 
d’une aire plane A, et qui est continue a l’intérieur de cette aire. Si 
Yon tient compte de la forme (13), p. 23 et de la note(?), p. 22, on 
voit que notre solution sera donnée par |’équation 


~a if {° (0) Fe tHE) 5 + eV \ g(e.ysgsn) df 


(81) 
= v(z,y) = 5 | {feng (x,y; 7) dé dn, 


(ou G est une fonction de Green — relative a l’équation AV = 0 — qui 
s’annule sur C), expression que nous transformons en 


alag) (0G) 
(82) vin +z ff “ + oghy (€,n)d& dn = (2,9), 


en intégrant par parties. 


La fonction de Green est comparable, a notre point de vue, a log , 
ou r est la distance entre les points (x,y) et (€,7). Donec le noyau de 
Véquation (82) est comparable a = et, d’aprés le n° 14, p. 64, elle 
possede bien, en général, une solution V. 


(2) Egalement par M. Hilbert : Nachrichten zu Gottingen, 1904, Heft 3, p. 248. 
(7) Rend. Circ. Matem. Palermo, juillet 1906. 
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Pour achever la discussion, il faut montrer que cette fonction V 
nous permet de remonter de (82) a (81) et de (81) a (80). C’est-a-dire, 
il faut montrer que V a des dérivées premieres et secondes. 

Nous écrivons le noyau de (82), K(2,y;£,7), et nous effectuons sur 
cette équation une «itération» (voir formule (7), p. 40). 

Il vient 


V(,9) - / fi Kol, ys€,n) V(Evn) dé dn 


(83) 
= V(2,y) — ff Kays.) HGn) a6 do, 


Or y(a,y) est comparable avec un potentiel : elle a donc des 
dérivées secondes, puisque nous avons supposé que f a des dérivées 
premieres ; la fonction Ky est comparable avec logr; de sorte que nous 
n’avons a nous occuper que du dernier terme. Il s’agit donc de voir 
par exemple que l’intégrale 


[[ 2 venacen 


a des dérivées premiéres; mais elle s’écrit 


(84) = ff age ae dn, 


et cette intégrale a des dérivées qui restent finies sur le bord. 
Il faut maintenant examiner les dérivées secondes. I] est évident 
que 


Jf Sole. VE om) a dy 


a des dérivées secondes, puisque V(€,7) a des dérivées premiéres. 
L’expression (84) a des dérivées secondes pourvu que wy ait des dérivées 
secondes. 

Donec l’existence des dérivées est établie. 

Nous pouvons dire que notre probleme a en général une solution 
unique. Mais il y a certains cas ot il existe une solution du probleme 
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quand nous faisons f égale a zéro : pour ces cas, — c’est-a-dire pour 
ces mémes fonctions a, b, ec — d’aprés la théorie de Fredholm (n° 21, 
p. 70), le probleme ot! f n’est pas identiquement nulle n’a pas de 
solution (A moins que w(x,y) fat précisément nul). 


NOTE A 


ITERATION DES NOYAUX INFINIS DANS LE CAS DES 
INTEGRALES DOUBLES (') 


1. — Nous savons que la résolution d’une équation de Fredholm 


b 
(1) gays / y(t) K(s,t) dt = f(s) 


a 


de noyau K(s,t) peut se ramener & celle d’une nouvelle équation de noyau 


b 
ey = i) K(s,t1) K(t, t) dts, 


ott Ko(s,t) est le noyau itéré une fois. 

Nous avons vu (n° 14, p. 64) V’effet d’une pareille itération pour une 
intégrale simple dont le noyau est susceptible de prendre des valeurs 
infinies de la forme 


H(s,t 
Ks, 7) = (s,#) 
(s—t)* 
(avec a < 1 pour que les intégrales aient un sens) — H(s,t) restant fini 


— et nous avons constaté qu’aprés un nombre suffisant d’itérations on 
obtient un noyau qui reste fini. 

Etudions au méme point de vue le cas d’une intégrale double. Soit 
Véquation 


(2) Mm) - f fe) K(M, P) dsp = f(M) 


M et P étant deux points pris sur une surface et la ff étant étendue a 
cette surface. 
Lorsque K(M, P) devient infini, si ’on a, par exemple, 


H(M, P) 


() D’aprés une lecon professée par M. Hadamard, recueillie par M. Pérés, 
éleve de l’Ecole Normale Supérieure. 
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(avec a < 2 pour que la ff ait un sens), on parviendra encore — aprés un 
nombre suffisant d’itérations — a un noyau restant fini auquel la méthode 
de Fredholm s’applique. 


2. — Prenons par exemple a = 1; nous allons démontrer que le noyau 
itéré une fois — qui peut encore augmenter indéfiniment pour M et P 
voisins l’un de autre — n’est plus que logarithmiquement infini. On peut 
sans restreindre la généralité supposer que la surface est un morceau du 


plan des €, 7 (*). 
Le noyau itéré est alors 


2(M, P) -[[ MQ, Be ata 


(Q étant le point de coordonnées €, 7). Considérons un cercle C de 
centre M, de rayon proportionnel 4 MP, de rayon 2MP par exemple et 


contenant P dans son intérieur. 
Décomposons l’intégrale double en deux parties, l’une relative aux 
points Q intérieurs a cette circonférence, l’autre relative 4 Q extérieur. 
La premiére intégrale est plus petite que 


2 ff &an 
ee os 


x étant la limite supérieure supposée finie de H(M, P). 


Quant a 
/[x1a-0F 


un changement de variables par similitude amenant le rayon du cercle C a 
Vunité ne la modifie pas, elle reste donc finie quand M et P tendent lun 
vers l’autre et il en est de méme de notre premiere intégrale. 


(‘)En effet, lorsque M et P sont voisins ’un de l’autre, on peut, entre la 
portion de surface (supposée réguliére) qui les contient et une portion de plan, 
établir une correspondance telle que le rapport de deux éléments de surface qui 
se correspondent soit compris entre deux limites finies et qu’il en soit de méme 
pour le rapport entre la distance de deux points quelconques et celle de leurs 
homologues. 
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Dans la deuxiéme intégrale, on peut remplacer QP par MQ, car le 
rapport de ces deux quantités reste compris entre > et — : elle augmente 


dés lors indéfiniment au plus comme log MP, ainsi qu’on s’en assure 
immeédiatement en passant aux coordonnées polaires de pole M. C’est. bien 
le résultat annoncé. 

Une seconde itération (que l'on pourrait méme éviter) nous ramenera 
a un noyau fini. 


3. — La circonstance précédente se présente dans la résolution du 
probleme de Dirichlet a trois dimensions par la méthode de Fredholm. 
Dans le cas de deux dimensions, le noyau 


dlog — 
K(M, P)= : 7 eg neem) avec r= MP 
T Mp 1 - 
cos(r, 2) 


reste fini quand M tend vers P, car tend vers a y étant la 
courbure en P. 


Dans le cas de trois dimensions au contraire 


1 
qd _ cos(n, 7) 
dnp rr? 


: : 1 
tend vers l’infini comme -—. 


r 
C’est le cas étudié précédemment : deux itérations nous donneront un 


noyau fini auquel s’appliquera la méthode de Fredholm. 
H(M,P 
Le cas ot K= oo) et ol a4 1 s’étudie comme le cas ot. a = 1. 
On partage comme précédemment l’intégrale 


Ko(M, P) -{{[= MQe Soe iP) dé dn. 


L’intégrale étendue a l’intérieur du cercle est inférieure en valeur absolue a 


* [xara 
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et Vintégrale 


i, i) d€ dn 
c MQ? QPo 
1 
vaut ————.—— fois l’intégrale analogue étendue & la figure semblable telle 
(2MP)20—2 S z 5 


que le cercle soit de rayon 1, cette derniére intégrale étant une quantité 
finie fixe. 
L’intégrale étendue a l’extérieur du cercle se comporte aussi comme 


(@ipy22—2 quand MP et par suite le rayon du cercle tendent vers 0. 


Si donc a < 1, une seule itération nous aura amené a un noyau fini. Si 
a > 1, nous avons remplacé l’exposant a par l’exposant plus petit (puisque 
a <2), 2a—2. 

Il n’est pas bien difficile de voir qu’un nombre fini d’itérations nous 
ameénera dans ce cas & un noyau fini. 


4. — Il est clair que la méme méthode s’appliquerait sans modification 
pour le cas de domaines & un nombre quelconque de dimensions, toutes 
les fois que le noyau deviendrait infini comme l’inverse d’une certaine 
puissance de la distance entre les deux points M, P mobiles dans un pareil 
domaine. 

Ajoutons que l’intégrale par l’étude de laquelle nous venons de 
commencer, savoir l’intégrale plane 


dg dn 
MQ. QP’ 
conduit a une formule particuliérement simple lorsqu’on considére la 
différence de deux intégrales analogues, soit par exemple 


1 1 
(3) BS || (sa-op 0.08) ot 


(€, 7 étant toujours les ordonnées du point Q, pendant que M, P, P’ sont 
trois points fixes du plan). Une telle intégrale reste finie lorsqu’on |’étend 
au plan tout entier, et la méthode employée plus haut en fournit, dans 
ces conditions, la valeur exacte. Si, en effet, on limite tout d’abord l’aire 
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d’intégration a Vintérieur d’un cercle déterminé C de centre M, on pourra, 
d’une maniére analogue a ce qui a été fait tout 4 ’heure, remplacer 


I [. 390-0 
par 


@ [] star 


étendue & un cercle C’ homothétique & C par rapport 4 M, avec le rapport 
/ 


de similitude ie L’intégrale (3) étendue a C, peut donc se remplacer 


par l’intégrale (4) étendue 4 la couronne comprise entre C et C’. En 
faisant croitre indéfiniment le rayon de C, on trouve ainsi, pour l’intégrale 
relative a tout le plan, 
MP 
I= 27lo : 
© MP’ 
Bien entendu, il en résulte, pour l’intégrale 


h = | f acan (sia. QP wT ap) 


(ot M, M’, P, P’ sont quatre points fixes), étendue 4 tout le plan, la valeur 


MP 
ly = 27 log wp" 


Enfin, on pourrait opérer de méme en remplacant par 


1 
MQ.QP 
n’importe quelle fonction homogéne de degré —2 de MQ, QP, non 
susceptible de devenir infinie le long d’une ligne, par exemple, 


1 
a.»MQ?+b.MQ.QP +c. QP?’ 


oti la forme quadratique (a coefficients constants) qui figure au dénominateur 


ne s’annule pour aucune valeur positive du rapport —. 


QP 
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On pourrait également écrire des formules analogues dans l’espace a 
trois dimensions (la fonction homogéne devant alors étre de degré —3), 
etc. 


NOTE B 


PROPRIETES DE LA RESOLVANTE DE L’EQUATION 
DE FREDHOLM 


Nous allons donner ici une étude un peu plus approfondie des propriétés 
de l’équation de Fredholm 


b 
(1) ols) -A f K(s,t) elt) at = f(s), 


dans le cas d’un noyau K quelconque. On trouvera intérét & comparer 
les résultats obtenus avec ceux du cas symétrique (p. 89 — 114), qu’ils 
généralisent en partie. Nous nous bornerons a développer les idées 
essentielles sans nous occuper de donner les démonstrations rigoureuses, 
pour lesquelles nous renvoyons le lecteur aux mémoires cités Cy, 


1. Les noyaux orthogonaux. — On appelle ainsi deux noyaux 
Ki (s,t), Ko(s,t) qui satisfont aux deux relations 


(2) i GHA ae. 


e 
(3) 7 GAG da =o: 


Nous allons démontrer les deux propositions suivantes : 
Si nous écrivons 


(4) K(s,t) = Ki(s, t) + Ko(s, ¢), 
Ki (s,t) et Ko(s,t) étant orthogonauzr, nous aurons 


(5) (1°) K(s,t,A) = Ki(s,t, A) + Ka(s,t, A), 
(6) (2°) D(A) = Di(A) x Da(A), 


(1) H-B. HEywoopb. — Comptes Rendus, 25 novembre 1907; Thése, 1908; 
Journal de Mathématiques, octobre 1908. 
E. GOURSAT. — Comptes Rendus, octobre et novembre 1907; Annales de la 
Fac. d. Sc. de Toulouse, t. X, 1908. 
M. GourSAT démontre ces résultats par un noyau intégral borné. 


L’EQUATION DE FREDHOLM 162 


ou nous avons adopté la notation habituelle. 
Pour démontrer le premier résultat, souvenons-nous que nous avons 
(formule (46), p. 61) les deux relations 


b 
(7) Ki(s,t, X) <4 Ki(s,¢) _ | Ki(s,7) Ki(7,t, A) dr, 
oh 
(8) = | Ki(7, t) Ki(s, T, A) dr, 
a 
et les mémes équations avec l’indice 2 — (71), (81). 


Nous multiplions l’équation (8) par Koa(t,t’) et nous intégrons par 
rapport at: il vient, en tenant compte de (2) 


(9) | "Ki(s,t, d) Ko(t,t’) dt = 0. 


De méme 
b 
(10) Fi Ki (s, t, A) Ko(t’, s) ds = 0, 
a 


de sorte que Kj(s,t, A), Ko(s,t) sont orthogonaux. 
Pareillement Kj(s,t), Ko(s,t,) sont orthogonaux. 
Ajoutons les équations (7) et (71) : nous aurons 


[Ki(s, t, 4) + Ko(s,t, 4)] — K(s, t) 
b 
=X J [Kils,7) Ki(rst,2) + Ko(s,7) Ka(r.t,2)] dr 
b 
= | [Ki(s, 7) + Ko(s,7)] [Ki(z, t, A) + Ka(r, t, A)] dr, 


b 
= af K(s,7) [Ki (7, t, A) + Ka(r, t, r)| dr. 


Il en résulte que la fonction [Kj(s,t, A) + Ko(s,t, A)] est déterminée par 
les mémes équations fonctionnelles (46, p. 61), que la résolvante K(s,t, A), 
et puisque leur solution est unique (n° 5, p. 43), nous avons l’égalité (5). 
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Pour démontrer la seconde égalité (2°), nous nous servons de la 
formule (50>) du n° 17 (p. 68) 


b 
(11) < (log D(A)) = - | K(s, s, A) ds, 


b 
= 7 [Ki(s, 8, A) + Ko(s, 8, )| ds, 


[log Di (A) + log Da(A)], 


siaele 


py tlog[Pi) x Da()]}; 


c’est-a-dire 
D()) =Cx Di (A) x Do(A). 


Pour achever, il faut démontrer que la constante C est égale a l’unité. 
Il suffit de faire A = 0 


D(0) = 1= D1(0) x Da(0). 


Donc nous avons Végalité (6) (2°). 


2. — De la derniére proposition il résulte que si A, est une constante 
caractéristique de Kj(s,t) ou de Kg(s,t), elle sera aussi une constante 
caractéristique de K(s,t). Cela posé, nous allons démontrer la proposition 
suivante : 

Les solutions de 


b 
(12) vis) — Aa / Ki (s, t) yi(t) dt = 0 


et les solutions de 


b 
(13) (p2(s) —, / Ka(s, t) p(t) dt =0 


sont aussi des solutions de 


b 
(14) Bia= af Keven: 
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Les solutions de (14) sont des fonctions linéaires des solutions de (11) 
et de (12). 

On n’oubliera pas que la solution de (12) ou de (13) doit se réduire 
(n° 21, p. 70) a zéro si Ay n’est pas constante caractéristique de K; ou 
de Kg. 

Nous multiplions (12) par Ko(r,s) et nous intégrons : 


i Ko(r, s) pi(s) ds = »1 [ i/ Ka(r, s) Ki(s, t) as yi(t) dt = 0. 


Donc 
b 
y1(s) — af K(s, t) yi(t) dt 


b b 
=a(s)-» f KGva@ at—>y | Ra, Oi BESO 


a 


Il en résulte que y1(s) est une solution de (14). De méme yo(s) est 
aussi une solution de (14). 

Réciproquement, soit y(s) une solution quelconque de (14). 

Posons 


b 
‘ei i Ki (s,t) y(t) dt = 5,(s) 
p(s) = 01(s) + Oo(s). 


(15) 


Je vais démontrer que ®;(s), ®2(s) sont des solutions de (12) et de (13) 
respectivement. Des équations (14) et (15) on tire 


b 
(16) Do(s) = Aq | Ka(s,t) y(t) dt, 
d’ot il résulte 


[x0 $o(t) dt = »1 [ [ [m0 Ka(t, 7) dt| y(t) dr = 0, 
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et d’aprés (15), 


b b 
(5) —r f Kite t) 5,0) a= | Kis Oued =t. 


c’est-a-dire que ®)(s) est une solution de (12). De méme en nous servant 
de (14), nous pouvons démontrer que ®9(s) est une solution de (13). 


3. Partie d’un noyau relative 4 une constante caractéristique. 
— La résolvante K(s,t,A) du noyau K(s,t), dans le voisinage du point 
A = 1, ot Ay, est une constante caractéristique, peut s’écrire comme au 
n° 21, p. 70, sous la forme 


a pr(s, b) ; Pr—1(5, b) ’ 1 (s,¢) 
(17) K(s,t, d) = Ox = yr T ont = rl Tr sosteyk On = ») 
= x(s,t, A) + yo(s,t, A), 


pols, t, X) 


ou Yo(s,t, A) est une fonction continue de X. 
Nous appelons y(s,t, A) la partie de la résolvante relative a la constante 
caractéristique A,; pour \ = 0, nous avons 


K(s, t) = x(s, t, 0) + yo(s, t, 0), 


soit 
= x(s, t) + vo(s, t), 
= (s,t) (s,t) (s,t) 
rls, Pr-1(S,t pls, 
i) = ..- oo. 
x(s, t) Gi ca +...4+ 7 


Nous appelons x(s,t) la partie du noyau K(s,t) relative a Ay. 

Nous allons démontrer que les deux fonctions x(s,t), yo(s,t) sont or- 
thogonales. De plus, si A1, A2,..-,An,--- sont les constantes caractéristiques 
du noyau K(s,t), & chacun de ces nombres correspond une partie 


p(s,t), p@(s,t), ..., p™(s,t), .... 
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Ces parties sont orthogonales entre elles. Si la série 
H(s,t) = p(s, t) + op? (s,t) +... + pM(s, +... 
est uniformément convergente, nous n’aurons pas en général 
K(s,t) = H(s,t). 


En posant 
K(s,t) = H(s,t) + p(s, t), 


nous appellerons p(s, t) la partie relative a l’infini : elle est biorthogonale 
a toutes les autres parties. 

Nous n’entrerons pas dans les détails, qu’on trouvera dans les mémoires 
cités, de la démonstration de ces propriétés. Nous nous bornerons aux 
idées principales. 

Il nous sera utile d’établir d’abord quelques relations entre les fonctions 
y1(s, t), po(s, t).... 

Rappelons dans ce but la formule (11) : 


b 
(11) i K(s,8,A)ds = -< [log D(A)]. 
Si 1 est un zéro de D(A) d’ordre p, nous aurons 
D(A) = (Ar — A)PEQ), 
oul 


E(\1) £0. 


Donc q d 
{2 = 
Da [log D(A)] = res + D [log E(A)]. 


d 
L’expression — [log E(A)| reste finie lorsque tend vers \1, et nous 


tirons de (11) l’équation 


b 
(18) lim (41 =A) f K(s,3,)ds= +p. 


(21) 
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e (17) et (18) il résulte que 


b b b 
i (ss) ds =, / eels side 0. wu, / (yo(s, s) ds = 0, 
(19) a a a 


b 
| ile ee 


Reprenons d’autre part, les équations (46) du n°11, p. 61, sous la 
forme suivante : 


— K(s, t) + K(s,t, d) 


b b 
(20) =(-») KG Kae ae rf K(s,7)K(r,t, d) dr 


= (A— 21) [s, T, A) K(t,#) dr + 4 i: K(s,7, A) K(z, t) dr. 


Nous substituons ensuite l’expression de K(s,t, A) tirée de (17) dans 
les équations (20) et nous égalons non plus seulement les coefficients de 
(Ay — A)~", mais ceux de (Ay — A)~4, (Ay — A)7?,..., (Ar — A)7777, dans les 
trois membres. Ceci nous donne comme & la page 70, 


b b 
oe a. f CTC eer ee af Keeler ar a0: 


b b 
er-als,t)— ar | Rie ate ar=- f Kise iar 
(22) a a 


b b 
~r—1(s, t) — af K(5t) @-=1(6,7 dr = a K(T, t) Yr (s, 7) dr, 


a 


b b 
y1(s, t) -a f K(s,7) yi(7,t ay K(s, 7) ya(7, t) dr 
(23) a b 
yi(s,t) — a. f K(7, t) yi(s,T) --| K(r, t) ya(s, 7) dr, 
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et 
b b 
yo(s,t, A) — af K(s,7) go(r, t, A) dt = K(s, t) -f K(s,7) yi(7, t) dr 
(24) rb 
yo(s,t, A) — | K(r, t) Yo(s, 7, A) dr = K(s,t) -{ K(z, t) ¢1(s, 7) dr. (1) 


On transforme ces équations pour leur donner les formes suivantes 


u b 
(25) i K(s,7) gr (7, t) dr a K(r1, t) gr(s, 7) dt = Belay t) 


b b 
[ Keneatsear= [ K(7, t) Gr_i(s, 7) dr 
_ prilsit) , Palsyt) 
My ee 


b b 
[ Koneurnar= f K(z,t) g1(s, 7) dr 


en 2° _ gilsst) , valst) , _, volt) 
ce T ”? iF as Sel x 
= x(s,t), 
b b 
| K(s,1) go(r.t,2)dr => f K(r, t) yo(s, 7, A) dr 
(28) pls) 4 rls!) 


= —K(s,t) + yo(s,t, A) 4 


yy fee uy 
= vols, t, d) = yo(s, t, 0). () 


On multiplie ensuite l’équation 


(20) K(s,t,4) —K(s,t) = af K(s,7, A) K(7,t) dr = af K(7, t, A) K(s, 7) dr 


a a 


ey On peut écrire simplement ces formules en remplacant les premiers membres 


par Sy,.,Sy,-1)+++;5y9, au moyen de la notation introduite au n° 5, p. 43. 
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par y,(t,u) et on l’intégre par rapport A t. En se servant de (25) on a, 
apres un léger changement de notation, 


b 
(29) | Kena elrpar= | K(r3t,) golsy7) dr = SAS. 


Et de la méme maniére on obtient les autres équations 


b b 
| K(s, T, d) Pals, t) dr = ji K(r, t, X) Pr—1(8, T) dr 
a a 


ee _ Grailsyt) , velsyt) 
M—d Oy =A)?’ 
[ Kor vacser= [ K(7,t, A) g1(s, 7) dr 
(31) = ¢1(s,t) y2(s, t) r(s, t) 
=o Oese Orar 
=i. te ANS 


On se sert de nouveau de l’équation (17) pour séparer dans les équations 
1 


1 1 
(29), (30), (31) les coefficients de = Gake’?  Oeae et lon 


obtient ainsi une série d’équations qui se résument sous la forme suivante : 


b 
a Bee He=0, 


b 
(32) | Ym(S8,T) Yn(T, t) dr simtn>rt+l, 
: a Pmin=1(s,t), 
sim+n<rtl, 


b b 
(33) / San idee i ARO O, 


pour toute valeur de m £0. 
Nous sommes maintenant en mesure de prouver les propriétés énoncées. 
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4. — A l’aide de (33) on peut démontrer que les deux parties de 


K(s, t) = LS; t) ot yo(s, t) 
sont orthogonales. On a, pour toute valeur réguliére de A, J, 


b b 
(34) i Man \eoteiidr = | COIN ESO. 


Il suffit pour le voir de remplacer y(s,7,) dans cette expression par 
son développement 


e(s.7) |, euler) 


x(8)7) 0) = Tyr 7 oT ty = 


et de tenir compte de (33). 

En faisant A = pp = 0 dans (34), on voit bien que x(s,t), Yo(s,t) sont 
orthogonales. D’autre part en y remplagant ’ par 0, par A, nous avons, 
d’apreés (28), 


b 
~$o(9;t)-+ volat,d) =A i yo(s,7) volt; t, A) dr 
(35) a 


b 
— af po(T, t) yo(s,7, A) dr, 
a 


et, d’aprés (20°) et (17) 


b 
— x(s,t) + x(s,t, A) = | x(s,7) x(7, t, A) dr 
(36) 


b 
= | x(7, t) x(s, 7, A) dr. 


On peut donc dire que les fonctions yo(s, t) et x(s,t) sont orthogonales 
et ont respectivement pour résolvantes les fonctions yo(s, t, A), v(s,t, A). Il 
en résulte que nous pouvons appliquer les considérations du n° 1, et que, 
pour calculer les solutions de Véquation intégrale homogéne (14), nous 
pouvons considérer a part la partie x(s,t) du noyau K(s,t) relative a 1. 
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Le déterminant E(A) de y(s,t) est donné par la formule (11) qui 
devient ici 

aS 

dr 


: =P P 
[log E(A)| [ x(8,8,A) d8 = > = yy 
la seconde égalité provenant de (19). 
Par conséquent, 


E(A) = (A — Ax)? X const. 


Mais 
(0) =1, 
d’ou 
(37) F(A) = (1 a >) . 


C’est le facteur du déterminant D(A) de K(s,t) relatif a la partie 
x(s, t). 

On peut démontrer enfin que les parties du noyau relatives a4 deux 
constantes caractéristiques sont orthogonales. 

Soient, en effet, y1(s,t), x2(s,t) les parties du noyau K(s,t) relatives 
& A1, Ag; soient y1(s,t,r), x2(s,t,A) les parties correspondantes de la 
résolvante. On peut poser 


K(s,t, Ad) 2 x1(s, t, A) + x2(s,t, ) ae wo(s, t, Xd), 


et, en appliquant les formules (34) — ot lon a permuté ju et A’ — pour 
w=0 


b 
| x1(8,7) [x2(7,#,A) + do(7, t, A)] dt 


b 
-| x1(7, #) [xa(s, T,A) + wo0(s,7, d)] dr = 0. 


La fonction ~o0(s,t,) reste finie lorsque » s’approche de 2; donc, 
en se souvenant de la forme de x2(s,t,A), en l’introduisant dans l’égalité 
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précédente et annulant les termes en ., il restera 


b b 
/ x1(s,7) x2(T, t, A) da = xa1(T, t) xa(s, 7, A) dr = 0. 
a a 


Faisant 4 = 0, on aura enfin 


b b 
(35) / x1(8,7) X2(T, t) dr =) x1(T, t) xo(s,7) dr = 0. 


a 


Nous avons donc démontré que yi(s,t), x2(s,t) sont orthogonales. 


5. Fonctions principales. — Envisageons maintenant la fonction 
yr(s,t). L’équation (21) indique qu’elle est une solution de chacune des 
équations associées de Fredholm sans second membre pour A = Aj. Celles-ci 
ont le méme nombre q de solutions linéairement indépendantes. Donc, 
considérée comme une fonction de s, y,(s,t) est la somme de q fonctions 
linéairement indépendantes au plus, les coefficients étant certaines fonctions 
de t. De méme, considérée comme une fonction de t, y,(s,t) est la somme 
de q fonctions indépendantes au plus. Nous pouvons écrire y,(s,t) sous la 
forme 


vr(s,t) = >> x(s) melt). 
k=1 


D’aprés l’équation (22) on voit que y,—1(s,t) est une solution de cha- 
cune de deux équations associées de Fredholm, avec second membre, pour 
» = 1. Les conditions pour l’existence des solutions sont nécessairement 
remplies, d’aprés (32). La forme des équations (22) démontre que y,_1(s,t) 
est aussi la somme d’un nombre fini de produits : 


k=m 
Gr—1(8,t) = S~ Ox(8) pr(t). 
k=1 
Nous examinons de proche en proche toutes les fonctions y,—2(s,t), 


Yr—3(s,t),...,¢1(s,t), et nous trouvons enfin que yj(s,t) est aussi la 
somme d’un nombre fini de produits : 


k=n 
(36) gi(s,t) = >— pe(s) vet). 
k=]. 
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Substituons cette expression dans une des équations (32) : 


b 
/ y1(s, T) yi(T, t) dr = g1(s,t). 


Il en résulte 


b 
YE vets v0 | es) velr) ar =F val) oul. 


k k 
On peut s’arranger pour que les fonctions y1(s),...,Qn(s) soient 
linéairement indépendantes, et que les fonctions a (t),...,W%n(t) soient 


aussi indépendantes. Dans ces conditions, égalons les coefficients de 
pr(s), W(t) dans les deux membres. Ceci nous donne les relations 


= si ) 
om awooufet 0 


D’ot, en intégrant l’expression (36) pour s =t 


b 
/ yi(s,s)ds =n. 


Rappelons une des équations (19), 


b 
/ yi(s, s) ds = p, 
a 


ou p est le degré du zéro A; de D()). 
Nous obtenons n = p, 


P 


(38) yi(s,t) = a pr(s) Wx(t). 


1 


Nous appellerons les 2p fonctions 


p1(s), ya(s), eae, p(s), 
os nae Ualt), 1 Uplt) 
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fonctions principales relatives a 4. 

Si un zéro de degré p de D(A) est considéré comme p zéros distincts, 
on peut dire qu’une paire de fonctions principales correspond a chaque zéro 
de D(A). 

Les fonctions principales relatives & A; satisfont aux équations (37). 
Elles constituent done un systéme biorthogonal. 

Si maintenant, on écrit l’une des équations (32) 


b b 
[ vteneurtyar =f rls.7) galrt) dr = als.) 


on voit que 2 est solution d’une équation intégrale ott le noyau ¢1(s,t) 
est de la forme (38) étudiée au n° 6, 2°, p. 45. D’oti l’on déduit que ye est 
de la forme 


yo(s,t) = S> ain vi(s) va(t) 


ot les a;z sont certaines constantes. On verra de méme que 


y3(s, t), nghely pr(s,t) 


sont d’une forme analogue. Dans un travail récent(!), M. Lalesco a 
montré comment en effectuant une substitution biorthogonale convenable 
sur (39), on peut réduire la résolvante 4 une forme canonique. 

Il y a lieu de ne pas confondre le systeme des fonctions principales (39), 
avec celui des fonctions fondamentales définies selon M. Schmidt au n° 43, 
p. 110; non plus qu’avec le systéme des 2q solutions des équations intégrales 
homogénes associées. Toutefois, on démontre facilement que ces derniéres 
fonctions sont des combinaisons linéaires des 2p fonctions principales (on a 
vu au n® 25, p. 81, que q < p). Dans le cas ot le pole est simple (r = 1) on 
démontre méme que ces deux systémes de fonctions coincident. Ce résultat 
est, en tout cas, évident pour nous quand le noyau est symétrique. 


@) Bulletin de la Soc. Math. de France, t. 39, 1911, p. 95. 
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